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significant contributions fo the 
art and science of rocketry. 


By developing new and better 

propellants (both solid and lig- 
vid)—by designing and build- 
ing improved power plants to 
utilize these fuels —by furnish- 
ing essential support equipment 
... Thiokol helps to strengthen 
the nation’s defenses, helps 
push back Qur spatial frontiers. 
Engineers, Scientists: perhaps there’s 
a place for you in Thiokol’$ @xpand- 
ing organization. Our new projects 


present challenging problems and 
a chance for greater responsibility. 
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i: ...For over twenty years—from the 


at Annapolis to the launching pads at Cape Canaveral— 
setting the standards for the control of high pressure fluids 
...the regulators which served as the critical control on 


World War II flame throwers and torpedoes...and made 


missile development work possible...Grove High Pres- 


sure Regulators...frequently imitated but never equaled. 


DROVE VALVE and REGULATOR COMPANY 
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This Journal is devoted to the advance- 
ment of astronautics through the dissemina- 
tion of original papers disclosing new scientific 
knowledge and basic applications of such 
knowledge. The sciences of astronautics are 
understood here to embrace selected aspects 
of jet and rocket propulsion, space flight 
mechanics, high-speed aerodynamics, flignt 
guidance, space communications, atmospheric 
and outer space physics, materials and struc- 
tures, human ee, overall system 

is, and possibly certain other scientific 
— The selection of papers to be prin ted 
overned by the pertinence of the topic 

to pag pe for of astronautics, by the current or 
probable future significance of the research, 
and by the importance of distributing the in- 
formation to the members of the Society and 


to the profession at large. a 


Manuscripts must be as brief as the proper 
presentation of the ideas will allow. [x- 
clusion of dispensable material and con: se- 
ness of expression will influence the Editors’ 
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Air Force 
‘““Sunday 
Punch” 


Boosted into space by the fiery thrust of three commentary, for all the world to heed, of the ne- 
fa huge rocket engines, the seven-story Atlas inter- cessity to maintain the peace. RCA’s Missile and 
continental ballistic missile roars upward from Surface Radar Department has been privileged 
its Cape Canaveral launching pad. Quickly it to design and develop ground check-out, launch 
sheds the frost encrusting the liquid oxygen control and cabling equipment as a major sub- 
tank and races to its predetermined destination contractor to Convair (Astronautics) Division 
in the far reaches of the globe. In its size and of General Dynamics Corporation, the Atlas 
7 range and capability, the Air Force es isa prime weapons systems contractor. 


RADIO CORPORA TION of AMERICA 


Taide) @ DEFENSE ELECTRONIC PRODUCTS 


CAMDEN, N. J. 
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The U.S. Army’s new Q-5 Kingfisher 
was designed by Lockheed’s Missile Sys- 
tems Division to provide our mighty 
arsenal of ground-to-air missiles with a 
realistic test of marksmanship—against 
high-altitude targets moving at super- 
sonic speeds over 1500 miles-per-hour. 
The Kingfisher is 38-feet long, 2( 
inches in diameter, has a 10-foot wing 
span and weighs more than 7600 pounds. 
As it flashes across the skies it electroni- 
cally simulates any desired size and type 
of “enemy” plane or air-breathing missile. 
The Kingfisher’s electronic Firing 
Error Indicator instantly and accurately 
tells ground controllers whether missiles 
fired at it are “hits” or “misses’”—and 
automatically evaluates each missile 
angle-of-attack, miss-distance, and othe 
- highly important technical data. 
Undamaged by “hits” scored on it 
electronic image, the Q-5 Kingfisher is 
parachute recovered after each flight. 
This Lockheed-developed “‘missile 
with 9 lives” will enable the U.S. Arm 
to achieve hitherto impossible proficienc 
in missile marksmanship against super 
sonic targets—at a saving to taxpayers of 
approximately half a million dollars on 


each recovery flight. - 


Q-5 is dropped by plane at 35,000 Q-5 is detected as ‘“‘enemy”’ by Missiles fired at Q-5 are like those 
feet ((A) in diagram). Then its twin ground radar (B), and its speed, alti- used in wartime—but lack high- 
rockets ignite, propel it to speeds re- tude, and course are fed into fire- explosive warheads. Nike missile (Q 
quired to operate its ramjet engine. control computer of Nike battery. scores “hit” on Q-5’s electronic image, | 
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Above: Entering oval flight pattern, Q-5 attains speeds over 1500 
mph. Second ground radar (D) and missile-launching battery (E) practice 
their marksmanship until Q-5 Kingfisher’s fuel supply is exhausted. 
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Left: Landing on its nose-spike in a remote, uninhabited area, after 
floating down by parachute (F), the Q-5 is recovered by U.S. Army 
ground crews—to be refueled and refitted for future flights. 
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Soon space probes will relay 
back to the earth pictures of 


the other side of the moon. 


< — * This scientific representation based on current knowledge was pre- 
pared under the supervision of Dr. I. M. Levitt, Director of the 
Franklin Institute Planetarium. 


At this time no one knows. But intricate electronic devices in | 
projected lunar vehicles will reveal this hidden surface. Instru- | 
mentation has extended the long arm of man to reach as far as 
the mind can project. With such devices as a key, science can a 
unlock the door to the future and to the very universe itself. 


At the Decker Corporation our sole occupation is instruments 
— instruments which range from a device to measure a millionth 
of an inch on earth to one recording the density of the most 
tenuous of the space atmospheres subject to man’s reach. 


= 
CORPORATION Bala Cynwyd, Pa. i 
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Kenefits of Fluorine in Rocketry 


LARGE rocket is largely a flying chemical tank or a 

group of such tanks. The release of the potential 
energy of chemical reaction residing in these tanks is controlled 
by metered flow through a finely engineered pumping system 
or pressure drive into the liquid rocket engine, or by the 
chemically controlled combustion rate of a solid propellant 
grain. 

The takeoff weight of a large rocket is predominantly that 
of the chemical propellant; the cost of a rocket, however, is 
mostly in its structure, in its operating parts, and in its guid- 
ance and payload (1).!_ An improvement in the energy storage 
value of the propellant system, therefore, reflects in magnified 
manner on the total cost of the device (2). 

Among all stable chemical materials, liquid fluorine oxidizer 
contributes the highest possible energy storage value, whether 
measured in terms of volume or weight (3). This is an out- 
standing reason for interest in fluorine chemistry. Such is the 
versatility of this chemistry, though, that fluorine chemicals 
other than liquid fluorine can substantially advance the art of 
high energy and storable propellants; furthermore, fluorine 
chemistry has given many compositions besides propellants, 
useful in the unusual chemical and rigorous physical environ- 
ments encountered in various parts of the rocket. 

Fortunately, fluorine chemistry has for some time been of 
interest in other major areas of technical advancement, in- 
cluding atomic energy generation, in steel making, oil re- 
fining, and in protection from corrosive or high temperature 
environments similar to those encountered in high per- 
formance rocketry. The advances which have been made in 
the last 10 or 15 years in fluorine chemistry have been nicely 
timed with the requirements of rocket development. 


Sources of Fluorine 


Requirements for chemicals used as rocket propellants can 
be very great, and therefore, since liquid fluorine and some 
fluorine compounds are being considered as propellant com- 
ponents, it is necessary to examine the present state of supply 
of the raw materials and intermediate products. 

As recently as 1952, a well-known report (4) stated that the 
United States’ known reserves of fluorspar would be sufficient 
for only some 15 years at the 1950 rate of consumption, and 
that reserves in other free countries would not be adequate to 


Received Dec. 15, 1958. 
* Numbers in parentheses indicate References at end of paper. 
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7 Pennsalt Chemicals Corp. 
Philadelphia, Pa. 


The author is Assistant Manager of Research and Development at Pennsalt Chemicals Corp., 
Since joining Pennsalt in 1938, he has carried out and directed chemical research in 
fluorine chemistry, electrochemistry, thermodynamics and the diverse special fields of interest of 
Dr. Gall was born in 1913 in Chicago and attended Northwestern University (B.S., 
1934) and the University of Chicago (Ph.D., 1939). 
American Chemical Society (Chairman, Division of Inorganic Chemistry, 1958) and the Electre- ' 
chemical Society (Chairman, Philadelphia Section, 1952-3). . 


He is a member of the American Rocket Society, 


supply U. S. deficiencies during the next 25 years. This con- 
clusion was based on an estimate from the Bureau of Mines on 
U.S. fluorspar reserves as of January 1944 (5). This estimate 
stated U.S. reserves of crude fluorspar to be 14,724,000 tons. 
Twelve years later, however, after the production of some 
8,500,000 tons of crude ore from U. S. mines, a new estimate, 
dated November 1956, set remaining reserves at 22,500,000 
tons crude ore (6). Even this large increase in estimated 
reserves remaining after 12 years of mining is considered by 
many to be conservative. Thus the history of fluorspar is 
similiar to that of other resources suspected to be in danger of 
exhaustion, for which new discoveries and better recovery 
methods have annually increased the estimated supply even 
while withdrawals are in progress. 

The outlook for fluorspar is further improved by large de- 
posits in Mexico and by rapidly rising imports from that coun- 
try and others. 

In this discussion, emphasis has been placed on fluorspar 
because it is by far the largest source of fluorine, and because 
so-called acid grade fluorspar is very much the largest source 
of fluorine chemicals. Table 1 shows total consumption of 
fluorspar both domestic and foreign in recent years in the 
United States, the portion of this distributed as acid grade 
spar, and the estimated production of hydrofluoric acid, which 
is the primary intermediate for nearly all fluorine chemicals 
(one ton of anhydrous hydrofluoric acid requires 2.5 tons of 
acid grade spar). 


Table 1 United States consumption of fluorspar and 


production of hydrofluoric acid, tons 


- Total fluorspar Acid grade spar Hydrofluoric acid 


(7) (7) 
1955 570,261 248,218 99,500 
1956 621,354 289 523 115,500 
1957 644,688 328,672 131,000 


Table 2 gives a perspective on the demands for fluorspar by 
its various uses; steel making and hydrofluoric acid manu- 
facture divide over 90 per cent of the spar. Anhydrous hy- 


drofluorie acid, in turn, is used about equally in chemical 
production and all other requirements. 

Anhydrous hydrogen fluoride is the essential chemical for 
the manufacture of fluorine products of interest to the missile 
program. 


1 and 2 will be com- 
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Table 2 Measures of fluorine products uses, 1956 

Fluorspar (8) Hydrofluorie acid (9) 

Use group % of total Use group % of total 

consumption consumption 

steel making 45.6 chemicals 49 
other metallurgy 0.4 aluminum 39 
glass and enamel 5.8 steel Ne 4 
hydrofluoric acid 46.6 petroleum 4 
all other 1.6 glass 2 
all other = 
100.0 100 


pared with the added burden which would be placed on the 
fluorine chemicals industry by the widespread use of fluorine 


propellants. The conclusion is that the burden is relatively 
small. 

Chemistry of Fluorine 


Recent advances in the basic chemistry of fluorine have re- 
vealed a system based on the following facts: 

The fluorine atom is small. This is shown in Table 3, where 
the ‘‘co-valent radius” measures the effective size of an atom 
in compounds where the atom exists in a nearly uncharged 
condition. The volume of a fluorine atom is only one-fourth 
that of a chlorine atom, and less than one-tenth that of an 
atom of iodine (27). 

The small size makes possible the packing of many fluorine 
atoms around a single atom of another element, so that many 
compounds rich in fluorine can exist. Examples include SF,, 
UF,, IF;, and the ions HF,~, BrF,~, SiFs=, TaF;~ and others. 

The fluoride ion is small. Table 3 shows that fluorine shares 
with oxygen the distinction of forming the smallest negative 
ion by transfer of an electron from an atom of another element. 
Fluoride and oxide ions are for this reason “hard” ions, bond- 
ing to other kinds of atoms with little deformation from elec- 
trostatic forces and thereby commonly retaining most of their 
ionic character. 

One consequence is that metal fluorides, in comparison with 
chlorides, bromides and iodides, melt and vaporize at higher 
temperatures and are usually white (15). 

A consequence of the similarity between fluorine and oxy- 
gen is ready replacement of one by the other. Thus fluorides 
lower the melting point of oxides by dissolving in them. More 
importantly, by replacing doubly-charged oxide with singly- 
charged fluoride, fluorides disrupt large polymeric molecules 
and so lower the viscosity of a melt, causing slags to flow 
away freely or fluxes to spread over and protect metal sur- 
faces (16). 

In these fluxing actions, so important in the manufacture of 
steel and aluminum and in brazing and welding: fluorine per- 


Fluorine will commonly substitute for oxygen or hydroxy] 
in organic compounds with little change in the major chemical 
properties of a compound; thus the replacement of oxygen by 
fluorine in minerals is common (17), and the interesting new 
compound, perchloryl fluoride ClO3F, possesses the stability 
and inertness of the perchlorate ion ClO,~ (18). 

Fluorine is the most electronegative element. Loosely defined 
as the vigor with which an element combines with the “posi- 
tive’ elements like hydrogen, sodium or magnesium, the con- 
cept of electronegativity is now subject to a variety of more 
exact, but not uniformly consistent, measures, such as con- 
parative ionization potential (reluctance of an atom to lose an 
electron), electron affinity (willingness to gain one), electrod» 
potential (decrease in free energy going from the element to its 
ion in solution), and many others (19). Table 3 lists the fre- 
quently-quoted electronegativity values derived by Pauling 
from the observed excess energy of a bond over that com- 
puted from an arithmetic combination of the energies of the 
diatomic molecules of the elements themselves (20). Almost 
regardless of which scale of electronegativity is used, fluorine 
stands above all of the other elements (21). 

This is another reason for the largely ionic nature of meta! 


- fluorides, and is further exemplified in the high heat of forma 


tion of many fluorine compounds. In addition, the hig): 
strength of fluorine-substituted organic acids as trifluoroaceti: 
acid, or of trifluoromethanesulfonic acid (one of the strongest 
known acids), results from this extreme electronegativity (22). 

Fluorine can also, being more electronegative than oxygen 
form the stable compound oxygen difluoride OF2, in whic! 
oxygen is the positive element. This is in sharp contrast with 
the other halogens, chlorine, for example, forming the un 
stable monoxide C1,0, in which oxygen is negative. Similarly 
nitrogen forms stable NF3, bearing positive nitrogen, but 
highly explosive NCl;, with negative nitrogen and _ positive 
chlorine (23). 

Since fluorine must always be negative relative to oxygen, 
oxyacids analogous to those of the other halogens cannot be 
expected (24). Almost certainly, perfluoric acid HFO,, or 
fluoric acid HFOs, or hypofluorous acid HOF, will not there- 
fore be prepared. 

Fluorine in its compounds is always univalent and negative. 
One or more characteristic numbers (valences) assigned to 
‘ach element form a consistent set indicating the number of 
bonds to other atoms which can be formed in stable com- 
pounds of the element. Fluorine, with only apparent and 
explainable exceptions, is thus limited to a single bond in its 
connections with other elements. 

This can be concluded from the electronic structure of the 
fluorine atom which would require excitation well beyond 
usual chemical energies to form bonds by donation of elec- 
trons. This is not entirely the result of high electronegativity, 
but occurs also because the energy difference between the high- 
est filled energy level in the normal atom and the next available 
level is large. With the other halogens, available energ\ 
states exist close above the ground state, so that these ele- 


forms its title role (Lat. fluere, to flow). — ments can show a series of positive valences (25). - 
Table 3 Some characteristics of the elements neighboring fluorine 
Property F Cl Br I O Ss N ig Reference 
Co-valent radius, A 0.72 0.99 1.14 1.33 0.74 1.04 0.74 1.10 (10) 
Ionic radius, A 1.36 1.81 1.95 2.16 1.40 1.84 i) ae (11) 
Relative electronegativity 4.0 3.0 2.8 2.5 3.5 2.5 3.0 2.1 (12) 
Heat of dissociation of the dia- 36.6 58.0 46.1 36.1 118.: (76.6) 225.1 (13) 
tomic molecule, gas at (14) for N 
25 deg C kcal/mole 
Heat of formation of the simplest —64.2 —22.1 —8.7 6.2 —28.9 —2.4 —3.7 0.7 (13) | 


hydrogen compound, gas at 25 
deg C, per bond kcal/mole 
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In consequence, fluorine is never a central element in a com- 
pound but must always be in terminal or exterior position, 
often, therefore, giving a strong dipolar character or negative 
external shielding to a molecular group. 

The heat of dissociation of the fluorine molecule is small. 
Some 25 estimates of the heat of formation of the fluorine 
molecule F, have been made (26). Because of the rigorous ex- 
perimental problems found in attempts at direct measure- 
ment, this important quantity has been much in doubt until 
recently. The temptation to assume a normal monotonic 
increase in heat of formation of the halogens from iodine 
through fluorine (see Table 3) has finally been resisted in the 
light of many good arguments, and a value close to 37.7 keal 
per mole, well below the 58.0 kcal per mole established for 
chlorine, must now be accepted (135). 

This somewhat surprising result has consequence mostly 
in our understanding of the high energy release observed in 
reaction of fluorine with fuels, for little heat is required to 
break up the fluorine molecule during the reaction. There is 
also significance in understanding the mechanism of reactions 
nvolving fluorine molecules. 

However, since calculations of energy release expected in 
fuorine-maintained combustion and hence of the performance 
to be expected from a liquid fluorine engine, are made from 
reliable measurements of reactions related to the combustion 
process, a change in the recorded value of heat of formation of 
fluorine molecule does not affect estimates of specific impulse 
or other rocket performance parameters. 

The low energy required to dissociate fluorine, together 
with the usually large heats of reaction, may account for the 
ready reaction, usually leading to immediate ignition (hy- 
pergolicity), of fluorine with most fuels. 

The heat of formation of bonds to fluorine is often large. For 
example, Table 3 gives the heat evolved when one gram 
equivalent weight of each element listed unites with one of 
hydrogen, and the large energy release with fluorine is note- 
worthy. These figures are heats of formation from the stable 
molecules of the elements, for example Fs, and result in part 
from, but do not separately measure, the low heat of forma- 
tion of fluorine. However, they directly demonstrate the very 
large energy release to be expected when fluorine combines 
with fuels, and from this the excellent engine performance 
which becomes possible. 

Conversely, the high heat of formation of bonds with 
fluorine indicates that many fluorine compounds will them- 
selves have high stability. In the rocket engine exhaust for 
instance, HF is a product existing as the intact molecule even 
at very high temperatures; it therefore carries to waste 
little energy except that residing in its high temperature. In 
fact, next to nitrogen, HF is the most stable diatomic molecule. 

Similarly, the compounds of fluorine with sulfur as SF, or 
with carbon as CF, and the fluorocarbons and fluorocarbon 
plastics, are highly stable to elevated temperatures. 

A shell of fluorine atoms around a central atom often shields 
the latter. Many highly fluorinated compounds, such as SF%, 
CF, or fluorocarbons, have low heats of vaporization. The 
small hard negatively-charged fluorine atoms peppering the 


surface of such molecules provide mutual repulsion between | 


the molecules and thus ease of separation as the substance 
goes from liquid to gas. 

Some of the noteworthy inertness of the compounds cited | 
must also result from this characteristic. Fluorocarbons are 
not thermodynamically stable towards water (27), the ulti- 
mate equilibrium situation being carbon dioxide and hydrogen 
fluoride. The extreme unreactivity of fluorocarbons toward 
water may be due as much to the remote chance for a basic 


water molecule to penetrate the picket fence of negatively- 


charged fluorine atoms as to the resistance of fluorine to leave 


the carbon after initial attachment of water oxygen. A defini- — 
tive explanation of the hydrolytic stability of fluorocarbons — 


has not yet been published. 


Hydrogen bonding is frequent in fluorine chemistry. The hy-— 


drogen atom may sometimes be attracted simultaneously to 
highly electronegative atoms, and thereby acts as a bridge 
between them. 
rather weak (e.g., 6.7 kcal/mole of FHF) they commonly 
have marked effects on the properties of compounds (28). 
Thus the boiling point of hydrogen fluoride is abnormally 
high because of extensive formation of groups of HF molecules 
in chains of various lengths, thus: F—H- -F—H- -F—H.- - 
(29). This effect persists into the vapor and causes large 
deviations from normal gas behavior. In fact, HF is perhaps 
the most imperfect gas thus far studied (30). 

Many other examples occur, such as the apparent weakness 
of HF as an acid in water solution (though its behavior in con- 
centrated solutions or in some catalyst activities show HF to be 
one of the strongest acids (31)), the stability of acid fluorides, 
such as K*(FHF)~ (32), and the excess crystal energy of 
ammonium bifluoride. 

Fluorine has only one stable isotope. This simplifies the 
nuclear chemistry of fluorine and renders volatile fluorides 
handy for isotope separation, hence, for example, the useful- 
ness of uranium hexafluoride (33). 

Fluorine has no long-lived radioactwe isotopes. This makes 
radioactive tracer studies with fluorine awkward. The long- 
est-lived isotope has a half-life of only 112 min (34). 


Recent Advances Pertinent to Rocketry 


In the years 1950 through 1957, about 3000 papers on 
fluorine chemistry or technology appeared in the technical 
literature. Out of this vast research activity the following 
areas of discovery are selected for the attention of workers 
in jet propulsion. 


Halogen Fluorides and Complex Fluorides 


Since 1931, when bromine pentafluoride was first prepared 
(35), no new halogen fluorides have been found; and one of 
those previously accepted, bromine monofluoride, can prob- 
ably not be isolated (36). It would be an achievement to 
prepare, for example, CIF; or BrF;; however, stereochemistry 
suggests that these will not be made. Table 4 lists the known 


members of the halogen fluoride series. vera 


CIF; CIF 
— 103.8" 
1.62 (—104 deg)? 
—13.3 (g)? 


12.02 
1.825 (20 deg)* 
—37.0(g)? 


boiling point 

density at T deg C 

heat. of formation, 
keal/mole, 25 C 


* Highly uncertain. 


Table 4 Properties of the known halogen fluorides — 


40. 76/ 
2.465(25 deg)° 


@ (37), °(13), © (39), (43), (38), 7(40), (41), (42), (44), (132), *(134). 


—122(at 0 deg —75* — 225.1 (g)° —196.0 


BrF; IF; 


5(subl.)* 100.54 


127.6" 
2.8(6 deg)? 3.19(25 deg)? 


2.843(8.8 deg)” 
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Though the bonds, even with fluorine, are | 
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In 1949, it was discovered that BrF; and IF;, but not the 
other halogen fluorides, have a high electrical conductivity 
and form a series of compounds indicating that the liquid halo- 
gen fluorides are appreciably self-ionized as (45) 


Thus when, for example, KF is dissolved in these materials 
and the excess halogen fluoride is evaporated, the new com- 
pounds KBrF, and KIF, remain. These compounds may be 
thought of as bases, in the same way that KOH is a base in 
water (133). 

Similarly, if a compound of an acid-forming element like 
antimony or tin reacts with BrF3;, the corresponding acid, for 
example BrF,-SbF¢, can be made. 

If a solution of KBrF, in BrF; is added to a solution of 
BrF;SbF, in BrF;, the expected acid-base reaction takes 
place; BrF; is regenerated, and the salt KSbF, forms. 

The ionizable halogen fluorides thus provide an engaging 
example of acid-base chemistry, but meie importantly they 
provide a medium for the preparation of a number of new in- 
organic fluorides (46). 

In another interesting preparative technique, elemental 
fluorine is reacted with a mixture of metal chlorides; the re- 
sult is a complex fluoride with the metal sometimes apparently 
in an unusually high oxidation state. Thusa mixture of KCL 
and MnCl, yielded K:MnF, containing 4-valent manganese, 
potassium and copper chlorides K;CuF., and potassium and 
nickel chlorides K,NiF (47). 

Another powerful method for preparing new inorganic 
fluorides employs a hypothetical intermediate, fluosulfinic 
acid HFSO.. This is accomplished by reactions of fluorides 
in liquid sulfur dioxide. The simple salts like KSO,F have 
been approximated (48). 

Out of these new ways of working, and through the now 
ready availability of elemental fluorine, many new inorganic 
fluorides are on hand. 


Perchloryl Fluoride and Other Fluorides of the Oxy Acids 


The principle of replacement, whereby the fluoride ion re- 
places the similar oxide ion, has been applied to the negative 
ions of the oxygen acids, such as sulfuric, perchloric and ni- 
tric. The resulting fluorine oxygen compounds are one less 
in ion charge, so that a singly negative anion becomes a neutral 
molecule, and the compound or ions so obtained often possess 
high chemical inertness and thermal stability while retaining 
the strong oxidizing value of the oxygen acids. These com- 
pounds can also be regarded as the acid fluorides of the oxygen 
acids wherein one or more OH groups from the parent acid 
are replaced by fluorine. 

An outstanding example of this type of compound is per- 
chloryl fluoride (49), the acid fluoride of perchloric acid 
HCI10,. Perchlory] fluoride may also be thought of as derived 
from the perchlorate ion ClO,—', in which replacement of one 
oxygen atom by one fluorine atom has given a neutral com- 
pound ClO;F. Like the perchlorate ion, perchlory] fluoride is 
stable and chemically rather inert. However, it possesses 
high oxidizing power contributed by the heptavalent chlorine, 
the three oxygen atoms and the fluorine atom. This oxidizing 
value appears at elevated temperatures, and perchlory! fluo- 
ride is therefore an excellent oxidizing substance to be em- 
ployed with hydrogen or carbon containing fuels. The prop- 
erties of this new liquid storable oxidizer are discussed more 
quantitatively in a later section. 

As by-products of perchloryl fluoride chemistry, two new 
classes of compounds have been discovered: One series is ob- 
tained through the strongly acidic perchlorylamide ClO;NH: 
which forms salts containing CIO;N™ ions. This, incidentally, 
is another interesting example of ion replacement chemistry; 
the nitrogen atom, though not replacing oxygen as smoothly 
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as does fluorine, here forms a reasonably stable species by re- 
placing oxygen in the ClO,~ ion (50). 

Another series of new compounds arises when perchlory| 
fluoride reacts in the presence of aluminum chloride with an 
aromatic ring, whereby the ClO; radical, new in organic 
chemistry, enters to form compounds like perchlorylbenzene 
and its derivatives. The perchlory! group is similar to the 
nitro group, but contributes twice as much to the oxygen 
balance of the compound. The perchlory! aromatics can be 
exploded and need to be investigated as components of ex- 
plosives or of solid propellants (51). 

Perchlory] fluoride is aiso an interesting selective fluorinat- 
ing reagent, whereby reactive hydrogen, as in compounds con- 
taining active methylene, such as C,H;OOCCH;COOC2H;, 
diethyl] malonate, are converted to the otherwise difficultly- 
prepared fluoro compounds (52). 

The principle of replacement by fluorine of oxygen acids 
and ions results in the existence of a number of other interest- 
ing fluorine compounds, such as ClO.F, I0.F, IO;F, SOF», 
8.0.F2, NO.F, Mn0O,F, CrO.F», ete. (53, 54, 55, 56). 


Hypofluorites 


Unlike chlorine, and for the reasons indicated previously, 
fluorine cannot exist in stable compounds except in a co-valent 
linkage or in a negative oxidation state. Therefore compounds 
analogous to the hypochlorites, as the familiar sodium hypo- 
chlorite NaOCl! in bleaching solution, cannot exist. Yet a 
number of compounds containing the OF group have recently 
been produced. There is still no reason to consider these as 
containing positive fluorine: Since fluorine is more electro- 
negative than oxygen, the bonding must leave oxygen the 
positive and fluorine the negative partner. 

Examples of hypofluorites include the highly explosive 
fluorine nitrate FONO, (57) and fluorine perchlorate FOCIOs; 
(58). We note that in structure as well as in properties, these 
compounds are radically different from those in which fluorine 
is directly connected to chlorine or to nitrogen, as in the 
highly stable perchlory] fluoride FCIOs, or nitry] fluoride FNO.. 
Another example is trifluoromethyl hypofluorite CF;OF 
which reversibly decomposes at elevated temperatures to car- 
bony! fluoride and fluorine, and reacts with carbonyl fluoride 
to form perfluorodimethy] peroxide, CF;OOCF; (136). 

Hypofluorites have not yet been controlled for high energy 
use, and this may be a worthwhile project. 


N-F Chemistry 


Much interest attaches to possible compounds of nitrogen 
and fluorine. In addition to the long known nitrogen tri- 
fluoride NF;, there is now difluorodiazine N2F, (59), and the 
very recently discovered tetrafluorohydrazine N2F, (60). 


Synthetic Fluoro-Minerals 


Crystal chemists unraveling the wide variety of lattice 
structures existing in natural minerals have long recognized 
the principle of fluorine replacing oxygen, but only recently 
have there been efforts to produce new synthetic crystal lat- 
tices in which this replacement is made. New structures can 
be visualized simulating the typical silicate lattices, which in 
silicate minerals are typically made up of tetrahedra produced 
by four oxygen atoms surrounding the silicon central atom, 
each tetrahedron being linked by one to four oxygen atoms to 
neighboring tetrahedra, and forming minerals which tend to 
be fibrous, plate-like or isotropic, respectively (61). Fluoride 
replacement also modifies the ionic charge of the array of the 
silicon-oxygen-fluorine tetrahedra and permits balancing with 
a variety of metallic ions. 

This science seems to be in an early stage of growth. A 
significant practical benefit has come in the synthesis of 
fluorophlogopite, a mineral closely resembling natural mica, 
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with improvement in high temperature properties resulting 
from reduction of anionic charge and elimination thereby of 
virtual water (62). This is the most promising synthetic 
mica now being investigated in large-scale production, which 
may reduce our dependence upon a somewhat uncertain 
natural resource. 


Fluorocarbon Chemistry 


The United States atomic energy program gave fluorine 
two strong impetuses. One came from the requirement for 
uranium hexafluoride as a volatile compound of uranium 
suitable for the separation of isotopes by gaseous diffusion; 
ihe other was the requirement for stable liquids and plastic 
inaterials to resist the violently reactive uranium hexafluoride 
und the element fluorine which was employed in its prepara- 
tion. Fluorocarbons, analogous to familiar hydrocarbons 
hut in which all or most all of the hydrogen has been replaced 
by fluorine, meet the latter need. The first methods to pro- 
duce fluorocarbons were cumbersome and yields were poor. 
Chus the direct reaction of fluorine with a hydrocarbon or with 
a chlorinated hydrocarbon or with a hydrocarbon partially 
fluorinated by less vigorous methods, completes the replace- 
ment of hydrogen by fluorine, but the energy evolved in the 
replacement is sufficient to damage the basic carbon skeleton 
(63). The method of direct fluorination, however, is still 
being studied and advances are being made (64). 

Better results were obtained by the use of a fluorine carrier, 
such as cobalt trifluoride, which in turn has to be made by the 
use of elemental fluorine (65). A rational procedure is poly- 
merization of short unsaturated fluorocarbons or chloro- 
fluorocarbons; out of this procedure came the liquid and solid 
polymers based on chlorotrifluoroethylene, CFCI==CF, or 
tetrafluoroethylene CF:—=CF, (66). 

Recently, a technique capable of producing fluorocarbons 
and a great variety of highly fluorinated organic and semi- 
inorganic compounds has been devised. This was developed 
by recognizing that a fluorocarbon iodide, in the simplest 
case CF;I, will at elevated temperatures or by activating 
short wave-length light, produce iodine atoms and a fluoro- 
carbon free radical. In the example given, CF;I yields CF; 
and I. The free radical fluorocarbon is now available for 
many reactions. One example is addition to an unsaturated 
molecule, and CF;I, for example, will add to tetrafluoro- 
ethylene F,C—=CF; to give CF;CF.CF.I, a new fluorocarbon io- 
dide. This in turn can be the starting point for a new fluoro- 
carbon free radical (67). Since the addition of a fluorocarbon 
free radical and the iodine atom to a double bond can take 
place many times in a single reaction sequence, a kind of 
polymerization results, which, since addition is to opposite 
ends of the molecule, is called telomerization. Many ingeni- 
ous variations of this technique have been carried out, giving 
a new method of tailor-making high molecular weight fluorine 
compounds (68). 


Organometailic Compounds 

The use of the fluorocarbon free radical technique noted in 
the preceding section has been valuable in the development 
of new chemistry in the semiorganic fields. For example, 
when sulfur is reacted with trifluoromethy] iodide, the com- 
pound CF;S-SCF; is obtained. This can be converted by the 
action of short wave-length light to the monosulfide CF3SCFs, 
which can be converted with chlorine to CF;SCI, and by oxi- 
dation and hydrolysis to CF;SO3H. This last compound, tri- 
fluoromethanesulfonic acid, is the strongest known carbon con- 
taining acid, and in some solvents it behaves as though it 
might be the strongest of all acids (22). 

Similar chemistry has led to a variety of trifluoroalkyl 
metallic compounds including the perfluoroalkyl magnesium 
halides analogous to the traditional Grignard reagents so use- 
ful in synthetic organic chemistry (69). New perfluoroalky] 
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derivatives of arsenic (70), nitrogen (71), phosphorous (72), 


mercury (73), antimony (74) and silicon (75) and other ele- | 
ments have been made by the perfluoroalky] free radical | 


technique. 


This powerful synthetic method, however, is not always | 
highly 
fluorinated derivatives are also obtained by electrochemical 
In this way, for example, the above mentioned © 


most convenient, and in the organometallic field, 


 fluorination. 


ene’ compounds have been duplicated, and the 


derivatives of sulfur hexafluoride CF;SF; and CF.(SF;)2 have © 
been made by electrofluorination of carbon disulfide (76). p 
The highly fluorinated sulfur compounds are of interest “fl ; 


alternates to sulfur hexafluoride in gaseous insulators for high | 
voltage equipment (77). — 


Highly Fluorinated Aromatic Compounds 


New interest has been raised in perfluoroaromatics, ex-_ 


emplified by hexafluorobenzene, through discovery of better 
methods of preparation. Thus simple pyrolysis of tribromo- 
fluoromethane CFBr; gives a quite useful vield of CeFs. 
The compound, though less reactive than the corresponding 
chlorine compound, should react with alkalies to give penta- 
fluorophenol. The latter compound will be rather strongly 
acid and can be the starting point for a new series of highly 
fluorinated organic compounds (78). 


Advances in Fluorine Chemical Technology 


This section describes materials recently available out of 
fluorine chemistry of especial interest to rocket engineering. 


High Energy Compounds 


Table 5 shows the pertinent characteristics of fluorine and 
fluorine compounds useful as liquid propellant oxidizers. 

Liquid fluorine participates in bipropellant systems of the 
highest energy storage value. Like oxygen, fluorine is a 
cryogenic oxidizer, which must be condensed and stored at 
liquid air temperature. It is, however, violently incendiary, 
igniting most organic materials, water and silicates, including 
concrete and asbestos (87). Fortunately, the ignition tem- 
perature of bulk metals in fluorine is generally high; this 
appears to be due to a protective self-healing film of metal 
fluoride. Common metals have not ignited under severe con- 
ditions in liquid fluorine (88), and in gaseous fluorine the 
ignition temperatures appear to be at least several hundred 
degrees C (89), so that accidental ignition is not likely except 
in contaminated situations. 

Although the problems of handling liquid fluorine are 
severe, the obstacles are one by one being eliminated, not so 
much by the use of unusual materials, but by recognition of 
the importance of great cleanliness, avoidance of conditions 
of hydrostatic shock, and test area design to reduce damage 
resulting from a leak or a spill (90). 

The advantages of fluorine are great. Thus, an immediate 
advance in rocket engine performance could be obtained 
through the use of fluorine in conjunction with a conventional 
fuel, such as hydrazine, when as is noted in Table 5 an 
increase in specific impulse of 46 see would be obtained over a 
conventional oxygen-hydrocarbon system. Because of the 
higher densities of fluorine and of hy drazine fuel, the compari- 
son in density impulse would show a further advantage. 

Much higher values of specific impulse are calculated for 
liquid fluorine using liquid hydrogen fuel. The very low 
density of liquid hydrogen seriously reduces the density im- 
pulse of systems using it, i.e., fuel tank volume becomes large. 
Fortunately the weight ratio of fluorine to hydrogen is large, 
and the density of fluorine itself is high, so that mean propel- 
lant densities of practical magnitude are realized. In this 
respect fluorine is considerably better than oxygen as an ox- 
idizer for hydrogen, as Table 6 shows. 
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Table 5 Fluorine containing liquid propellant oxidizers compared with oxygen 
Performance* 
Oxidizer Formula Melting Boiling Density Heat of formation Weight ratio, Specific 
ow point, C point, C (T deg C) 25 C keal/mole gas oxidizer: impulse 
bree: Fuel fuel sec 
fluorine F, —217.96 —187.92 15.4(—196 deg) 0 hydrazine 2.08 318/ 
hydrogen 4.49 3819 
perchlory] fluoride ClO;,F —146 —46.8  1.392(25 deg) —5.12 UDMH 2.55 274" 
hydrazine 1.1 274" 
chlorine trifluoride CIF; —83 12.0 1.825(20 deg) —37.0" hydrazine ae | 268” 
bromine penta- BrF; —62.5 40.76 2.465(25 deg) — 122° hydrazine 2.19 225" | 
fluoride (at 0 deg K) | 
nitrogen trifluoride NF; —206.65 —128.87 1.54(—128 deg)¢ hydrazine 2.64 276**" 
oxygen difluoride OF, — 223.8 —146.5  1.65( —190 deg) 7.6 UDMH 1.8 317! 
oxygen Oo —218.76 —182.97 1.14(—183 deg) 0 JP-4 2.3 272/ 
* Specific impulse calculated for frozen equilibrium at 600 psia chamber pressure, 14.7 psia exhaust pressure; propellants injected 
as liquid at 25 C, except F., NF;, OF. and O, at their boiling points. 
** Value for heat of formation of NF; should probably be more negative, reducing calculated specific impulse (85). 7 a | 
@(79), °(13), (highly uncertain (41)), “(86), ©(41), (80), 9(81), "(82), '(83), 7(84). | 
| 


The logistics of liquid fluorine require study. The problem 
of transport has been solved by mobile tankage comprising : 
central zero loss container maintained below the boiling point 
of liquid fluorine by a jacket of liquid nitrogen and the whole 
surrounded by vacuum insulation (91). Such a container per- 
mits ground storage of liquid fluorine at economies comparable 
with that for storing liquid oxygen and has been demonstrated 
for cross country movement of liquid fluorine. The economics 
of shipping liquid fluorine from central manufacturing plants 
to several engine developments or launching areas has not, 
however, been worked out in comparison with the benefits of 
establishing fluorine generating and liquid nitrogen facilities 
close to the point of use. 

An engineering program to attain liquid fluorine engines of 
about 12,000-lb and about 80,000-lb thrust has been an- 
nounced (93). We can quickly estimate the order of magni- 
tude of fluorine demand for this development. The 80,000-Ib 
thrust engine would require some 250 lb per sec of propellant 
combination in a fluorine-hydrazine system; about two- 
thirds of this would be fluorine, or about 167 lb per sec during 
full thrust operation. Approximately one ton of fluorine 
then would be required for a 12 sec run, which may be 
taken as a typical test. With 20 runs per month on a given 
test stand and perhaps two parallel programs running simul- 
taneously, fluorine consumption at the rate of 480 tons per 
year would be required for this phase. The smaller engine 
would add another 70 tons. Overlapping of related programs 


Table 6 Comparison of fluorine and oxygen as oxidizers 
for hydrogen fuel* 


Mixture ratio, Mean 
Combination oxidizer: fuel propellant Specific Density 
density Impulse, impulse 
g/cc Isp sec Ispd 
F.—H:, 
max.Isp 4.5 0.32 381 123 
max. [spd 19.0 0.75 345 258 
max. Isp 3.5 0.26 372 97 
O.—Hb, 
max. I spd 8.0 0.43 326 139 


* Frozen exhaust composition, pressure ratio 600: 14.7; 
computed from data in (92). 


and later, full duration runs of perhaps 100 sec, could boost 
the fluorine requirement to perhaps 1000 tons per year. We 
note that such a program, while certainly requiring construc- 
tion of additional fluorine generating and liquefaction facili- 
ties, and perhaps also of some anhydrous hydrofluoric acid 
‘apacity, would in no way embarrass the United States’ 
fluorine chemicals industry, which is based on over 300,000 
tons per year of fluorine equivalent. 

Of the six known halogen fluorides, chlorine trifluoride and 
bromine pentafluoride are the best candidates for propellant 
use. The physical properties and the energy content of 
chlorine trifluoride make it interesting as a storable liquid 
oxidizer, though its violent incendiary nature is objectionable. 
Chlorine trifluoride already has field use, however, in a com- 
mercial metal cutting tool employed in deep wells (94). 

The place of bromine pentafluoride among storable pro- 
pellants is not yet certain. Its very high density indicates a 
superior density-impulse figure, but calculations depend on 
very uncertain estimates of the heat of formation of the com- 
pound, for which measurements are not yet reported (41). 

Superior characteristics are expected for perchloryl fluoride 
as a storable liquid oxidizer. The thermodynamic and 
physical data are accurately known, and specific impulse 
calculations show a performance approaching liquid oxygen. 
The compound is permanently storable in common materials 
of construction, is nonincendiary and nonexplosive (49). 
The toxicity is low (95), and there is no corrosive action on 
skin. The somewhat high vapor pressure and the thermal ex- 
pansion at higher temperatures may sometimes be objection- 
able (49). First discovered in 1952, this compound is now 
being made available in large quantities (96). 

Oxygen difluoride is a cryogenic liquid oxidizer with handling 
problems similar to those of liquid fluorine, but milder in re- 
spect to incendiary and corrosion properties. It is, however, 
highly toxic (97). The oxygen content is helpful when carbon- 
containing fuels are considered. Since the cost of oxygen di- 
fluoride will probably not be less than that of liquid fluorine, 
it is not yet clear in which situations the compound would be 
advantageous. 

Nitrogen trifluoride, in contradistinction to nitrogen tri- 
chloride, is a stable, inert, noncorrosive gas which requires 
low temperature liquefaction and storage. Good specific 
impulse can be calculated for this compound, but not much 
advantage over liquid oxygen is evident. 

Solid propellants: It is tempting to seek the benefits of 
fluorine in solid oxidizers. The energy of the carbon—fluorine 
bond is so high that the usual fluorine polymeric compounds 
would have limited use. The nitrogen—fluorine bond appears 
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to be less energetic, but no solid nitrogen—fluorine compounds 
have yet been reported. The highly endothermic oxygen—— 
fluorine bond would be an energetic component, but most 
compounds with this bond are explosive, and no solid com- 
bination has been reported (98). 

There are also the complex metal fluorides derived from 
the halogen fluorides, as for example KBrFy, BrF».-SbI’s. 
These are stable and tractable solids, but the presence of 
relatively heavy atoms detracts severely from calculated per- 
formance (99, 100). 

Other high energy uses: Although no further applications of 
the high energy storage value of chemical systems containing 
fluorine compounds have been reported, it is evident that 
opportunities in formulating explosives, monopropellants, in- 
cendiaries and the like, should be sought, and their effective- 
ness in energy conversion systems like electric primary and 
secondary batteries should be investigated. 


Chlorofluorohydrocarbons 


The commercial products known variously as Isotrons 
(101), Freons (102) or Genetrons (103) have recently ex- 
tended their usefulness in refrigeration and to other heat 
pumping cycles, with expanded ranges to higher (104) (e.g., 
Isotron 21) or lower (105) (e.g., Isotron 13) temperatures, and 
in pressure packaging of a very great variety of commodities 
(106). The packaged energy of these materials is also used in 
small mechanical drives (107). The use of Isotrons for foam- 
ing insulating fillers for temperature controlled spaces is de- 
veloping (108). An interesting new use for Isotrons is in sol- 
vent cleaning of systems to be used for liquid oxidizers, taking 
advantage of the wide selection of boiling points and the non- 
combustible nature of these materials (109). Some very 
promising work indicates that lubrication of moving parts at 
high temperatures can be effected by enclosing the parts in 
Isotron 22 (110). Effective Mach numbers in wind tunnel 
testing can be much increased by replacing air with the denser 
Isotron vapors (111). 

An important advance in control of aircraft fires has re- 
sulted from introduction of bromotrifluoromethane in fire ex- 
tinguishers (112). 


Fluorine-Containing Plastics, Elastomers and Oils 


The two older fluorine-containing plastics based on tetra- 
fluoroethylene (CF.==CF.) (113) and on _ chlorotrifluoro- 
ethylene (CFCI=CF,) (114) have gained in extent and variety 
of uses, especially in situations involving highly corrosive en- 
vironments or elevated temperature. Better fabrication 
methods and lower cost permit the manufacture of larger 
objects, such as diaphragms or plastic bags (115). Ad- 
vantage is taken increasingly of the low friction character in 
sliding parts, and of the water and grease repellency in surface 
coatings (116). By their nature, fluorocarbon plastics resist 
adhesion to other materials, but a new method using a surface 
modifying treatment now permits bonding of fluorocarbons to 
themselves or to other surfaces (117). 

These two basic polymers lack resiliency and freedom from 
creep under sustained strain, and successful research was 
carried out to develop highly fluorinated elastomers. Typical 
of these is one based on a copolymer of vinylidene fluoride 
(CHs=CF,) and perfluoropropylene (118). 
Another is a fluorinated butylacrylate polymer (119), and 
still another is a copolymer of vinylidene fluoride and chloro- 
trifluoroethylene (CFCI—=CF,) (120). A fluorinated silicone 
polymer, with superior solvent resistance, is also elastomeric 
(121, 122). A fluorinated polyester shows excellent low tem- 
perature elastomeric properties (123). 

New fluorocarbon oils are of much interest to the missile 
program, being thus far the only useful nonflammable lubricat- 
ing, sealing and hydraulic fluids other than water-based 
products. Present compositions, in use in liquid oxygen and 
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nitric acid handling equipment, are low polymers of chloro- 
trifluoroethylene (CCIF=CF,) (124, 125, 126). 

Recent reports describe oils derived by thermal telomeriza- 
tion of perfluoropropylenes (127, 128). Improvements in the 
temperature-viscosity relation are suggested by very recent 
work on vinylidene fluoride (129). 

Materials of medium range molecular weight, bearing func- 
tional groups, as sulfonic and carboxylic acids, alcohols, ete., 
are being promoted for imparting water and grease repellent 
characters to susceptible surfaces (130, 131). 


The Future 


Fluorine chemistry, once frustrating to the experimental 
chemist, is now well in hand. The principles governing reac- 
tions and properties of fluorine compounds are understood. 
A variety of techniques is available for conducting reactions, 
for analysis and structure determination, and for large-scale 
preparation and transportation of fluorine products. Fluorine 
chemistry is mature. 

But discovery continues. Especial advances can be ex- 
pected in tailor-making organic fluorine polymers, in synthe- 
sizing new and still relatively simple inorganic structures, in 
better understanding the kinetics and mechanism of reactions 
involving the fluorine atom, and in best employing fluorine 
chemical products in the new exciting surge of extra- 
terrestrial exploration. 
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At times of solar flares, bursts of extremely penetrating x-rays reach down to the 
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When photodetectors in rockets were pointed 


output of a great flare may reach a level as high as one 
thousandth of the solar constant. Flares generate a sequence 
of events in the Earth’s ionosphere from its lowest fringe to its 
outermost reaches that profoundly affect ratio communica- 
tions, distort the geomagnetic field and create the brilliant 
Those disturbances of the ionosphere which appear 


‘Superintendent, Atmosphere and Astrophysics Division. auroras. 

Fellow Member ARS. instantaneously are associated with flashes of short wave 
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electromagnetic radiations, ultraviolet and x-rays. Delayed 
effects are produced by the later arrival near the Earth of 
clouds of solar gas blown off by the flare. 

To understand the nature of the forces and processes in- 
volved in the flare mechanism, one must study all the charac- 
teristic flare emissions from cosmic rays to radio noise bursts. 
The Rocket Flare Patrol (1),2 carried out on San Nicolas 
Island from July to September 1957, was one of many IGY 
solar flare study programs. It was designed specifically to 
discover whether x-ray and ultraviolet radiations are pro- 
duced in flares and are responsible for sudden short wave 
radio fadeouts. Since these rays are absorbed at the altitude 
of the D-region ionosphere, rockets were needed to carry 
the detectors to heights greater than 75 km. Furthermore, 
since flares are short-lived transients that give no advance 
warning of their development, constant readiness and push- 
button launching within a matter of a minute were prime 
experimental requisites. 

The major accomplishment of the Rocket Flare Patrol 
program was the detection of intense x-ray flashes at suf- 
ficiently short wave lengths to produce the low lying D-region 
ionization required to explain radio fadeout. From the astro- 
physical standpoint, the x-ray production can be understood 
as a consequence of heating the coronal atmosphere above the 
visible flare to temperatures greater than 10 million deg K. 

Although the x-radiation of a flare can traverse the sun to 
Earth distance with little absorption or scattering, clouds of 
particulate matter ejected by the flare interact strongly with 
the interplanetary and geomagnetic fields. Basic to any 
description of the motion of charged particles from sun to 
Earth is a knowledge of the density and spatial distribution of 
the interplanetary ionized gas. One of the most important 
developments in the understanding of solar-terrestrial rela- 
tionships during the past few years has been the accumulation 
of evidence that interplanetary space contains a much higher 
density of ionized gas than was previously supposed. Rather 
than the old astronomical estimate of one particle per cm‘, 
it appears that densities of about 100 to 1000 per cm are more 
appropriate. Theories of whistling atmospherics, the extent 
of the solar corona, and estimates of the emission of matter 
from the sun require the higher densities. 

Recently, a promising new field of “Rocket Astronomy” 
has been opened up (2) in which photometric surveys can be 
made of the invisible ultraviolet emissions of interplanetary 
gas, galactic nebulosities and hot stars. The first rocket 
observations of the light of the night sky in the far ultraviolet 
appear to provide new evidence for higher particle densities 
in interplanetary space. Above the atmosphere the entire 
sky is diffused with a glow of ultraviolet light scattered back 
to the night hemisphere of the Earth from interplanetary gas. 
It appears that the original source of this light is solar radia- 
tion in the resonance wave length of excited hydrogen atoms, 
d 1216 A, known as Lyman a. The sun radiates this wave 
length into space where cold, neutral hydrogen atoms absorb 
and instantly reradiate some of it back to the dark side of the 
Earth. From the intensity and angular distribution of the 
nighttime Lyman a flux, rough estimates of the neutral 
hydrogen content can be made. The indications are that the 
neutral atom density may be as great as 1 per cm* and that 
the ionized content should be a few hundred times as great 
near the Earth. 

From a height of 146 km the rocket’s photoelectric sensors 
measured a bright albedo of the atmosphere in this same wave 
length, indicating the presence of a layer of free hydrogen 
atoms above 80 km. The combined fluxes of extraterrestrial 
ultraviolet light and the albedo exceeded the input of total 
visible starlight. The energy contributed to the nighttime 
ionosphere by Lyman a@ was found to be comparable to the 
entire energy dissipated in the visible night airglow. As a 
source of nighttime ionization at ionospheric altitudes, the 


2 Numbers in parentheses indicate References at end of paper. 
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observed Lyman a@ flux is capable of producing as much 


ionization as is derived from meteor trials. 


The Solar Corona 


A few decades ago, astrophysicists still viewed the sun as a 
glowing sphere of hot gas, radiating into space like a black 
body at a temperature of 6000 K. This picture was inade- 
quate to explain the very short wave length ultraviolet or 
x-ray emissions required for the formation of the ionosphere. 
Only with the discovery that temperatures in the coronal 
atmosphere reach a million deg K or more, did the mystery 
of solar control of the ionosphere begin to unravel. These 
high temperatures also explain the steady radio noise of the 
sun on meter wave lengths. 

At the time of a solar eclipse, the corona is visible as a thin 
white halo with a slightly greenish cast spreading millions of 
miles out into space. The simple fact that the corona 
extends to such great distances requires that its temperature be 
of the order of a million deg K to balance the tremendous 
gravitational pull of the sun. The white light halo may |e 
compared to the halo around a street lamp on a foggy night. 
In the latter case, the lamplight is scattered by a mist of 
water droplets, whereas in the corona the light of the sun’s 
photosphere is scattered by fast electrons. The fact that tlie 
dark Fraunhofer lines do not appear in the corona means that 
the electrons are moving so fast that the Fraunhofer lines are 
completely washed out by Doppler effect. The required 
temperature of the electron gas is close to a million deg K. 

In addition to the white light corona, spectroscopic studies 
reveal an emission line corona whose integrated brightness is 
about 2 per cent of the coronal light. Some 31 lines are ob- 
served between wave lengths of 3328 and 10,798 A. For 
almost half a century after their discovery, the origin of these 
lines was not understood, and they were attributed to an 
element “coronium” unknown on Earth but supposed to 
exist on the sun. We now know that these lines originate in 
highly ionized metal atoms of the coronal atmosphere. The 
brightest line is the green wave length \ 5303 A of Fe XIV, 
which is an iron atom stripped of half of its normal quota of 
26 electrons. Theory requires a temperature of a million 
degrees to produce maximum brightness of this line, which is 
the source of the greenish tinge of the coronal light. Other 
prominent lines are the red line of Fe X and the yellow line of 
Ca XV. The coronal red line of Fe X appears uniformly 
distributed over the sun and changes little with the 11 year 
solar cycle. It is representative of the base level of coronal 
temperature, about 500,000 K, characterizing a quiet sun. 
The green line of Fe XIV is distributed in patches which 
brighten over active sunspot regions and reflects conditions of 
higher local temperatures up to a few million deg K. In eases 
of extreme excitation, such as solar flares, the yellow line of 
Ca XV appears, often indicating temperatures in excess of 
5 million deg K. 

The ionization potentials of the metal ions are 233 ev for 
Fe X, 355 ev for Fe XIV and 814 ev for Ca XV. Quanta of 
these energies have wave lengths of 48, 35 and 15 A respec- 
tively. It is immediately obvious that recombination of ions 
and electrons, plus electron impact excitation of line emissions 
will produce a spectrum of wave lengths extending 
well into the soft x-ray region. Furthermore, the intensities 
of the visible coronal emission lines should be an index 
of the quality and intensity of the invisible x-ray emis- 
sion. This has been verified by direct n 
rockets (3). 7 


X-Ray Emission of a Quiet Sun 


The ‘Rocket Spectrum” of a quiet sun extends to x-ray 
wave lengths as short as 20 A which can penetrate to about 
100 km with the sun overhead. These x-rays can originate 
only in the million deg K gas of the solar corona. To a first 
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approximation the shape of the spectrum may be described 
as a Planckian distribution with a maximum near 50 A, even 
though the coronal gas is so tenuous that its emissivity is only 
10-"*, if we treat it as gray body radiation. At times of 
coronal activity, evidenced by glowing patches of green and 
yellow line emission, the total x-ray emission increases, and the 
short wave length limit may reach 6 or 7 A. This enhanced 
x-ray emission may be attributed to local hot spots or ‘coronal 
condensations.” By observing the increase in decimeter 
wave radiations from a condensation, it is possible to map its 
extent and derive an estimate of the local electron density. 
Radio interferometers show that the condensation areas cover 
several thousandths of the solar disk. The electron density 
appears to increase from 10 to 20 times normal, and the x-ray 
emission should accordingly increase several hundredfold. 
When integrated over the solar disk, the x-ray contributions 
from these local hot regions raise the total x-ray output by a 
fxctor of two or three. 

The behavior of a quiet sun over a solar cycle is adequately 
explained by the combination of a uniformly distributed 
general x-ray emission characterized by a half million degree 
temperature on which is superposed intense emission from 
coronal condensations, corresponding to local temperature in- 
creases of 5 X 10° deg K or higher. These x-ray emissions 
ionize all constituents of the atmosphere equally effectively 
and may account for the major portion of the E-region iono- 
sphere and much of F-region. X-ray measurements from 
rockets (4), although still comparatively few and scattered 
over the past 10 years, confirm a general correlation with 
coronal activity and with slow variations in the E- and F-re- 
gions of the ionosphere. For a period of two to three years 
about the maximum of the 11 year sunspot cycle, radio com- 
munications on short wave bands improve as the electron 
density of the E- and F-regions increases by a factor of 1.5 to 
2 compared to sunspot minimum. The source of this in- 
creased ionization is the increase in short wave length ultra- 
violet and x-ray emission from the corona. 


Lyman a Radiation and D-Region 


Although it may seem to be contradictory at first thought, 
it is possible for increased ionization of the upper atmosphere 
to decrease the strength of radio reflections. If the ioni- 
zation is produced above 100 km, in the E- and F-regions, the 
air is so thin that the electrons behave to a first approxi- 
mation as free electrons. They vibrate in response to the 
incident radio waves and reradiate the same frequencies 
without attenuation. At altitudes below 100 km, the air 
density becomes high enough to increase the collision fre- 
quency to the point where an electron cannot radiate the 
energy it picks up from the incident radio waves before it 
collides with an air molecule. The collision damps out the 
electron vibration and absorbs the radio energy. The D- 
region of the atmosphere, 75 to 90 km, normally contains 
only a few thousand electrons per cc near 90 km, but these are 
sufficient to produce detectable absorption of megacycles 
radio waves and cosmic radio noise. At the time of a flare, 
the D-region ionization increases to the extent that fre- 
quencies of 5 to 20 megacycles per sec may be completely 
absorbed, a condition known as radio fadeout. 

When the sun is quiet, rocket measurements rarely show 
any x-ray penetration below 100 km. In the far ultraviolet, 
only Lyman @ (A 1216 A) is detected in D-region (5). Al- 
though it cannot ionize the major constituents of the atmos- 
phere it is a very efficient ionizer of nitric oxide. The flux 
of Lyman a is so great that less than a part per million of NO 
suffices to provide the D-region ionization. 

Only on a few rare occasions have flares been seen in white 
light. Ordinarily, a flare is undetectable unless viewed with 
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a spectroheliograph or through a birefringent filter that pasess 
only the hydrogen red line (Balmer a) at 6563 A. It then 
appears as a local brightening of the solar surface in the 
neighborhood of a sunspot group. Within a matter of 
minutes, the flare may expand over 100 million to a billion 
square miles and reach a brightness 10 times that of the nor- 
mal photosphere. After flashing to its peak the flare decays 
more slowly in half an hour to several hours, depending on its 
size. Flares are classified according to area covered as 
Class 1, 2 or 3 in order of increasing size, with the addition of 
+ or — to indicate larger or smaller than average in any 
class. 

Simultaneously with the appearance of the visible flare, 
there occurs a fadeout of short wave radio reception over the 
sunlit hemisphere. Radio fadeout and other Sudden Iono- 
spheric Disturbances (SID) generated by flares can all be 
traced to increased ionization of the D-region. The flare 
must therefore generate Lyman a or x-rays of wave lengths 
shorter than 10 A, either of which can ionize the D-region. 
Radio absorption measurements during an SID require a five 
to tenfold increase in total ionization of D-region for Class 2 
and 3 flares. It is possible that Lyman a alone could produce 
the required increase in absorption, but other D-region phe- 
nomena are more difficult to account for by Lyman a. 

Very low frequency radio waves (16 keps are reflected from 
the base of the ionosphere where the electron density is about 
200 per cm*. The height of this level of electron density can 
be measured by observing the propagation of 16 keps waves 
between a transmitter and receiver about 100 miles apart via 
the ground and via reflection from the ionosphere. The phase 
difference between the two waves arriving at the receiver will 
vary with vertical movement of the reflecting region. During 
a flare, a Sudden Phase Anomaly (SPA) can be observed as a 
result of a drop in the height of the 200 electron per em* 
level. During a large flare, the SPA may correspond to a 
drop of 16 km. To account for this drop as the result of 
ionization by Lyman a requires a 10 thousandfold increase in 
solar flux (6). Since this radiation must come from the flare 
region, which is about 10~* of the area of the disk, the flux 
from the flare region must increase by 107. Such an increase 
is impossible to account for by any astrophysical mechanism. 
On the other hand, an x-ray flux of only 107° erg em~*sec™! 
at a wave length of 1 or 2 A could produce the 16 km 
drop indicated by the most extreme SPA. 


The Rocket Flare Patrol 


The Rocket Flare Patrol was undertaken primarily to settle 
the question of the relative importance of x-rays and Lyman 
a in the radiation of a flare. Since flares occur without ad- 
vance warning, it was essential that instrumented rockets be 
constantly in readiness for instant activation and push- 
button launching. A two-stage combination of solid pro- 
pellant rockets, a Nike booster plus a Deacon fulfilled the 
requirement of lifting a 30 lb payload to about 350,000 ft. 
The rocket base was established on San Nicolas Island, 65 
miles off the California coast on the Navy’s Point Mugu 
Range. Decisions to fire were based on observations of 
radio fadeout and on visual observations reported via radio 
and teletype links to the Island from observatories on Mt. 
Wilson, Calif., Sacramento Peak, N. Mex. and Climax, 
Colo. In charge of the operation on San Nicolas Island were 
NRL scientists, Dr. T. A. Chubb, Dr. J. E. Kupperian Jr. 
and Dr. J.C. Lindsay. Engineers of the Cooper Development 
Corp. handled the preparations of the two-stage rockets for 
launchings. 

Each Deacon rocket was equipped with a photon counter 
to measure x-rays (1 to 8A), an ion chamber to measure 
Lyman a@ (1216A) and a pair of visible light photocells to 
provide rocket aspect throughout flight. Information was 
telemetered continuously throughout flight on four channels 
of FM-FM telemetering. 
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Eleven launchings were carried out at times when solar 
flares were observed. In spite of varied misfortunes, suf- 
ficient data were obtained on all three classes of flares to 
answer the basic scientific questions. The flares produced 
marked increases in x-ray radiation of shorter wave lengths 
than are produced by a quiet sun. The largest flare produced 
the shortest wave length x-ray emission, about 2 A. The 
total x-ray flux varied from about 10~* erg cm~*sec™ for a 
Class 1+ flare to about 10~? erg cm~*sec™ for a Class 2+ 
flare. Although the total flux for a Class 3 flare was not 
determined because that particular rocket barely reached the 
D-region, it appears that the hardening of the spectrum rela- 
tive to the smaller flares may play a more important role than 
the increase in total flux in accounting for the associated 
SID. No evidence was obtained for large increases in Lyman 
a. However it is possible that Lyman a may have increased 
earlier in the flare process than the timesof firing of the rockets. 
In each shoot an attempt was made to time the rocket firing 
to coincide with the maximum phase of the flare. In each 
firing, therefore, the first few minutes of flare development 
were bypassed. 

These rocket measurements of x-ray flashes in solar flares 
can be explained as a consequence of intense heating of the 
coronal atmosphere above the visible flare. As was explained 
above, the temperature of a quiescent corona is about half a 
million deg K. At that temperature, the M shells of all 
atoms are stripped of electrons, and the corresponding ion- 
ization potentials fall in the range of a few hundred electron 
volts. The x-ray emission associated with such an electron- 
ion plasma lies in the wave length range from about 20 to 
190 A. (The innermost shell of an atom, containing two 
electrons when filled is designated the ‘‘K”’ shell; the “ZL” 
shell next surrounding the “‘K”’ shell contains 8 electrons when 
filled; and the ‘/” shell, with at most 18 electrons, sur- 
rounds the “7” shell.) At 5 million deg K, even the K shells 
of the lighter elements up to oxygen may be stripped cor- 
responding to ionization potentials of the order of 1500 ev 
or a wave length of 8 A. The 2 A x-rays measured in a large 
flare indicated that much higher temperatures are generated 
by the flare. At temperatures above 10 million deg K, even 
the comparatively heavy iron atoms lose all their electrons. 
The resulting recombination radiation has a series limit at 
1.4 4. Since the abundance of elements heavier than iron is 
negligibly small, thermal excitation of coronal x-ray emission 
should exhibit a short wave length limit close to 1.4 A. 

In support of the above picture of thermal! excitation of the 
x-ray flare, are the observations of Sudden Phase Anomalies. 
The largest flares have produced drops as great as 16 km in 
the level of the base of D-layer. Sixteen km is about the 
difference. in penetration of the atmosphere by Lyman a 
and by 1 A x-rays. 

For the future, we may anticipate that much more informa- 
tion about x-ray and ultraviolet flares will be obtained from 
satellites. Although the IGY Vanguard program is limited 
in its scope, increased payload in subsequent programs should 
greatly enhance the yield of data. At the same time, there 
is a great deal more to be learned from rocket experiments. 
Satellites provide a continuous observing platform that is 
ideal for the detection of infrequent transient events, but the 
vertical rocket flight offers a simple means of obtaining spec- 
tral distribution by observing the variation in rate of absorp- 


tion vs. altitude. ad 

Since the beginning of upper air research with rockets, as- 
trophysical measurements have been confined almost ex- 
clusively to solar spectroscopy. In view of the success 
achieved in measuring the sun’s ultraviolet and x-ray emis- 
sion with photon counters, we were encouraged to search for 
these radiations from celestial objects in the night sky. The 
sun is a comparatively weak star; it is classified as a yellow 
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Fig. 1 Lyman a directional intensity contours when the 
detector looked at the upper hemisphere. These data were ob- 
tained from the portion of the flight above 130 km and are un- 
corrected for the 1225 to 1350 A discrete sources. The small- 
est intensity contour circle contains the antisolar direction 


dwarf and is long past its most brilliant phase. Early stars 
such as the blue white are so much hotter than the sun that 
maximum brightness ought to occur in the far ultraviolet. A 
simple rocket experiment was attempted in November 1955 
to determine the possibility of obtaining a map of stella: 
magnitudes in wave lengths of Lyman @ and the adjacent 
ultraviolet near 1300 A. The observed fluxes were unex- 
pectedly high and saturated the photon counters after they 
had penetrated only a few kilometers into the D-region. 
From that experiment it was possible to conclude that the 
major component of nighttime ultraviolet fell within a wave 
length band of 1050 to 1225 A ine luding the Lyman a wave 
length (1216 A) and a lower limit of 3.4 X 10-4 erg em~*see~! 
was assigned as the flux (2). At wave lengths between 1225 
and 1350 A discrete celestial sources were detected. 

With these early results as a guide, a much more elaborate 
experiment (7) was carried out on March 28, 1957 with NRL 
Aerobee 31. The rocket was launched at 2156 MST and reached 
an altitude of 146km. The ultraviolet telescopes consisted of 
bundles of hypodermic tubings through which the sky could 
be viewed over an angular field of 3 deg by sensitive photon 
counters. With less sensitive detectors, such as ionization 
chambers, wider collimators were used providing fields as 
broad as 0.5 steradian. The detectors were mounted with 
their collimators looking outward in directions normal to the 
rocket’s long axis. A spatial scan of the sky was produced by 
the combined spin and yaw motions of the rocket. The 
rocket behaved in almost ideal fashion, approximating an 
arrow trajectory with horizontal altitude near the peak of the 
flight and rolling with a period of 16 sec. Almost a complete 
scan of the sky was therefore obtained. 

Numerous discrete sources were observed in the 1225 to 
1350 A region. The determination of the direction of view of 
each detector at each instant of flight i is an exceedingly tedious 
job. It involves fitting together in a self-consistent manner 
the star signals from visible light photometers and positions 
of the visible airglow horizons. The end-product of this 
analysis is a solution of the aspect problem with an accuracy 
of 1 deg, permitting the construction of a map of the 
ultraviolet sky near 1300 A with a resolution of 3 deg. The 
brightest single source observed was the Orion complex which 
delivered about 2 X 10~* erg cm~*sec~! at the top of the at- 
mosphere. A complete description of this work is being pre- 
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pared for publication the Astrophysical Journal (9) and will 
not be discussed here, except to compare the radiation flux 
observed with other energy sources in the night sky. 


Interplanetary and Atmospheric Hydrogen 


The detectors whose bandwidth included the Lyman a line 
showed a very different picture of the sky. At an altitude of 
75 km an intense glow, diffused over the entire sky, was seen 
above the rocket. At 85 km, an additional source became 
evident, originating in the atmosphere below the rocket. 
Both fluxes increased up to 120 km and then remained con- 
stant to the peak of the flight, 146 km. Both the atmospheric 
and extraterrestrial sources of radiation have been identified 
with the Lyman @ emission of hydrogen. Furthermore, this 
radiation appears to be of solar origin; the flux from outside 
tle atmosphere is scattered by atomic hydrogen in space; 
the flux from below the rocket is scattered a second time by 
hydrogen in our atmosphere. The observed albedo was about 
42 per cent. 

Fig. 1 is a polar plot of the intensity isophotes obtained 
with a detector collimated to view a field of 16 deg half width. 
The center of the plot is the zenith. The minimum contour 
ovcurs in the antisolar direction, i.e., looking directly away 
from the sun. This a plot of the raw data and the effect of the 
unsubtracted 1300 A radiation from Orion, and the galactic 
plane can be seen in the West. From the presence of the 
minimum in the antisolar direction we can conclude that the 
source of this radiation is the sun. 

That the influx of Lyman a is scattered to the night side of 
the Earth by neutral hydrogen in space is deduced from the 
high albedo of hydrogen in the Earth’s atmosphere. Lyman 
a from the sun is a broadened emission line characteristic 
of the high temperature of the solar atmosphere. The scatter- 
ing cross section of the cold hydrogen for the solar line is small 
except for the very center of the line. As a consequence no 
detectable albedo has been found in daytime measurements 
with the solar flux of Lyman a directly incident on the at- 
mosphere. On the other hand, if Lyman a in the night sky is 
once scattered by neutral hydrogen atoms, it has been mono- 
chromatized by the scattering process to include only the core 
of the original solar line. When this highly monochromatized 
radiation enters the denser hydrogen region of the Earth’s 
atmosphere near 80 km, it is backscattered with the observed 
albedo of 42 per cent. 

Estimating the neutral hydrogen density of interplanetary 
space from these observations is complicated by our lack of 
detailed knowledge of factors, such as the solar line profile 
and particularly its central intensity. From the mere fact 
that the distribution still shows a minimum when viewed 
from the earth (Fig. 1) it is possible to conclude (10) that the 
upper limit of neutral hydrogen density is about one atom per 
cm’, In the vicinity of the Earth, the ionizing solar flux and 
the rate of recombination are both known. These figures lead 
to the conclusion that the neutral atomic hydrogen concentra- 
tion is between 10~* and 10~* of the proton density. We con- 


Table 1 Radiation sources in the day and night skies 


Source Ergs per cm~*sec~! 
sun a! 1.4 X 108 

full moon 3000 10-8 
total starlight 1:3 x< 16-* 
airglow (visible) 16 X 10-3 
OH (infrared) 10-* 
Lyman a 10 10-* 
celestial sources (1230- 0.1 X 10-3 
cosmic rays 3.8 


clude therefore that the ionized particle density may be a few 
hundred protons per ce. 

Future experiments are being planned to determine the 
contour of the solar Lyman a line with high resolution, to 
determine the daytime albedo, and to complete the measure- 
ment of the angular distribution of scattered Lyman a@ over 
the entire sky by making observations during an eclipse. 

The total flux of Lyman a at night on a horizontal plane 
surface is 10~? erg em~'sec~! when integrated over the entire 
hemisphere. This flux plus the contribution of the albedo is 
capable of producing at night almost one tenth the electron 
density of D-region produced by the daytime Lyman a flux, 
if the same ionization and recombination processes are opera- 
tive both day and night. As a source of ionization at alti- 
tudes near 100 km, Lyman a is capable of ionizing the Ca 
atoms accumulated at that level from the influx of meteorites. 
Accepting Nicolet’s value (8) for the Ca atom density at 100 
km, it can be shown that the rate of ionization by nighttime 
Lyman a@ is comparable to the ionization contributions from 
the ablation of meteorites as observed by radar echoes. 

Table 1 is a list of the solar constant and of the various 
sources of energy in the night sky. It is interesting to note 
that the newly discovered Lyman a contribution exceeds the 
sum total of visible starlight and is comparable to all the 


visible airglow. 
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is presented. 


HE STRUCTURE of the upper atmosphere is influenced 

by photochemical dissociation of its molecular con- 
stituents. This, as well as chemical reactions and diffusion 
among the species, gives rise to the complex composition of 
the layer above 90 km, known as the chemosphere (1).? 
Fig. 1 shows some of the pertinent data (2, 3, 4) on com- 
position, pressure and temperature of the lower portion of 
this layer. To the extent of our knowledge, the atmosphere 
still contains nitrogen and oxygen in a 4:1 ratio (on the basis 
of molecular concentration); the concentration of ozone, NO 
and atomic nitrogen is negligible, but the concentration of 
atomic oxygen, which is dissociated by absorption of the 
ultraviolet solar radiation, increases rapidly from about 90 
km. At 150 km almost all the oxygen is in free radical form. 
The heavy line in Fig. 1 shows the values used in the calcula- 
tions that follow; it is thought to represent a reasonable high 
limit of the degree of dissociation of oxygen with altitude. 

The temperature of the chemosphere is relatively low. 
Consequently, the atomic oxygen is not in thermochemical 
equilibrium with its surroundings. If the atmosphere were 
shielded from the solar radiation, the radicals would recom- 
bine. However, the density of the gas is so low that the rate 
of chemical reaction is extremely slow. The layer of dis- 
sociated oxygen is quite permanent in spite of variations in 
the radiation field. 

The energy stored in this photochemical system is sub- 
stantial, and the possibility of using it for propulsion of mis- 
siles and satellites has been proposed repeatedly. An 
“oxygen recombination ramjet” (5, 6, 22) would induce and 
accelerate the recombination process within a controlled 
thermodynamic system and convert the internal energy re- 
leased by the reaction into kinetic energy of the exhaust jet. 

The evaluation of the potentials of such a propulsor is made 
difficult by the relative uncertainty about the composition 
of the atmosphere, the kinetic reaction rate constants and the 
aerodynamics of flight through rarefied gases. Accordingly, 
the assumptions made in the following discussion are clearly 
stated, and the basic arguments are presented so that the 
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An analysis of the oxygen recombination ramjet for propulsion in the Earth’s upper atmosphere 
Its performance depends on the chemical kinetics of the oxygen recombination reac- 
tion, the aerothermodynamics of the ramjet cycle and the rarefied gas aerodynamics. 
a limited flight range, these phenomena imply the existence of mutually exclusive conditions. Near- 
orbital speeds appear to be unattainable, and only low supersonic flight is possible, and only if 
appropriate catalysts to accelerate the reaction can be found. The performance of the engine 
depends primarily on the flight Mach number, and it is difficult to justify speculations that « 
better knowledge of atmospheric structure or of reaction kinetics may reverse these conclusions. 


A. F. CHARWAT! 


University of California 
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Outside of 


reader can form his own conclusions. A basic thermodynami:: 
analysis, about which there is not much doubt, shows that 
the range in which an oxygen recombination ramjet can 
operate is rather narrow and restricted to low supersonic 
flight Mach numbers; flight at near-orbital speeds (6) is 


200 


L 
130 150 
ALTITUDE KM. 
Fig. 1 Temperature, pressure (from (4)) and degree of oxygen 
dissociation of the upper atmosphere from 80 to 160 km 
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where W is the molecular weight of the mixture (gr per 


proved not to be possible. In the range in which this sys- 
tem will produce a net thrust, its performance is very mar- 
ginal. Although one cannot positively disprove its useful- 
ness or its feasibility, this author doubts that it can be built 
successfully. 


The basic oxygen radical recombination requires three-body 
collision. Not enough is known about the relative efficiency 
of various third bodies to warrant differentiating among 
them. Consequently we assume 


Analytical Formulation* 


where the third-body concentration M is 


(M) = (Nz + O: + O) 


[2] 


A wall or any solid surface can also serve as the “third body”’; 
this will be further discussed in a subsequent paragraph. 

In addition to the direct recombination, several chain 
mechanisms which lead to the recombination of oxygen 
atoms can be written. Their importance is emphasized 
when they increase the speed of the process of recombina- 
tion. In particular, the intermediate formation of NO, is 
often quoted as such a mechanism. This chain requires the 
existence of NO in the gas. However, since there is no evi- 
dence of it in the atmosphere at altitudes under consideration 
here, it must be considered in conjunction with the synthesis 
of NO in the engine proper. It is then found that these 
secondary mechanisms of recombination cannot be important 
in comparison with the basic three-body recombination, at 
least at gas temperatures around 1000 K, which are of in- 
terest here. Accordingly, only the basic three-body process 
given by Equation [1] is considered in subsequent calculations. 

The degree of dissociation of the oxy _— be defined as 


The mole fractions of (O), (Oz) and “third bodies” (M) are 
2a/(5 + a), (1 — a)/(5 + a) and 1 mole of constituent per 
mole of mixture, respectively. The ideal gas law is assumed 
to hold 


4 
G+a) 


mole). Dissociation of nitrogen is negligible at the thermo- 


dynamic conditions under study, and the concentration of | 


the other possible compounds is assumed to be small, as ex- 
plained above. 
The rate law is (see, for instance, Eq. 3.2 of (17)) 


It can be expressed in terms of a as 


da Kpa 4a? 

— = k,(—) (1 — a)| ———___| [6 

dt (=) (1 = + 
where Kp, is now the dissociation equilibrium constant based 
on partial pressure, and it is related to the other quantities as 


ka 4a,? P 
RT)- = 7 


= (RT) — 


* A list of definitions of symbols appears at end of paper. 
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The equation can be integrated if a functional relation be- 
tween the state variables P, T and a is specified (in addition 
to the equilibrium and the rate constant). We shall restrict 
our attention, for the present, to recombination at constant 
pressure. Then dP = 0, and the relation between T and a 
is given by the energy equation 


dh = qda = C,dT [8] 


C, is an average specific heat of the reacting mixture. For 
the particular case in question, it varies so little between the 
extreme limits of composition that we may consider it con- 
Numerical 


stant. gq is the specific energy of dissociation. 
values used in subsequent calculations are 
Cy = 0.25 cal/gr of air deg K _ [9] 


q = 815 cal/gr of air 


It remains to specify the values of the two constants. The 
equilibrium constant is fairly well established. Since the 
subsequent analysis is concerned with the temperature range 
between 1000 and 2000 K, some of the most rec ent ds lata (8) 

(1956) were used to evaluate it in that range > 


= 26,512 6-6 
= 7.944 x 10° xX 10 T atm [10] 


The rate constant k, is known with considerably less pre- 
cision. Experimental data are lacking, and estimates vary 
by two orders of magnitude. In (9) shock tube experiments 
(high temperatures) in a 1 per cent oxygen, 99 per cent argon 
mixture indicate for k, the magnitude 10! (cc/mole)? sec. 
Feldman infers (from Prandtl-Mayer expansion observa- 
tions) an upper limit between 10!8 and 10 ((10 and 11) re- 
spectively). Other estimates are made (12, 13, 14, 15) based 
on theory and on measurements of the rates of similar reac- 
tions, notably bromine and iodine recombinations. 

The temperature dependence of the recombination constant 
is even less well-known. Available analyses (a brief discus- 
sion and assessment of the arguments is reported in (15)) 
predict dependencies from the inverse square root to the 
5/2 power. Experimental data for oxygen in air do not 
cover a wide enough range to permit a good evaluation of this 
function; the authors of (9) feel that the inverse square root 
correlates their data best. Somewhat better known are the 
similar iodine and bromine recombinations, for which the in- 
verse 3/2 power dependence is usually assumed. 

For the present discussion we assumed the following rela- 
tion for the recombination rate constant 


1 \ 
This formula yields the higher rate in the range of tempera- 
ture covered by the most recent experiments (9) and exhibits 
the strongest reasonable increase in rate with decreasing tem- 
perature. With this upper value, the slowness of the process 
cannot be ascribed to the uncertainty in the choice of the rate 
constant. 
The aerothermodynamics of the propulsive duct must also 
be considered. The stagnation pressure and enthalpy rela- 
tive to the moving ramjet are uniquely related by 


¥, 
-( +4 2 at) 


where ne is the inlet compression ratio, a function of the geom- 
etry of the inlet (and the resulting shock structure in super- 
sonic flow). If partial compression is considered, then the 
internal duct Mach number M, and the local pressure and 
temperatures in the duct are similarly related 
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Equations [12 and 13] assume a ‘frozen’’ compression. 

The expansion in the exhaust nozzle is also assumed to be 
chemically frozen. Subsequent discussion will show that 
these assumptions are reasonable in view of the slowness of 
the recombination process in comparison with the time of resi- 
dence of the gas in a diffuser or a nozzle of reasonable length. 

The nonequilibrium expansion presents, incidentally, none 
but computational difficulty. A solution can be obtained by 
eliminating a between the rate Equation [6] and the energy 
equation for (externally) adiabatic flow.* The jet velocity 
can be computed from the relation 


The recombination will be considered to proceed isobari- 
cally, and it will be assumed that the areas of the chamber are 
adjusted to maintain this condition. The complicated rela- 
tions governing the aerothermodynamics of flow in a duct with 
heat addition need not be included. This will be justified in 
further discussion. 

For frozen compression and expansion and “isobaric com- 
bustion”’ the well-known relation for the specific impulse (net 
thrust per unit weight of airflow) of a ramjet is, in terms of the 
degree of recombination that took place in the chamber 


[14] 


(7 3260(a — ar) 
VT = 4.56 T. To 


7 


Here nz is again theram compression recovery ; ng is the burner 
pressure ratio, and 7, is the exhaust nozzle thrust coefficient. 


Self-Limiting Character of the Recombination 


Recombination of the oxygen radicals can proceed only 
until the composition of the gas reaches equilibrium at the 
existing temperature; Fig. 2 shows the characteristic plot of 
equilibrium composition at various pressures and tempera- 
tures in the range of interest. As an extreme example, as- 
sume the ambient atmosphere is totally dissociated (ao = 1), 
and that recombination begins at the lowest tempera- 
ture, that is, at ambient temperature. This implies that there 
is no ram compression (the engine consists simply of a super- 
sonic duct which shields the air from radiation for a sufficient 
time to allow the gas to reach thermochemical equilibrium). 
As recombination proceeds, the gas temperature is raised by 
the release of chemical energy. If the pressure were 1077 atm 
(corresponding to 90 km, the lowest altitude at which disso- 
ciated oxygen exists at all), the maximum fraction that can 
recombine before the reaction is stopped is A = (a) — az) = 
0.62, and not unity. To cause further recombination the 
temperature must be decreased. This cannot be done by ex- 
pansion because the pressure in the duct is ambient. Compres- 
sion, on the contrary, leads to redissociation because of the 
associated temperature rise. Dissipation of heat to the en- 
vironment causes recombination only as a ‘‘make-up”’ of the 
energy lost, and the propulsive gas can in no way acquire 
more thermal energy transformable to thrust work than the 
limit above. 

Thus this ‘“combustion”’ is limited not only by composition 
(the “fuel” available in the atmosphere) but also by tem- 
perature of the combustion products. This is a new con- 
sideration in the experience of the powerplant engineer. 

Similarly, consider the opposite extreme whereby the gas is 
compressed to stagnation conditions (frozen flow compression) 


3A detailed treatment of this, for a similar chemical recom- 
bination, can be found in the article by S. S. Penner and D. 
Altman, J. Franklin Inst., vol. 245, 1948, p. 421. 
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Fig. 2 Equilibrium degree of dissociation of oxygen at low 
pressures as function of temperature 
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before allowed to recombine. It is easily seen from Fig. 2 
that the temperature rise associated with adiabatic compres- 
sion will soon (for instance, M = 6.5 at 100 km and isentropic 
compression) permit no recombination at all. 

The complete solution of this problem is a function of am- 
bient pressure, temperature, composition and the degree of 
ram compression. One may suspect an optimum point at 
some partial compression of ram air: The pressure rise favors 
recombination, while the associated temperature rise favors 
dissociation. 


Gas Phase Recombination 


It seems possible to eliminate the use of partial ram com- 
pression and the internally supersonic duct as a practical 
solution on the ground of the slowness of the reaction. The 
“residence time’ required to recombine a fraction of the dis- 
sociated oxygen can be calculated from Equations [6 and 8] 
for a given state of the gas at the chamber inlet. It is dif- 
ficult to present the results in general form; one example 
should suffice. This is shown in Fig. 3 for a specific altitude, 
stagnation pressure in the chamber (with an assumed diffuser 
performance as indicated) and with only 50 per cent of the 
total possible recombination taking place. (Note that the total 
possible recombination is in turn less than all the dissociated 
oxygen in the atmosphere, as discussed in the preceding sec- 
tion.) Minimum residence times of the order of 10? sec are 
indicated. But under these conditions only a very small frac- 
tion of the “fuel” is used: At Mach 5 (ao — ar) & 0.08 re- 
sulting in a total enthalpy increase of 65 cal/gr air (Eq. [8]). 
This heat release is not sufficient to produce a net thrust 
(Eq. [15]). Larger fractions of oxygen must be recombined 
and the Mach number must be decreased, both causing the 
required residence time to increase inordinately. These mag- 
nitudes of residence times preclude the use of partial compres- 
sion, because, with the resulting chamber Mach number, ex- 
cessive lengths are needed to achieve a practical degree of re- 
combination. Examples covering the whole range of possibili- 
ties were reviewed, leading always to the same conclusion. 
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Catalytic Recombination 
To shorten the excessive residence time required to com- 


plete the reaction, surface catalysis has been suggested. 
Catalysis cannot shift the equilibrium, but it can accelerate the 
reaction rate. 

Let it be assumed that surface catalysis is possible and that 
through its use the recombination can be affected very rapidly. 
To employ a catalyst, each atom must be made to collide with 
the surface. This is equivalent to bringing the gas to rest at the 
surface, and consequently to stagnation conditions. 

This phenomenon is essential to the understanding of the 
limitation of this engine. It is impossible to accelerate the 
reaction by surface catalysis without raising the temperature 
of the surface and of the adjacent gas layer to stagnation con- 
ditions as well as involving the “self-limiting” character of 
the recombination as the flight Mach number increases.‘ It 
is clear that a catalytic recombination chamber should 
actually be preceded by the best possible diffuser (to obtain 
as high a stagnation pressure as possible). The power cycle 
is then a standard ramjet-cycle with the catalytic surface 
acting as a sort of ‘‘flameholder.”’ — wks & 


Ramjet Performance 


By considering, on the one hand, the residence time for 
volume recombination, and on the other hand, the aero- 
thermodynamic implications of surface catalysis, we demon- 
strated that only a “standard” ramjet where the combustion 
chamber is essentially at stagnation conditions should be 
analyzed further. The very best heat release attainable under 
this condition occurs when all the possible recombination, from 
the composition of the ambient atmosphere to equilibrium, is 
achieved. Combine Equations [7, 8 and 12] in the form 


= Po ) 
(1 — + az) 06 (Pome 


log Kpa — log 


For given ambient parameters, the left-hand side can be 
plotted against (a9 — az). The right-hand side contains the 
product of the ambient pressure and the inlet pressure re- 
covery; the reader can select the value of this term according 
to the degree of his optimism. Figs. 4 and 5 show the solution 
for two characteristic atmospheres. We have designated the 
points corresponding to flight at an altitude of 100 km (Fig. 4) 
and 120 km (Fig. 5). The limits of ram pressure recovery, 
isentropic and normal shock, are shown. Fig. 6 shows the 
same calculation as function of the final (equilibrium) “‘com- 
bustion” temperature. 

It is readily seen that neither the degree of dissociation 
existing in the atmosphere (as long as it exceeds 0.4), nor the 
ambient temperature (nor, for that matter, the exact value of 
the ambient pressure or inlet recovery) bear much influence 
on the heat release. The phenomenon is almost entirely con- 
trolled by the Mach number of flight. 

In further calculations we shall assume an atmosphere at 
To = 240 K; ap = 0.4, and Py = 5.6 X 1077 atm, which cor- 
respond approximately to 100 km in altitude. Within the 
accuracy of our knowledge of the atmosphere, this altitude 
seems to be the lowest before the disappearance of atomic 
oxygen begins to reduce the “heating value” of the fuel. It 


‘This can be avoided only if one postulates some geometric 
characteristics of the catalytic surface such that the particles 
hitting it are all re-emitted downstream. This would take us en- 
tirely out of the field of “thermodynamic” engines; the machine 
would be simply a catalytic screen without diffuser or nozzle, 
traversing a fully free molecule flow. However, it is hard to 
imagine the required directional re-emission properties, and it is 
not even known whether “‘catalysis’’ is possible at all when the 
kinetic energy of the oncoming stream is very large relative to 
the surface. 
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Fig.3 Residence time required to recombine oxygen noncatalyt- 
ically to 50 per cent of the maximum (local equilibrium) possible. 
Stagnation conditions in recombination chamber; flight altitude 

100 km. Po 


= 5.6 X 1077 atm, ao = 0.4, 7) = 240 K) 


= 
= 


is clear that at such low Mach numbers optimum flight per- 
formance will occur at the lowest possible altitude. 

Equation [15] can be used to calculate the performance of 
the ramjet. It is necessary, however, to specify the compo- 
nent efficiencies. The average mean freepath of the atmos- 
pheric particles® at 100 km is around 10cm. The compression 
shock waves of potential supersonic flow can no longer be con- 
sidered discontinuities, but instead, they are shock transitions 
of considerable thickness. At a Mach number of 2, a normal 
shock may be 40 em thick (18), while an oblique shock (the 
thickness of which depends on the Mach number component 
normal to it) may be as much as 100 em thick. It is not clear 
whether a multiple shock external compression inlet can be 
designed at all unless the geometry of the inlet were scaled to 
sufficiently large dimensions, so that the shock transitions re- 
mained comparatively small. 

Moreover, the growth of the boundary layers on an inlet- 
spike (or on the walls of an internal diffuser) at these low 
Reynolds numbers would affect the overall efficiency of the 
inlet severely; the calculations presented below regarding the 
viscous losses in the expansion nozzle can serve as an example 
of this point. The performance of the inlet, where the bound- 
ary layer is subjected to an adverse pressure gradient, would 
suffer even more than that of the nozzle. The aerodynami- 
cist’s experience indicates that even in continuum flow there 
is little overall gain from the use of oblique-shock inlets over 
Pitot inlets at Mach numbers below 2. The latter would also 
be considerably lighter and simpler. Consequently, we feel 
justified in assuming that the overall inlet pressure recovery 


5 Throughout this discussion it is assumed that the gas is in 
kinetic equilibrium and that temperature, mean freepath, etc., 
are defined, at least on an average. The mean freepath is 


L = 1/V/2xno? where o = 3.75 X 10-8 em for air. 
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would not exceed the total pressure recovery across a normal 
shock.® 

If a catalytic surface is used, then we assume that the re- 
combination chamber is very short and neglect the losses, ex- 
cept the shear against the catalyst proper. An efficient 
geometry of the catalytic surface calls for a grid of flat plate 
louvers aligned with the flow. Since each atom must strike 
the plates, but should strike them only once, they should have 
a chord length and also a spacing eoual to one mean freepath 
of the oxygen atoms. (Note that this is longer than the mean 
freepath of the mixture, because of the smaller collision cross 
section of the atom.) The total surface area would equal the 
cross-sectional area of the duct. The shear would be (19) 


[17] 


6 Recent experiments in the low-density tunnel of the Uni- 
versity of California, Berkeley (personal communications) have 
shown an amazingly good diffuser performance at pressures as 
low as 20 4. At Mach numbers around 2, the recovery assumed 
herein, however, is a good approximation. 
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Fig. 4 Maximum possible recombination (to local equilibrium) 
of oxygen at stagnation conditions. Atmospheric data typical of 
100 km 
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Fig. 5 Maximum possible recombination (to local equilibrium) 
of oxygen at stagnation conditions. Atmospheric data typical of 
120 km 
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where p, is the mean upstream density of the gas (all species 
strike the plates), and M, is the Mach number upstream of the 
grid. Setting this shear equal to the pressure drop across the 
grid (the areas are equal), the “burner” pressure ratio is 


=1- Me = 1 — 0.47M, [18 


M, is determined by the normal shock strength at the inlet 
and by subsequent subsonic diffusion. Now, the diffusion of 
the internal flow requires an area expansion which, in turn, 
necessitates increasing the frontal area of the engine over 
that of the inlet opening. The external pressure drag mor: 
than overbalances the decrease in burner pressure loss, par- 
ticularly since diffusion at Reynolds numbers as low as th 
is very inefficient, if at all possible. (By the same token n: 
external pressure recovery over a boattail or in the base regio: 
can be expected.) Therefore, it would seem appropriate t» 
use Equation [18] with M, equal to the Mach number dow: 
steam of the normal shock, that is, assume no subsonic dii- 
fusion. The engine would then be a cylinder of uniform cros 
section (with walls as thin as possible). 

The nozzle efficiency depends on two factors: The complete- 
ness and uniformity of the expansion at the exit of the nozzle 
and the growth of viscous boundary layers on the wall. Th: 
first consideration can be related to the ratio of nozzle lengt! 
to its diameter; the second can be estimated from know 
boundary layer theories once the geometry and the pressur:: 
gradients are known. That the boundary layer losses ar: 
substantial in a rarefied gas expansion can easily be verifiec! 
(20). 

Fig. 7 shows the internal performance of the oxygen recom- 
bination ramjet at a typical altitude based on the calculation: 
and the values for yz and 7» discussed thus far, with 7, as a 
parameter. The calculation of 7, and also of the external 
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Fig. 6 Final combustion temperature for recombination to 
local equilibrium (same conditions as Fig. 4) 
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drag of the engine vehicle is less precise and will therefore be 
considered separately as follows. 

The supersonic portion of the nozzle must be sufficiently 
long to contain the expansion of the flow. Experience indi- 
cates that a 15 deg divergence half-angle of the supersonic 
nozzle is a practical maximum consistent with good nozzle 
performance. Since we assumed a normal shock inlet without 
subsonic diffusion and a short recombination chamber (order 
of one mean freepath) the nozzle length will essentially charac- 
terize the engine length. The combustion chamber area, the 
nozzle throat and exhaust areas vary only little (in the flight 
regime of practical interest), so that for large diameters, the 
engine dimension can be accurately characterized by a single 
diameter D, say, the diameter of the nozzle exhaust. We shall 
assume that the overall engine length is 50 per cent larger than 
the divergent portion of the exhaust nozzle. Thus 


L , 1 — 
D sin 15 deg 


where A*/A, is known from the thermodynamic cycle as a 


[19] 


temperature, which is in turn related to the flight parameters, 


function of altitude, flight Mach number (and io, mere v2 


it has the value of 1.4. For favorable pressure gradient it is — 
less, but not much below unity. (R/L) is the Reynolds num- 
ber per unit length. 
In exactly the same manner, the external drag can be 
written 
Ty 
DR. P\M,? (L/D) (R/L).: “wp” 
The drag due to lift for a gross weight of the vhicle (Y lb) is 
expressed in terms of the lifting efficiency € (the overall lift-to- 


drag ratio). 
Equating thrust to drag and solving for gross weight one 


finds 
F/A 
Po 4 


F/A\w 
p2 


As before, we shall illustrate this for the specific altitude of 
100 km, but, in this generalized form the results are not too 
sensitive to altitude (provided a remains above 0.3). 

The external specific Reynolds number is then 


[24] 


(R/L)ex = 5M, 


and the internal specific Reynolds number (for viscosity 
varying linearly with temperature) is 


y-1 3/9 
1 — M,? 
( ae ) 


(R/L)ez Tr My 
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Mo 
Fig. 7 Specific impulse and thrust per cross-sectional area of 
oxygen recombination ramjet. Normal shock diffuser and drag 
loss against catalytic surface included. 100 km altitude (re- 
duced parameters reasonably insensitive to altitude provided 
ao 0.3) 
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CATALYST STRIPS 


Fig. 5). at 
The net thrust of the engine is » oaeh 
1000;- 
~\AJ, 4 D = 
where (F'/A), is the (isentropic) thrust/exhaust area, and 6 a z 
is the momentum thickness of the internal boundary layer. > 6 = 
Although both the internal and the external flows are in the > ; =e 
“slip flow” regime (20), the growth of the boundary layers is i 
well approximated by the usual dependence on Reynolds a a, 
number 
6 (L/D)'/ PAYLOAD © LB 
- = — = [21] SKIN 0.00! IN. THICK 
R, (R/L)n Mo *2 ALT. 100 KM. 
The subscript (m) indicates internal flow in contrast to the — | 
external flow (subscript ex); K can be evaluated from 
laminar boundary layer analyses. For zero pressure gradient Mo 
a b 


Fig. 8 Flight performance of ‘‘minimum’’ oxygen recombina- 

tion ramjet: (a) Gross weight capability vs. Mach number at 

100 km altitude; (b) minimum engine size for zero payload 

cruise at 100 km and M, = 2 as function of density of the struc- 

tural material. Assumed configuration shown in insert to scale. 
No wings or jet deflectors included 


The values of K, and K,, are set equal to 1, which cannot 
cause a large error (and on the contrary, probably under- 
estimates the external drag). € is also taken to be unity. It 
is doubtful that aerodynamic lift-to-drag ratios higher than 


[23] 


this can be obtained in rarefied gas flows. (For free molecule 
flow lift-to-drag ratios around 0.2 can be expected.) Further- 
more, the wing loading would be so small (around 500 ft? per 
lb of lift at Mop = 2 and 100 km altitude with a 30 ft chord) 
that the structural weight and practical design problems would 
probably rule out the use of wings in favor of a direct thrust 
vector lift, which is what we suspect would be most efficient. 
In that case the best € would be 1. 

The results are shown in Fig. 8. As we implied, the best 
performance occurs around My) = 2, and it is not very spec- 
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The engine weight has not yet been taken into 


tacular. 
account. If the skin is assumed to be 0.001 in. thick and only 
the external skin plus the catalytic surface are considered, an 
estimate of the net payload capacity of the vehicle can be 
made as indicated on the figure.’ An all aluminum engine 
could not fly at all, regardless of size; if built partly of 
aluminum foil and partly of “mylar’’ plastic, only incredibly 


large engines would sustain themselves. A mylar engine is 
also shown, although mylar melts at about 800 K, a consider- 
ably lower temperature than required. 


Conclusion 


The oxygen recombination ramjet is limited to flight below 
Mach number 4 and has its optimum performance around 
Mach number 2. Internally supersonic recombination is im- 
practical because of the excessive recombination times re- 
quired to complete the reaction, whereas catalysis necessarily 
implies stagnation conditions in the reaction zone. 

Optimistic but reasonable estimates show that only very 
large and extremely light engines can sustain flight at all. 
Even then, their payload capacity is extremely small. Opera- 
tional problems, primarily the launching and “unfolding” of 
the engine and then boosting it to speed, would seem to be ex- 
tremely severe. The statements that the engine would be 
ready to fly in two years (21) appear to be somewhat prema- 
ture. 
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Post Scriptum 


Shortly after this article was submitted for publication, 
another analysis of the oxygen recombination ramjet by 
Baldwin and Blackshear (22) became available. This study 
is concerned only with volume recombination (no surface 
catalysis) of the ambient atmosphere in a duct which con- 
sists of a frozen flow diffuser followed by an equilibrium flow 
expansion nozzle. The engine is assumed to fly at orbital 
speeds with only partial compression, i.e., the internal flow 
is totally supersonic with velocities in the order of 20,000 
ft per sec. 

Baldwin and Blackshear study an extremely idealized 
system. In a sense, they seek a value of the recombination 
rate constant k, that could make the engine feasible. In the 
text they do show that even for k, as high as 10"” (cc/mole)? 
(one order of magnitude higher than the most recent data of 
(9)) the engine is no longer feasible. A careful reading of their 
analysis seems to complement our findings that the orbiting 
recombination ramjet is not marginal but impossible in view 
of the best current knowledge of the reaction kinetics, even 
without assuming reasonable aerodynamic losses. 


Nomenclature 

A = cross-sectional area, ft?, Ag = nozzle exit area a 
D = diameter of nozzle exit, ft a 

“om = average specific heat of mixture, cal/gr deg K a 
F = net thrust, lb _ 
h = specific enthalpy, cal/gr 

AH = enthalpy of dissociation of molecular oxygen eo 
I; = specific impulse, sec 


7 When the engine size exceeds about 100 ft in diam., the hoop- 
stress begins to require a thicker skin. 
8 Product of Dupont Company. a4 vee 


k = reaction rate constant; k, recombination rate constant 
for oxygen, (mole/cc)? 

Kya = equilibrium constant for oxygen dissociation based on 
partial pressures, atm 

L = characteristic length of engine, ft 

if = pressure, atm 

R = Reynolds number, also universal ideal gas constant in 
atm-cc (mole) K) =! 

T = absolute temperature, deg K 

t 

q 


= time, sec 
= heat release of recombination of oxygen, gal/gr of mix- 
ture 
a = degree of dissociation of oxygen 
Y = average isentropic exponent for frozen flow 
p = density of mixture gr/cc 
n = component efficiencies: nz ram recovery, ng recombin:i- 
tion chamber pressure ratio, 7, nozzle efficiency 
Subscripts 
0 = ambient atmospheric conditions 
8 = stagnation conditions 
E = chemical equilibrium at local thermodynamic condi- 
tions 
1 = conditions at inlet to isobaric recombination chambe 
F = conditions at end of recombination chamber 
c = recombination chamber ad 
n = internal surfaces in expansion nozzle os: 
ex = external surfaces 


Square brackets denote concentrations, moles/ec_ 
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HE CHARACTERISTICS of an optimum aerodynamic 

configuration for a beam rider missile, which obtains al- 
most all of its lift directly from deflection of the control sur- 
faces, are discussed. The performance improvements over 
conventional missile configurations identified by a three- 
dimensional simulator study are as follows: 

1 Improved beam riding characteristics and increased 
guidance accuracy in the presence of large beam rotation rates 
and beam jitter. 

2 Improved beam riding performance at high altitudes. 

3 Large reduction in induced rolling moments. 

4 Elimination of requirements for artificial pitch damping. 

This paper has been prepared to present the results of a 
study of an optimum aerodynamic configuration for a beam 
riding missile. The study was originally concerned with the 
aerodynamic changes required to extend the high altitude 
performance of a beam rider missile. 

A beam rider guidance system is defined as a system for 
guiding missiles which utilizes a beam directed into space, such 
that the center of the beam axis forms a line along which it is 
desired to direct a missile. The beam, which may be fixed or 
moving in elevation and azimuth, may be a radar beam or a 
light beam. Equipment is contained in the missile that can 
determine when the missile is in the center of the beam, or can 
determine the direction and magnitude of the error when the 
missile has deviated from the center of the beam. This equip- 
ment involves suitable electronic circuits, servo mechanisms, 
aerodynamic control surfaces and other equipment such that 
the missile will automatically operate to reduce the error dis- 
placements from beam center to zero. 
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Selection of an Aerodynamic 
_ iguration for Improved Beam 
Rider Guidance 


rates, and radar noise during the remainder of flight. 
configuration which performs well at low altitude tends to have inadequate response at high altitude. 
The necessary improvement in response can be achieved either by an aerodynamic configuration 
which provides a very high missile pitching frequency or alternatively by a configuration which 
obtains all the lift by wing deflection alone. 
of attack) missile, is particularly interesting because it eliminates the dependence of missile re- 
This configuration yields substantial improvements in beam riding perform- 
As secondary advantages, system simplification is achievable through the 
elimination of the need for rate gyro feedback to damp body osc ‘illations and through the consider- 7 
able reduction of induced rolling moments. 
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The beam rider type air-to-air missile system basically requires a missile with a fast response to 
This characteristic is necessary for good synchronization with the beam 
at the beginning of missile flight and for good beam riding, in the presence of large beam rotation 


A missile with a conventional wing control 


The latter configuration, essentially a no-a@ (no angle 


In air-to-air applications, the beam, generally a radar beam, 
is directed from the launching aire raft toward the target air- 
craft. The beam must contain provisions for continuously 
tracking the target in both elevation and azimuth, regardless 
of any relative motion between the launching aircraft and the 
target aircraft. Further, the beam must contain coordinate 
information to allow the missile to interpret its error, and the 
sense or direction of error from beam center. This character- 
istic can be achieved in a number of ways, one of which is 
conical nutation of the beam center, to provide error mag- 
nitude, combined with pulse coding of the beam during the 
nutation cycle to provide direction sensing. 

Beam rider missiles are launched as nearly as possible along 
the axis of the beam. This procedure generally is required to 
allow the missile to enter the beam at a sufficiently small 
crossing angle to allow capture of the missile by the beam. If 
this angle is relatively large, the missile will experience a large 
transient in its flight path in its attempt to ride the beam, or 
may pass through the beam completely, and be lost to it. 

For air-to-air missile systems, the flight path of the launch- 
ing aircraft ideally should be one which produces a minimum 
of lateral motion of the beam and, consequently, the missile. 
The requirement for launching the missile along the axis of 
the beam, combined with the tactical objective of insuring 
attacks off the target’s nose or tail, produces radar beam rota- 
tion rates. These beam rates introduce an environment of 
continuous acceleration demand within which the beam rider 
missile must perform. 

In practice, radars do not track complex targets smoothly, 
but jitter their position in space among the number of reflect- 
ing elements contained within the complex target. This jitter 
contributes to the noise of the system. For airborne radars, 
requiring radar beams space-stabilized and slaved to the tar- 
get aircraft, it is impossible to filter out all the jitter without 
affecting the radar’s tracking capability. This radar noise is 
an environment in which the beam rider missile must perform. 
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missile with a fast acceleration response to acceleration com- 


The beam rider type air-to-air missile basically requires a 


mands. This characteristic is necessary for good synchroniza- 
tion with the beam at the beginning of missile flight and for 
good beam riding (during the remainder of flight) in the pres- 
ence of beam rotation rates and beam jitter, which are sig- 
nificant compared to the accuracy required. A missile which 
develops a large fraction of its lift from angle of attack per- 
forms well at low altitude, but tends to have inadequate re- 
sponse at high altitude. The necessary improvement in re- 
sponse can be achieved either with an aerodynamic configura- 
tion which provides high missile pitching frequency at high 
altitude or, alternatively, by a configuration which develops 
lift independently of the missile angle of attack response. The 
latter can be achieved witha configuration which obtainsall the 
lift from wing deflection alone, essentially a no-a missile. The 
use of this configuration has been investigated as part of an 
improvement program in connection with a _ bes am riding 
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Fig. 2 System stability characteristics—conventional missile 


System Control Considerations | 


The beam rider missile control loop requires a relatively 
high bandwidth in the forward part of the loop from error 
command to missile acceleration. This requirement lies in 
the nature of the beam rider. It must fly along one particular 
path in space—the beam—in order to complete its mission 
successfully. Since missile control surface motion produces an 
acceleration of the missile with respect to the path, the open 
loop transfer function of the missile basically is a 40 db/decade 
system. Therefore, it is necessary to provide derivative com- 
pensation in the missile to obtain a stable loop, and the result 
is the large missile bandwidth. If there are delays in the loop, 
such as may be caused by the control surface actuator or aero- 
dynamic response, they either have to be compensated or the 
design must allow several octaves of frequency between the 
critical control loop frequencies and the frequency where the 
delays have effect. 

Thus, the open loop transfer function of the beam rider 
missile may be written as 


20 T 80° € s? 
AERODYNAMIC TRANSFER FUNCTION 
Lm, —|MAG 
-10 > -40° € = missile to beam error 
PHASE \ = beam rider loop gain 
-20 | fs) = compensation frequency functions 
> -30 T wo 3 f(s) = all other frequency functions including the missile con- 
2 4 trols and aerodynamics 
ne s complex frequency variable 
> 20 ‘eS er -100° z The gain X is chosen to provide adequate beam rider loop re- 
< 10 ee eee eer eae — sponse and tolerable steady-state error when the missile ex- 
=] as ee PHASE periences a steady acceleration demand. These considerations 
oe -1s0° establish the zero gain crossover frequency of the open lcop 
vA N g. function. In order to have a well damped system, it is neces- 
=e MAG. ie sary to provide a phase margin on the order of 30 to 40 deg at 
220 N _— this frequency. The function f.(s) is chosen to supply this 
* margin without producing major changes in the zero gain 
- 300° - crossover frequency. It is not practical to compensate for 
oz os seas 2 5 ad ‘wel the additional delays f,(s), so they must be designed to have 
ATIVE effect at least several octaves above the zero gain crossover 
Fig. 1 System stability characteristics conventional missile frequency. 
ae low altitude POET retina Sy The primary delay is the aerodynamic response, a typical 
a4 plot of which for a conventional missile at low altitude is 
shown in Fig. 1. The corresponding beam rider open loop 
20 ‘ soe frequency plot, including the aerodynamics for low altitude 
AERODYNAMIC at FUNCTION a conditions and compensation for stability, is shown in the 
10 1 | MAG. 40° same figure. There is adequate phase margin for stability at 
\ this altitude. A similar plot using high altitude aerodynamics 
laa ‘2 BS is shown in Fig. 2. Here the phase margin is practically non- 
-10 \ + -40° existent. Several solutions are possible to improve the phase 
2 ‘ ] VA margin. One of these solutions is to lower the beam rider 
a \ _ loop gain and shift the compensation to lower frequencies. 
\iy -120° However, the resulting performance is not acceptable, par- 
PHASE ticularly that performance associated with missile lag behind 
3 0 0? a rotating radar beam. Another solution, and perhaps a more 
desirable one, is to modify the missile aerodynamics. This 
= MISSILE OPEN LOOP TRANSFER FUNCTION a modification can be done in two ways: Increase the aerody- 
< 10 f namic natural frequency of the missile, or eliminate the angle 
= Local / \ of attack response of the missile, and consequently the delay 
SLPHASE | 7 associated with this response. 
Any significant increase in aerodynamic frequency would 
tt me require a large increase in wing and tail area, since the aero- 
dynamic frequency is proportional to the square of these 
a ia vane) areas. Not only would size alone prohibit such a change, but 
2 5 10 - 5 the increase in drag and possible increase in induced roll 


would make it intolerable. The only acceptable solution then 
appears to be the elimination of the angle of attack response of 
the missile. The aerodynamic configuration presented in the 
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next section essentially eliminates the delay by obtaining 
nearly all lift from the movable surfaces. In addition it offers 
other advantages. Since the angle of attack is eliminated, 
induced roll is reduced to a very small percentage of that ex- 
perienced by a conventional configuration. The need for rate 
gyro feedback for damping of body angular motion is eliminated. 
With the aerodynamic delay eliminated, the bandwidth of the 
missile can be reduced by subjecting the error signal to addi- 
tional filtering. The decrease in bandwidth reduces the usage 
rate of control surface actuator fluid and the effect of noise in 
the presence of system nonlinearities. In addition the missile is 
less susceptible to noise caused by jamming or other inter- 
ference. 

The negligible angles of attack resulting from wing de- 
flection eliminate the need for body attitude stabilization. 
Since the body attitude tends to follow the relative wind in a 
turning maneuver rather than pitching to develop lift from 
angle of attack, feedback of body pitch rate to wing deflection 
would have relatively little effect on the missile flight path. 
Conversely, wing deflection also has relatively little effect on 
body attitude. Thus, body attitude stabilization is obtained 
primarily by the aerodynamic stability and damping provided 
by the missile tail. The static aerodynamic stability for the 
no-@ configuration described below is estimated to be slightly 
greater than for the conventional configuration. This level of 
stability is ample for body attitude stabilization, as was 


demonstrated by the analog computer results. 


Missile Aerodynamics = 


The basic characteristic of the no-a@ system is the utilization 
of an aerodynamic configuration in which nearly all lift is ob- 
tained directly from deflection of the control surfaces. Maxi- 
mum lateral acceleration is obtained immediately upon de- 
flection of the wings, thereby eliminating the delay associated 
with development of lift due to angle of attack. This in- 
crease in speed of response is shown in Fig. 3, which compares 
the lateral acceleration responses of a conventional missile and 
a no-@ missile to a step wing deflection. 

The appearance of a configuration which achieves a fairly 
good approximation to the no-a concept is not greatly different 
from a more conventional typical missile configuration, as 
shown on Fig. 4. Essentially, the desired aerodynamic 
characteristics are obtained with wing and tail planforms 
whose geometry is chosen with proper attention to the down- 
wash interference between wing and tail. 

The practical realization of a no-a missile depends on some 
means of fixing the center of pressure at the center of gravity 
throughout the Mach number range. 

An investigation of test data, summarized on Fig. 5, shows 
the center of pressure shift over a limited Mach number range 
as a function of the ratio between wing and tail span for a 
specific set of planforms. In accordance with theory, the 
farther the tail is removed from cores of the wing trailing 
vortices, the less the wing-tail interference and the less the 
shift of center of pressure with Mach number. 

For the conventional configuration, with almost equal 
spans, the shift is about four body diameters, whereas for a 
configuration with a span ratio of 1.4, the shift is reduced to 
about one body diameter. Further increase in the span ratio 
to 2.4 reduces the shift to negligible values. 

The shift of center of pressure for the chosen no-a@ configura- 
tion is about one-half body diameter over the Mach number 
range from 1.2 to 3.5. Although this is a great reduction, the 
reduction in stability resulting from the relatively smaller 
tail yields some trim angle of attack at the lower Mach 
numbers. 

A direct result of the reduction of trim angle of attack of the 
no-@ missile compared to conventional configuration is a 
loss of body and tail lift. This is offset by the reduction in 
negative lift induced on the tail by wing deflection in the no-@ 
configuration, and in increase in the no-@ wing area in order to 
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equal the trim lift capabilities of the conventional configura- aerodynamic characteristics of the conventional and no-a 
tion. The lift due to wing deflection is from 70 to 90 per cent missiles are shown in Fig. 6. 

of the total lift for the no-a missile, compared to 40 to 60 per An additional benefit obtained from the revised wing and 
cent for the conventional configuration. Comparisons of the tail planforms is a substantial reduction in rolling moment in- 


by the lifting surfaces on missiles experiencing sideslip ir 
eh combination with pitch accelerations. This effect is also in 
agreement with theoretical considerations concerning induce: 


rolling moments on such bodies. Further, since the no-@ mis- 
sile essentially eliminates angles in either yaw or pitch, the in- 
duced moments are attenuated almost to the point of non- 
existence. 

A study of the hinge moments of the no-a wing indicate: 
that, with proper choice of hinge line, the supersonic hing: 
moments are of the same order of magnitude as for the con- 
ventional configuration. 

Missile zero lift drag is essentially the same for both the no 
a and the conventional configurations. However, drag due t: 


CONVENTIONAL —4 
CONFIGURATION 


“NO CONFIGURATION 


T 
“NO @" CONFIGURATION 


rede, lift is larger for the no-@ configuration because larger wing de 
LOAD CONFIGURATION flections are required for equivalent lateral load factor; thu 
a ed the no-@ missile has somewhat more drag deceleration during 
4 7 gliding flight than does the conventional configuration. 
res 
a 1 2 3 System Performance 
MACH NO. 


A simulator study was conducted to evaluate beam riding 
performance with the no-a@ configuration. The study was pari 
j= < of a much broader investigation where every effort had bee1 


Fig. 6 Comparison of aerodynamic characteristics of conven- 
tional and no-a missile 


— made to simulate missile dynamics and contro] system char- 
om LAUNCHING AIRCRAFT TURN-O5g acteristics as accurately as possible. This completeness in- 
Es} ALTITUDE - 50,000 FT. cluded simulating missile motion with six degrees of freedom 
° 
gat The simulator results showed a large degree of improve- 
<A ca . . 
et 3 ment in performance with the no-a@ configuration over a sys- 
g ys 
oll, tem with the conventional configuration. The improvement is 
12 demonstrated by Fig. 7 which shows the decrease in varianc 
of missile to average beam position error caused by beam 
Weer fe a jitter. There also was a significant decrease in the syn- 
— est Sater ie chronization times at the start of missile flight and a sharp 
a reduction in the amount of missile roll. The perturbations in 
= =< | angle of attack were quite small, so that nothing was sacrificed 
ag by eliminating the rate gyros. Not only was performance im- 
ok proved at the altitude represented by Fig. 7, but the operating 
° 6 8 10 l2 range was significantly extended. 


RANGE- THOUSAND FEET 


Conclusion 


The ability to eliminate the major detrimental aerodynamic 
effects from the beam rider missile control loop offers the 


LAUNCHING AIRCRAFT TURN-Q5g opportunity of obtaining a large improvement in performance 

ALTITUDE-50,000FT| | | | | and a reduction in missile system complexity. Being able to 

6 8 10 2 4 accomplish all filtering in the error coordinate system allows 

RANGE- THOUSAND FEET a maximum choice of control function frequency shaping, 

Fig. 7 Comparison of system performance of conventional and thereby minimizing noise problems and control energy re- 
no-a configurations quirements. 
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of the dynamic pressure phe- 
nomenon within the thrust chamber of a rocket motor, 
cperating in an unstable manner, is a difficult, but necessary 
task in the research of rocket combustion instability. The 
transducers used in this application must be capable of opera- 
tion in a very extreme vibration and temperature environment. 
In addition, these transducers must also accurately provide an 
analog signal of the pressure oscillations. 

Several well-known pressure transducers, such as Photocon, 
Norwood Controls Model 107 and the early Li Liu (now known — 
us the Norwood Controls Model 107) were evaluated through — 
the use of a shock tube to obtain the transient response to 
step pressure excitation. The resultant transient response is 
analyzed to obtain characteristics, such as rise time, per cent 
overshoot and per cent damping. Further, the frequency re- 
sponse was obtained by the use of the ‘“Guillemin Impulse 
Approximation Technique” and digital computer assistance. 
A discussion of the analysis inherent approximation errors and 
concluding comparisons of the characteristics thus obtained 
for each transducer is given. 


General Testing Procedure 


A shock tube was used to generate a step pressure excitation 
to the input of the pressure transducer. The input step gen- 
erated has a rise time of less than 10 microsee with a duration 
of 10 millisee. The pressure step magnitude can be varied by 
adjusting the initial chamber pressures and diaphragm ma- 
terial over a range of 0.1 to over 500 psi. Because a number of 
excellent papers have been written covering operational and 
theoretical problems in the use of shock tubes, this paper will 
not include this discussion. 

The transducer response signal is displayed on an oscillo- 
scope and recorded by a Land camera. A typical transient re- 


sponse is shown in Fig. 2 — 


Direct Analysis* 


Considerable information is available directly from the 
transient response. In general, transducers capable of fre- 
quency response to the range of 10 ke have by the nature of 
the construction of the sensing element very little inherent 
damping. Thus, the typical response is similar to that of the 
classical second order system with a damping ratio of 10 per 
cent or “ps 

Fig. 2 indicates the method of obtaining the following 
characteristics by direct measurement: 

Damping Ratio 


62 A, — As 
| D.R. = loge {2 ac = 6 


Received May 5, 1958. 

1 Instrumentation Development Engineer, Facilities Engineer- 
ing Department. 

* A list of definitions of symbols appears at end of paper. 
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Transducer Frequency 
Response Evaluation for Rocket 
Instability Research 


Azusa, Calif. 


CHARLES R. TALLMAN! 


Aerojet-General Corp. 


For small damping D.R. & 6/2r. 
Natural frequency of the sensing diaphragm or suspension 
system 


Wo 
— D.R2)” 


21 
< 
= 
12 


Wn = 


For small damping Wn = Wo 
Per cent overshoot, Mp 


Rise Time = time required to reach the 10 to 90 per cent 
steady-state; Tr = 7}. 

If the transient response has one dominant ring frequency, 
it can be treated as a second order system with good approxi- 
mation and the following relationships hold: 

Time to peak, Tp = 1/Wo. 
Per cent overshoot, Mp = e~?7?; o = WnD.R. 
Decay envelope on ring frequency has the form f(t) = e~”. 


Fig. 1 Shock tube facility used for transducer frequency 


response evaluation at Aerojet-General Corp. 


Ay x 100 


%OVERSHOET: 


SCRITICAL DAMPING= Log. 


Ring FREQUENCY Wo, 


RING FREQUENCY (Jo. 


Treicat Response Curve 
to the step input pressure magnitude Ay 


Fig. 2. A typical pressure transducer’s response to a step pres- 
sure excitation showing the measurements that can be made 
directly 
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TRANSDUCER 
MPLITUDE 
E (s) G(s) =— Cs) STEP EXCITATION 0 


Fig. 3 Four terminal block diagram in- 
dicating the input and output terminals 


System Function 


The transient response to a step input of the transducer can 
be analyzed to give the complete description of the trans- 
ducer’s system function. 

The system function is equal to the ratio of the output/in- 
put as a function of the complex frequency (s). 

Thus G(s) = C(s)/E(s) for a linear system, where E(s) = 
Le(t), and C(s) = Lc(t). The input pressure excitation 
generated by the shock tube is then e(t) and is idealized as a 
perfect step. Thus the Laplace transform of the step excita- 
tion is simply 


Frequency Response 


The frequency response of the transducer is a plot of the 
magnitude of the system function as a function of frequency 


where s = jw. 

Since G(s) = C(s)/E(s), and E(s) = (1/s)(Pi — Po), then 
G(s) = sC(s)/(P: — Po), and frequency response = |G(jw) |. 

The technique for obtaining the value of |G(jw)| from the 
transient response data is best described by an example. The 
example to be used is of practical worth, as will be shown later, 
for quick estimations of the frequency response. 

Assume the output C(t) for a pressure transducer is simply a 
ramp type function as shown in Fig. 5. The derivative of this 
output function is simply C’(t). The second derivative C’t is 
also directly obtained. 


{1 — 


Ts 


= Gs) | G(jw) | = 


9 
: “tite result is useful in itself for quick estimation of the fre- 


Cyst 
=O Ttgeo 


Chy= Seq - 
-+ 5(e-1) 


Fig.6 First derivative of the 
ramp response function C(t) 
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Fig. 7 Second derivative of 
ramp response function C(t) 


Fig. 4 Idealized pressure step excitation 
to the transducer input generated by the 


1 2 
Bn cos wr.) + ( Bn sin WT», 
w n=0 n=0 


Fig. 5 Example response function of a 
system to an ideal step excitation 


quency response of a transducer, if the response approaches 
the ramp function. Similar results are reported in (1).? 

The Guillemin’s Impulse Method of Approximation (2) de- 
scribed by the example above can be applied to any oscillator) 
transient response by approximating the oscillation by) 
straight line segments. 

If the transient response g(t) is given, and the correspond- 
ing transform G(jw) is desired, the following steps are used: 

1 G(é) is approximated by a sequence of straight lines, say 
gx(t). 

2 The approximated function g,(t) is differentiated 
enough times (n) to make g,”(¢) a sequence of impulses. For 
straight line approximation two differentiations are sufficient 
to make g+"(t) an impulse chain. 

3 The transform of g,(t) then can be written by inspection 
as a sum of p terms, where p is the number of break points in 
the straight line segments. 

The generalized solution for the equiv alent frequency re- 
sponse for the transient response C,(t) is then 


| G(jw) | = 


where 


—An 
tn4i — tn 


)-( 


Ay Art) 
tn tn-1 


A,, terms are defined as the magnitudes in the generalized 
transient response, shown in Fig. 8. 

Application of this technique to a very oscillatory type of 
transient response involves approximations in the time do- 
main and, of necessity, some error is introduced. To evaluate 


1 
1! 
oa 
ttt TIME 
Fig.8 A generalized transducer’s response to a step excitation 


showing the straight line approximation. The sampling times 7. 
and the corresponding amplitudes 4. 
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Fig. 9 A second order system response to a step excitation 
showing the straight line approximation C, (¢) 


the magnitude of error introduced in sampling the transient, a 
second order system with a transfer function of 


1 


( s \? DR. 
Wn “Wn 


natural frequency of the transducer’s sensing 


+1 


where Wn = 
element. 
For the mechanical system F(t) = M# + Cpz + Ker 


= 


The transient response to a step input of the second order 


system is 
Wn Wo 
— e-Wn D.R.t Wot — 
Error Analysis 
Set Wn =m X 10% rad/sec; D.R. = 1/x, and evaluate 


the C(t) above for various values of t; the resulting transient 
response is shown graphically in Fig. 9. 

If this response is approximated by straight line segments as 
shown, where the amplitude values are obtained by evaluating 
C’(t) = 0 for each maximum and minimum, sampling then 
occurs at a period of 

The resulting frequency response obtained for the ap- 
proximated transient C,(t) then deviates from the actual fre- 
quency response as shown in the Fig. 10 where 


1 
2D.R.[1 — D.R.2]'/2 


Mr = = 1.51 


Mrs (approximated) = 1.40. Per cent error, 7.8 
Wr = Wn[1l — D.R.2]'/? = 4475 eps 
Wr« (approximated) = 4200 cps. Per cent error, 6.27 


With further sampling of the given transient response, the 
error is reduced as shown in Fig. 11. 

The foregoing analysis was then applied to the three trans- 
ducers mentioned, the Photocon Model 342, Norwood Control 
TY and the Norwood 107, all of 2000 psig range with transi- 
ent responses as shown in Fig. 12. The equivalent fre- 
quency then obtained is in Fig. 13. 

The dynamic characteristics for the transducers are shown 
in Table 1. 

In conclusion, a recently developed mathematical tech- 
nique has been utilized for determining the frequency re- 
sponse from a given transient response. This technique, when 
using the most simple approximation of straight line segments 
at a sampling rate equal to the 3 period of the ring frequency 
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T 
F (w) | 
w | 
O3- SAMPLE FREQ | 
} 
| 
“a 1000 We 10000 


FREQUENCY, cps 1 


Fig. 10 The calculated frequency response for the example 
second order system showing the difference obtained when addi- 
tional sampling is made in the step time response. As the 


sampling is increased the calculated curve approaches the 
theoretical values 


° 4 1 
Wo 2Wo 
TT 


SAMPLING FREQUENCY 


Fig. 11 The effect of sampling in the time transient response is 
shown. The error in the calculated frequency response com- 
pared to the theoretical is plotted vs. the sampling frequency 


Table 1 Computed dynamic characteristics for three 
a high frequency response transducers 
Freq. flat 
+10 
Tp Tr Mp Wn Mr per cent 
sec sec percent keps keps 
Norwood 
ETY 33 10 =110 30 11 2.3 
Norwood 
107 20 6 36 35 2 10 
Photocon 
342 18.2 40 28.6 55 10 4 


in the transient, results in an error of less than 10 per cent in 
the amplitude of the amplification in the frequency response. 
Further sampling at a smaller period of }7> results in errors 
not greater than 5 per cent. 
The results of this analysis applied to the transducers used 
in rocket instability study indicate that some correction might 
be required to the rocket oscillatory data if it exists above 3 
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500 usec. 


Fig. 12 The actual response of three transducers to a step pres- 
sure excitation of approximately 350 psi. 
range transducers. 


These are 2000 psi 
The amplitude was adjusted for each to 


improve the readability. 


ke, and this correction can be obtained from the a 
response curves given. 


Nomenclature 


damping ratio to critical damping, unitless 


D.R. = 

6 = logarithmic decrement, unitless 

Wn = undamped natural frequency of a resonant system, 
rad/sec 

Wo = damped natural frequency of a resonant system, 
rad/sec 

To = period of oscillation of a damped resonant system, sec 

Mp = per cent overshoot, unitless 

Tr = rise time from the 10 to 90 per cent value, sec 

Tp = time to the first peak in the response to a step input of 
a resonant system 

o = damping exponential constant, sec 

E(s) = LaPlace transform of the input signal 

e(t) = input excitation signal to the device under investigation 

g(t) = time equivalent of the system function 

G(s) = system function of the device being investigated 


FREQUENCY 
RESPONSE 


T 
wn 


RELATIVE AMPLITUDE 


4 


NORWOOD ETY 
NORWOOD 107 — 


10 80 
FREQUENCY (KC) 


Fig. 13 The computed frequency response from the transien: 
responses for three transducers 


C(s) = LaPlace transform of the output of the system or de 
vice being investigated 
c(t) = output signal of the device being investigated 
= pressure in the shock tube, Py ambient expansion cham 
ber pressure, P; shock step pressure 
6(0) = dirac delta function at time 0 
P = number of points taken in sampling of step response 
s = complex frequency = o + jw 
Mr = amplitude ratio (maximum value) 
Wr = frequency at which Mr occurs in frequency response 
M = effective mass of the sensing element 
Kk = spring constant of the mechanical system 
Ca = damping constant 
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Laminar Heat Transfer on 


Three-Dimensional Blunt Nosed Bodies 


transversal pressure gradients. 


derived. 


URING the last few years several theoretical investiga- 

tions of three-dimensional boundary layers have been 
published; many of these have been reviewed in (1).’ Only 
more recently, quantitative results have been presented for 
some particular problems as the stagnation line flow over a 
yawed infinite cylinder (2), the windward streamline in the 
plane of symmetry of a yawed cone (3), and the general three- 
dimensional stagnation point in the presence of a very cold 
wall (4). Apart from these results, very little information is 
available that may lead to reasonable engineering estimates 
of the boundary layer characteristics over more general 
bodies. 

Of particular current interest is the determination of the 
heat transfer on axisymmetric blunt nosed bodies at incidence 
and on sweptback wings with blunt edges in hypersonic flow. 
For these problems one may contemplate the possibility of 
extending methods previously applied in the analysis of 
three-dimensional boundary layers, namely: 

1 Methods based on linearization of the crossflow. 

2 Methods that consider the nonlinearity of the three- 

dimensional viscous flow, but are limited to particular situa- 
tions (e.g., the “initial” boundary layers of (2 and 3)). 
For example, in the case of a blunt nosed body one may cor- 
respondingly investigate the small incidence problem for some 
categories of body shapes (5), and the viscous flow on the 
windward side of the pitch plane for shapes not included in 
the aforementioned categories and for large angles of attack. 

In general, a three-dimensional boundary layer can be rep- 
resented by perturbations of the pertinent two-dimensional 
ir axisymmetric solution under the following conditions: 

1 The three-dimensional effects in the inviscid flow can 
linearized. 

2 The curvature of the outer edge streamline has a small 
‘omponent in the plane tangent to the surface of the body (6). 
Both requirements are satisfied in many instances of flow 
ibout rounded noses at small incidence. The perturbation 
inalysis may then be applied, provided the results are used in 

Received May 27, 1958. 

' This research was conducted under the sponsorship of con- 
ract AF 33(616)-3265 monitored by the Aeronautical Re- 
search Laboratory, Wright Air Development Center, USAF. 

* Associate Professor, Department of Aeronautical Engineering 
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transfer over the nose portion of arbitrary three-dimensional blunt bodies in hypersonic flight is 
The approximate results thus obtained compare very favorably with known exact solutions. 
The method in question reduces the boundary layer problem to the determination of the inviscid 
streamline pattern on the surface of the body; a procedure for obtaining these streamlines from a 


connection with a rn al coordinate systen m bh: aving the 
axis oriented normal to the surface at the actual stagnation 
point (5); indeed, with this reference system only, the 
meridional component of velocity is everywhere larger than 
the circumferential component. A perturbation analysis is 
also applicable in the three-dimensional portions of the flow 
field about a wing with sharp leading edges at moderate angles 
of attack. However, such is not the case if the leading edge 
of the wing is blunted and the streamlines exhibit a large 
curvature in the neighborhood of the stagnation line. Also 
the extension of the approximation to more general three- 
dimensional flows is not apparent. 

It is shown in the present paper that the small crossflow 
approximation is valid for an arbitrary streamline pattern in 
the presence of highly cooled surfaces and of moderate Mach 
number outer flows; hence, the laminar heat transfer over 
blunt nosed bodies in hypersonic flight can readily be deter- 
mined once the inviscid flow field is known on the surface. It 
is shown also that the results obtained from the proposed ap- 
proximation compare favorably with the exact predictions of 
2 and 3). 

The subject method requires that the streamlines of the 
inviscid flow be known. In the general case of an arbitrary 
three-dimensional body a complete theoretical analysis of the 
flow field may not be available; however, the streamlines 
may be determined if the pressure distribution on the body is 
known either from experiments or from approximate theoreti- 
‘al estimates. A. procedure for describing the streamlines and 
for evaluating the metric of the streamline coordinate system 
from pressure information is outlined in two appendixes. In 
view of the approximate invariance of the ideal gas flow field 
about blunt bodies in hypersonic flow, pressure measurements 
performed at comparatively low hypersonic Mach numbers 
can be used in the analysis of the streamline pattern for a 
wide range of free stream conditions is ise beyond those 
attainable in experimental facilities. 


Analytical Considerations* 


With reference to an orthogonal curvilinear coordinate 
system (21, 22, v3) having the surface x3 = 


* A list of definitions of symbols appears at end of paper. 
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It is nee that in the presence of highly cooled surfaces and of moderate Mach numbers at the 


outer edge of the boundary layer, the crossflow in the boundary layer is negligible even for large 
From this consideration a rapid method for estimating the heat 


0 coincident with the 
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ase 1 Op + (22 ogi [1b] 

hy Oz 2 hs hyhe hihe O22 


= 0 =a =90 H = dH, 


surface of the body, the equations governing the steady mo- 
tion in the laminar boundary layer in the absence of body 
forces are 


hs 


If the lines (22 = constant, x3 = 0) are chosen coincident with 
the streamlines of the inviscid flow, the boundary conditions 
are 


Also the momentum equations governing the inviscid flow 
become 


Ge Ore _ Op 
hy Ox 
If we denote by 
_ 1 
hyhz O22 


the component of the outer streamline curvature in the plane 
tangent to the surface of the body, the momentum equation 
in the 2-direction can be rewritten in the following form 


hy he O22 hz hy ho Ox) hs O23 


2 


When either x — 0 or [(q:?/qie2) — (p-/p)] ~ 0, Equation [5] 
is homogeneous in gq: and with the homogeneous boundary 
conditions (2) admits of the solution gq. = 0. The condition 
x — Ois satisfied by quasi-two-dimensional flows as discussed in 
the preceding section. The condition [(q:?/q.?) — (p-/p)]—~0 
is closely approximated when the wall is highly cooled and 
ne/(H. — Hw) <1; for example, this is the case of the 
laminar boundary layer on the nose portion of an arbitrary 
three-dimensional blunt body in hypersonic flight. The 
neglect of the right-hand side in Equation [5] when x ¥ 0 
can be justified on physical grounds by extending Lees’ 
observation that when p./p. > 1, the velocity and, even 
more, the enthalpy profile near the surface are affected only 
slightly by the pressure gradient, in contrast to the more 
familiar case of moderate temperature differences across the 
boundary layer (7). Under these circumstances a small 
crossflow is to be expected; hence a complete solution can be 
determined conveniently by expanding in powers of q. Only 
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the first term in this series (corresponding to the zero crossflow 
approximation q: = 0) is considered in the present paper. The 
pertinent similarity variables and results are presented; also 
the range of validity of this first approximation is investi- 

gated by comparison with the exact results of (2 and 3). 
If the crossflow is negligible, and the length elements for 
the streamline coordinate system are chosen so that 
1 


(pm) + (pq3*) = 0 [6a | 


p dx p Ox; Oz: 


where 


Introducing the new independent. variables* 


s= f, 


die 
=- 
(2s)'* Jo par: 


the stream function = (2s)'/*f(m) defined by 
oy oy 
9 
ors 
and the stagnation enthalpy ratio g defined by 
[10 


together with the assumptions that: 

1 At any station (2, pu = = constant across the 
boundary layer, and 

2 H., = constant or H,.< H,., Equations [6] can further 


_ be reduced to the following 7 


9: 
[? - =0 


Gie OS 
with the boundary conditions 
n =0 f(0) = = =O 


Equations [11] and the boundary conditions [12] are the same 
as those derived by Lees (7); their solution can be approxi- 
mated by the Blasius profile if (H../H.) < 1, q.2/(H. — Hw) 
< 1, and the effects of the vorticity in the inviscid flow are 
neglected. Under these conditions the boundary layer at a 
station distant / (actua] are length measured along the stream- 
line) from the stagnation point on the body of interest is 


4 The choice of the reference condition p;u, is discussed at the 
end of this section. 
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from the leading edge.® 


equivalent to that sc on a flat plate at a station distant x 


ma 


Evidently the solution loses physical 
significance when z < 0. Conversely, the condition zx = 0 
can be used to establish the limit of validity of the boundary 
layer approximation at the corresponding point along a stream- 
line; as suggested in (6), one may conjecture the formation of a 
wake beyond the point in question. This consideration is of 
particular interest for the analysis of the viscous flow over the 
leeward side of an axisymmetric blunt nose at large incidence. 
When the condition z > 0 is satisfied, the heat transfer at a 
general point on the surface of the body is® 


kw / 
(H. — H.)(0.47 Pr'/) Prbtrhe [13] 
Mul pw 
= 0.47 — He) hs 
with s given by 


and h, determined according to the procedure outlined in Ap- 
pendix B. Thus the heat transfer at any point on the surface 
of the body can be obtained readily when the streamlines of 
the inviscid flow are known. 

The product p,u, appearing in Equation [13] is to be selected 
suitably case by case; for example, an appropriate value for 
evaluating the heat transfer under hypersonic conditions, in- 
cluding also the effects of dissociation, is (peMe)??(Pwéw)? (8). 

If the inviscid flow is highly rotational, a twofold effect on 
the heat transfer arises (9). For one, the dimensionless stag- 
nation enthalpy profile and, in particular, its slope at the wall, 
are affected by the new boundary condition on the velocity 
profile at the outer edge; hence, the factor 0.47 in Equation 
[13] should be changed accordingly. Secondly, the entropy 
at the outer edge may differ markedly from the inviscid flow 
prediction at the surface of the body, due to entrainment of 
mass flow in the boundary layer; under these circumstances 
the relations between the local values of pressure, density and 
velocity at the outer edge cannot be established a priori but, 
instead, must be determined from an estimated variation of 
entropy with mass flow while proceeding forward step-by-step 
and evaluating the entrained mass flow at any station. This 
second effect influences the values of p,u,, gic and s to be used 
in Equation [13]. When these modifications are included, the 
subject method of analysis can be applied also in the presence 
of highly rotational outer flows. 

In the following section the accuracy of the proposed 
method of analysis is tested by comparing the pertinent ap- 
proximate results against some known solution; it is shown 
that the comparison is satisfactory. 


Comparison With Known Exact Solutions 


The boundary layer at the stagnation streamline on a 
yawed infinite cylinder (2), and the boundary layer on the 
windward side of the pitch plane for a cone at large incidence 
(3) will be considered. 

On the yawed infinite cylinder, the streamline divergence at 


5 Notice that z depends essentially on the inviscid flow pattern 
and, to some extent, on the wall temperature ratio because of 
the factor prpr. 

6 The dependence on the Prandtl number (Pr'/s) is suggested 
by the flat plate analogy established above. 
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a 


the axis is independent of the axial coordinate and, using the 


notation of (2), is given by 


dy 


Hence 


=e 
when the boundary condition h, = 1 at y = 0 is imposed. 
Correspondingly, the similarity variable 7 defined by Equa- 
tion [8] becomes 


[15] 


1 du. z 


It can be shown that this variable is the same as that used by 


2p [ hey | Jo 


Reshotko and Beckwith (whe also Pr Pwltw); 
namely 
\ ~ [16] 
Vo 
a 
dX 
w to a Ae 
O po tw Qo Po Aopo 
Le 
Indeed, at the stagnation line 
dX Mw tyo Pe dx 


Substitution of this expression into Equation [16] reproduces 
Equation [15]. Hence a direct comparison can be established 
between the slopes of the nondimensional enthalpy profiles 
at the wall as predicted by the present method (g..’ = 0.47) 
and by the calculations of (2). The results of this comparison 
are presented in Fig. 1, where the ratios between Reshotko and 
Beckwith’s values (gp’) and the value suggested here are 
plotted for several wall temperature ratios against the Mach 
number M, at the outer edge of the boundary layer. For an 
ideal gas, the Mach number M, is related to the parameter 


ty/txo of (2) by the equation 
t 
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Fig. 1 


Comparison between exact and approximate heat transfer 
rates at the stagnation line on a sweptback cylinder (y = 


1.4) 


= 
1e 
: 
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Inspection of Fig. 1 shows that predictions satisfactory for en- 
gineering purposes can be obtained by the method of this 
paper for low wall temperature ratios (t/t < 0.25) and 
moderate outer flow Mach numbers (M, S 2). These condi- 
tions are encountered in the neighborhood of the blunt leading 
edge of a sweptback wing in hypersonic flight if the sweep 
angle A is below a certain value, A» depending on the free 
stream Mach number (Ap, ~ 45 deg for M.~ = © in an ideal 
gas). 

The boundary layer on the windward side of the pitch plane 
for a cone at large incidence will be considered now. With 
the notation of (3) the streamline divergence for ig —_ 


flow is given by dn 
1 1 Owe 1 1 


where z is the radial distance from the apex, ¢ is the angular 
coordinate around the cone, ® is the sine of the cone half- 
angle, u. and w, are the radial and circumferential com- 
ponents of velocity, respectively. Hence 
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Fig.2 Comparison between exact and approximate heat transfer 
rates on the windward side of the pitch plane for a cone at large 
incidence (y = 1.4) 
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Fig. 3 Cone at large incidence—comparison between exact and 
approximate enthalpy profiles on the windward side of the pitch 
plane 
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when the boundary condition hz = 1 at x = 2 is imposed. 


Consistent with Equation [17], the expression of the similarity 


variable » defined by aia [8] is 


It can be shown that this variable is related to the variable 
used by Reshotko through the equation 


“wre 


If the approximation proposed in the present paper is ac- 
ceptable, the following relation should be satisfied u unde r hy 
personic conditions 


[20 
0.47 k+1 

where gp’ is the slope of the nondimensional enthalpy profile 
predicted by Reshotko for various cone angles, Mach num- 
bers and angles of attack. The actual values of the ratio ir 
question are plotted in Fig. 2 against the Mach number M 
at the outer edge of the boundary layer. For an ideal ga: 
this Mach number M, is related to Reshotko’s parameter t.” 
by the equation 


Ls & 


* 


~ 
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Again it can be seen that under hypersonic conditions (i.e., 
low wall temperature ratios and moderate outer edge Mach 
numbers) the approximate predictions of the present analysis 
provide satisfactory accuracy. As an additional term of com- 
parison the enthalpy profiles calculated by Reshotko for the 
cases t.,/to = 0, 1/t.* = 2.5, k = 0.60 and & = 1.20, respec- 
tively, are plotted in Fig. 3 together with the Blasius profile 
suggested here. 


Conclusions 


Analytical considerations and comparisons with known 
exact solutions indicate that in the presence of highly cooled 
surfaces and of moderate outer flow Mach numbers the cross- 
flow in the boundary layer is negligible even for large trans- 
versal pressure gradients. Therefore, the boundary layer 
characteristics and, in particular, the heat transfer at the wall 
may easily be estimated by using an orthogonal curvilinear 
coordinate system that has the streamlines of the inviscid 
flow as one family of coordinate lines. Expressions for the 
_ evaluation of the heat transfer are provided depending on the 
_ geometrical characteristics of the coordinate system. If a 
_ detailed theoretical analysis of the inviscid flow field is not 

available, the geometric characteristics of the streamlines can 

be determined from pressure measurements; the details of 
the pertinent analysis are described in the two appendixes. 
The effects to be encountered in the presence of highly rota- 
ional outer flows have been reviewed; the pertinent modifica- 


been indicated. 

The subject results indicate that the heat transfer over the 
nose section of either blunt three-dimensional bodies or air- 
foils in hypersonic flight can be predicted within engineering 
accuracy by the simple method proposed here. Inspection of 


: : tion of the final equation for evaluating the heat transfer has 
> 


the reduced data obtained from (2 and 3) also suggests that 
an even wider range of wall temperature ratios and outer flow 
Mach numbers may be covered by including a linearized term 
for the crossflow in the analysis based on streamline coordi- 
nates. 
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Appendix A 


In the general case of an arbitrary three-dimensional body, 
a complete theoretical analysis of the inviscid flow field may 
not be available; however, the streamlines at the outer edge 
of the boundary layer may be determined if the pressure dis- 
tribution on the body is known either from experiments or 
from approximate theoretical estimates (e.g., Newtonian plus 
centrifugal theory in hypersonic flow). The pertinent proce- 
dure is outlined in this Appendix. 

If the inviscid flow is rotational, the entropy at the outer 


\? dx 
i [Fzy(1 + F,?) (2) + F,%) — F.A1 + [FyyF Fy 


(1 + F,? + F,?)K? — (1 + Fz? + F,?)[Fr2(1 + Fy?) — 2F2F + Fy(1 + + (FezF yy — = 0 


edge of the boundary layer over the body in question must be 
determined first. This distribution of entropy can be ap- 
proximated with the surface value dictated by inviscid flow 
considerations’ only if the mass flow entrained in the boundary 
layer is a small fraction of the mass flow over which large 
entropy changes are distributed; if such is not the case (e.g., a 
station sufficiently far downstream on a flat plate with a blunt 
leading edge), estimates based on the entrained mass flow 
and on the variation of entropy with mass flow are necessary. 

When the entropy distribution is established, the local 
values of the velocity and of the density may be obtained from 
the pressure distribution; however, the velocity direction re- 
mains undetermined. Under these conditions each stream- 
line must be constructed by a forward step-by-step process, 
beginning at a point (say P) where the velocity direction is 
known. The unit problem pertaining to each step of the 
procedure involves two operations, namely: 

1 The determination of the plane osculating the stream- 
line at the initial point (P). 

2 The determination of the adjacent point (Q) on the 
streamline, which is defined as the intersection between the 
aforementioned osculating plane and the surface of the given 
body. 

In connection with operation 1 we observe that the com- 
ponent (Oh;/hih2d22) of the streamline curvature in the plane 
tangent to the surface of the body can be evaluated upon sub- 
stitution of the available local values of pressure, density and 
velocity into the z--momentum equation 


[A-1] 


The osculating plane is defined by the condition that the | 


curvature of its intersection with the surface satisfies Equa- 
tion [A-1]. The equation of the plane can then be derived 


7The surface value of the entropy in the inviscid flow can 

—e be determined with good accuracy. It may be constant 
., blunt or pointed bodies) or variable (e.g., a wing with sharp 
sapcanaie edges of variable thickness). 

* A locus of initial points P, where the direction of the velocity 
vector is known, may easily be determined in general; for ex- 
ample, consider the case of a blunt body where the locus is con- 
veniently chosen as a curve enclosing the stagnation point, and 
the case of either an elongated body at incidence or a wing-like 
body with a blunt leading edge, wherein initial values of the 
velocity direction can be determined on a curve adjacent to the 
stagnation line. 
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from the following classical geometrical considerations: If the 
surface is described in Monge’s form® 

= F(x, y) [A-2] 
the equation of the plane tangent at point P = (xp, yp, zp) 


is 


(x — xp)F. + (y — yp)F, — (2 — zp) = [A-3 ] 


the inclination w of the parametric curves (x = 
constant) is 


constant, y = 


F.F, 
cos ® = [A-4] 
V1 + + Fy?) 
the equation of the lines of curvature is 
Fey (1 + Fy?)] =0 [A-5] 
and the equation for the principal curvatures is 
[A-6] 


Let a and 6 be the two principal directions defined by Equa- 
tion [A-5]; let K, and Ky be the pertinent principal curva- 
tures defined by Equation [A-6]; let 6 be the angle between 
the a- and the z-directions defined by (Fig. 4) 


sin(w — | B) _ dz 
sin sin 


If @ is the angle between the streamline and the z-direction 
in the plane tangent to the surface (Fig. 4), the curvature of 
the normal section in the streamline direction is 


[A-7 ] 


Kn = Ky cos? (a — B) + Ky sin? (a — B) [A-8] 


The Equations [A-1 and A-8] require that the normal to the 
surface and the principal normal to the streamline at P make 
an angle such that (Fig. 5) 


[A-9] 


® Naturally one may select any suitable parametric representa- 
tion of the surface, and follow the procedure described below 
modifying the equations accordingly. 


Fig. 4 Detail of geometry in plane tangent to the surface 
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with m determined by the equation 


Hence, the equation of the plane osculating the streamline at P is 


y — yp k — up 
m 


sin (w — @) sin a sin @ 


sin aF, + sin(w — @)F, 


F, 
sin(w — a@) sin @ 
sin(w — @) sin? @ 
n = sin aF, + sin(w — a)F, 


As for the second operation of the unit seis ie ob- 
serve that the continuation of the streamline beyond point P 
is obtained by selecting a point Q adjacent to P on the line of 
intersection between the surface [A-2] and the osculating 
plane [A-10]. The trace of the osculating plane on the plane 
tangent at Q there determines the direction of the velocity 
vector. 

With this information available one can proceed along the 
streamline repeating the procedure described above. ) 


Appendix B 


Evaluation of the heat transfer at a general point on the 
surface of the body (see Eq. [13]) requires that the law of 
variation of the length element he(x, x2) (in the direction 22 
norma] to the streamlines) be known over the entire surface. 

One evident alternative for the determination of the func- 
tion hs(x, x2) is that of proceeding to the actual geometrical 
construction of a family of lines x2 = constant, and of measur- 
ing subsequently the are length included between adjacent 
streamlines at various stations 2,. However, the accuracy of 
such a procedure relies, to a large extent, on a suitable choice 
of the streamlines between which the intercepted arc length is 
measured. The other alternative, based on a stricter ap- 
proximation, is described below. 

Essentially he can be obtained along each streamline by in- 
tegration of the quantity Oh2/hOl; in turn, this quantity may 
be determined at selected points on the streamline by repre- 
senting the local flow as the flow in the neighborhood of the 
pitch plane on a body of revolution at incidence. The locally 
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7 Fig. 5 Schematic for geometric considerations 


K,, = Kqsin{a — B) + Ky cos(a — B) [B-2 


6 


equivalent body of revolution is defined by the followiny, 
conditions: 

1 The actual streamline coincides with the generating 
curve. 

2 The axis lies in the plane osculating the streamline (de- 
fined by Eq. [A-10]), and is oriented so that the parallel sec- 
tion through P has the same curvature as the homologous sec- 
tion of the actual body.” 

The condition 2 can easily be cast in mathematical terms 
by using the classical relations between the curvatures of gen- 
eral plane sections through a point on an analytic surface. 
Let ¢ be the line of intersection between the aforementioned 
parallel plane 7; (as yet unknown) and the plane (a, 22) tan- 
gent to the actual surface; let y be the angle between the ¢- 
and x-directions in the latter plane; let 6 be the angle betwee: 
the plane 7; and the plane zm, defined by the normal n, to th: 
surface and by the direction t; let « be the angle between th: 
plane 7, and the plane 73 determined by the normal ng to 
the surface and by the z-direction; let 


K, = Ka — B) + Ky sin*y — B) [B-1 


be the curvatures of the normal sections cut by the plane 7, 
and by the plane 73, respectively (Fig. 5). The condition 2 
can then be written as 


Kz 
cos€ = 


The Equation [B-3] can be made explicit by observing that, 
with Monge’s representation of the surface as used in the Ap- 
pendix A, the equations of planes 7, 72 and 73 are, respec- 
tively: For the plane 7 (the angle y is used as a parameter) 

— zp) + — yr) + Ci(z — zp) = 0 [B-4] 
where 


A; = sin (w — a){m sin a sin y + (m — 1)[F, sin(w — a)+ 


F, sin a][F, sin(w — y) + F,sin y]} [B-5a] 
B, = —sin a{m sin (w — sin(w — y) + [Fr sin(w — a) + 
F, sin sina — y)+ F,sin y]}  [B-5b] 
C, = [F. sin(w — a) + F, sina]: 

{sin sin y — (m — 1)sin (w — @) sin(w — y)] [B-5c] 

and m is the root of Equation [A-11]. 

For the plane me 

A(z — xp) + Bly — yp) + — zp) = 0 


10 Notice that the normal to the actual body is not the same 
as the normal to the ideal body of revolution; the latter coincides 
with the principal normal of the streamline. 
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A, = sin y + F, IF: sin(w — y) + F, sin y] [B-7a] 


B, = —sin(w — y) — F.[F. sin(w — y) + F, sin y] [B-7b] 


C, = F, sin y — Fy sin (w — y) [B-7¢ ] 


For the plane 73 
Az — zp) + By — yr) + C2 — = [B-8] 
where 
A; = cos a(F,? + 1) — cos (w — a)F:F, [B-9a 
B; = cos (w — a)(F,2 + 1) — cos aF,F, [B-9b ] 
C; = F, cosa + F, cos (w — a) [B-9e | 
The Equation [B-3] then becomes, in extenso 


+ + CC, _ Key + BiB: + 
+(A;? + C;2)'/2 Kit +(A,? + + C.2)'/2 


[B-10] 


Solution of Equation [B-10] for y identifies the position of the 
axis for the equivalent body of revolution at incidence, which, 
therefore, is completely defined. If the flow about this body 
is referred to a cylindrical coordinate system z, r, ¢ (x coin- 
cident with the axis of symmetry) with corresponding velocity 
components u, v, w, the streamline divergence in the pitch 
plane is given by 


1 ow 


From the previous geometrical considerations one obtains 


Fig. 6 Schematic of section with parallel plane 


= pressure 

Pr = Prandtl number 

1, Y Ys = velocity components in curvilinear coordinate sys- 
tem 

q3* = transformed velocity component (Eq. [7]) 

s = transformed coordinate along streamlines (Eq. [8] ) 

t = temperature 

u,v, W = velocity components in the notation of (2 and 3) 

22,23 = orthogonal curvilinear coordinates 

t, Uy 2 = space coordinates in the notation of (2 and 3) 

a, B = angles defined in Appendix A 

7, 6, € = angles defined in Appendix B 

n = transformed coordinate normal to body surface 
(Eq. [8]) 

v = angle defined by Equation [A-9] 

k = component of streamline curvature in plane tan- 


gent to surface 


lor = K..tane = K [(AiB; = A;B,)? (A.C; ont A;C;)? + (B,C; 
z 
r ol A\A;3 + + CC; 
The remaining term in Equation [B-11] can be evaluated ies 
readily in the following way. Consider the intersection be- uw = absolute viscosity Ste 
tween the actual body and the parallel plane through the defined in Appendix B 
point P of interest (Fig. 6). From the known distribution of ° = density . ‘ 
velocities, determine over the entire contour C' the components 
w = angle defined by Equation [A-4] 
of velocity v,; tangent to this contour. Plot v, vs. a cireum- 
ferential variable and determine Subscripts 
r 4 & = at the outer edge of the boundary layer 
ou ow \ 0 = at stagnation 
r = reference conditions 
w = at the wall 
Since (u? + v*)'/: is the known resultant velocity at the point ay = partial derivatives (in the Appendixes) 
P, the quantity Oh2/h» Ol can thus be determined. Repetition 
of the calculation outlined above for several points along a So 
streamline and subsequent integration along the same yield 
the law of variation of the length element h. to be used in 1 Moore, F. K., ““Three-Dimensional Boundary Layer Theory,” in “Ad- 
Equation [13]. Applied Mechanics,” vol. 1V, Academic Press, New York, 1956, pp. 
2 Reshotko, E. and Beckwith, I. E., ‘‘Compressible Laminar Boundary 
‘f ——) —— = Layer over a Yawed Infinite Cylinder with Heat Transfer and Arbitrary 
Nomenclature Prandit! Number,’’ NACA TN 3986, June 1957. 
= principal di ior ‘ a 3  Reshotko, E., ‘‘Laminar Boundary Layer with Heat Transfer ona 
a,b a directions on the surface of the body Cone at Angle of Attack in a Supersonic Stream,” NACA TN 4152, Dee. 
A;, B;, C; = coefficients defined by Equations [B-5], [B-7] and 1957. 
[B-9] 4  Reshotko, E., ‘‘Heat Transfer to a General Three-Dimensional Stag- 
Cp = specific heat at constant pressure nation Point,’ JET PROPULSION, vol. 28, no. 1, Jan. 1958, pp. 58-60. 
d = streamline divergence (Eq. [14]) 5 Vaglio-Laurin, R., ‘‘Three-Dimensional Laminar Boundary Layer with 
F 4 = function related to the stream function, y/(2s)'/? Small Cross sheng about Blunt Bodies in Hypersonic Flight,’’ GASL Tech. 
F = equation of body surface in Monge’s representation 
— ayes, e nree-Uimensiona oundary Layer, 
4 i — enthalpy ratio (Eq. [10] Report No. 1313, U.S. Naval Ord. Test Station, Inyokern, May 1951. 
= enthalpy mn P 7 Lees, L., ‘‘Laminar Heat Transfer Over Blunt-Nosed Bodies at Hyper- 
4 he, hs = length eee “on curvilinear coordinate system sonic Flight Speeds,” Jet Propuusion, vol. 26, no. 4, April 1956, pp. 259- 
= stagnation enthalpy 269. 
k = thermal conductivity; also quantity related to 8 Fay, J. A. and Riddell, F. R., “Theory of Stagnation Point Heat 
circumferential gradient of circumferential ve- noe in Dissociated Air,’’ J. Aeron. Sci., vol. 25, no. 2, Feb. 1958, pp. 
rin th 
K a locity in on — of @) rf ? the body 9 Ferri, A. and Libby, P. A., ‘‘Note on an Interaction Between the 
= curvature of a section on the surface of the body Boundary Layer and the Inviscid Flow,” J. Aeron. Sci., vol. 21, no. 2, Feb. 
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An Experimental Investigation of 
Blunt Body Stagnation Point 
Nelocity Gradient 


BOISON, 
and 


HOWARD A. CURTISS? 


Aveo Research and Advanced 
Development Division 
Wilmington, Mass. 


An experimental investigation of the stagnation point velocity gradient of five axisymmetric 
spherical-segment shapes ranging in bluntness from concave to the equivalent of hemispherical is 
reported. The tests were conducted in the California Institute of Technology Jet Propulsion Labora- 
tory 20-in. supersonic wind tunnel at Mach numbers of 2.01, 3.01 and 4.76, and at Reynolds numbers 
based on body diameter from 1 to 1.5 million. Pressure distributions were obtained to a high degree 
of accuracy by micromanometer techniques and by the use of numerous pressure orifices. A corre- 
lation of stagnation point velocity gradient with a bluntness parameter based on the sonic point 
coordinates was obtained. Modified Newtonian theory was found to be reasonably accurate for pre- 
dicting the stagnation velocity gradient of near-hemispherical shapes, but current theories for a 
truncated cylinder predicted values considerably greater than the experimental results. As a by- 
product of this investigation, two shock detachment correlations were obtained, one using an 
effective nose radius parameter to compare bodies of varying bluntness with shock density ratio, and 
the other using a bluntness parameter to obtain an almost linear variation of shock detachment dis- 


tance and shock sonic point location. 


tion of ton personic flight has placed 
emphasis on accurately predicting aerodynamic heating. 
Design of hypersonic re-entry vehicles such as a re-entering 
satellite requires reasonably accurate predictions of the stag- 
nation point heat transfer to obtain optimum configurations. 

Stagnation point heat transfer of blunt axisymmetric 
bodies of revolution at high velocity has been found from 
theoretical considerations, by Sibulkin (1)4 for example, to be 
dependent on the square root of the streamwise velocity 
gradient at the stagnation point. Further theoretical analysis 
by Lees (2) and Fay, Riddell and Kemp (3), and supported 
by shock tube experimental investigations (4), has shown 
that this stagnation point velocity gradient dependence also 
exists for hypersonic conditions with proper corrections for 
dissociation effects. 

Heretofore, theoretical (du/ds); values for nearly hemi- 
spherical bodies have been based on the pressure distribution 
predicted by the modified Newtonian approximation (5). 
For a true hemisphere, the pressure over a considerable por- 
tion of the body is reasonably well predicted by the New- 
tonian approximation. It is well-known, however, that for 
very blunt bodies, such as a truncated cylinder with the body 
axis aligned with the flow, the Newtonian prediction is not 
accurate. Several theoretical solutions of stagnation point 
velocity gradient for a truncated cylinder have been pre- 
sented, Hayes (6) and Zlotnick (7) among others. Unfor- 
tunately, none of these theories have been evaluated experi- 
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mentally and, in addition, the important region of body shapes 
falling between the truncated cylinder and the hemisphere 
has not been investigated by either theory or experiment. 

This investigation was initiated, therefore, to determine the 
effect of body bluntness on the stagnation point velocity 
gradient.© This included an investigation of the absolute 
value of (du/ds),; for both a truncated cylinder and the 
equivalent of a hemisphere cylinder to check the current 
theories. Since the bow shock wave shape and location are 
important for predictions of shock layer flow, the effect of 
body bluntness on shock detachment distance and shock sonic 
point locations were also investigated. 


* 
Experimental Technique 


Models 

A series of five, 3-in. diameter, blunt shapes consisting of 
spherical segments of varying radii with a common cylindrical 
afterbody were used for the investigation. In Fig. 1, the es- 
sential geometric details are shown including the relative 
bluntness in terms of nose radius (R,/D) and a second form 
based on the coordinates of the sonic point (x*/r*). The 
models span a range of bluntness in roughly equal increments 
from a concave body to a shape that should essentially cor- 
respond to a hemisphere. The spherical-segment-cylinder 
junctions were all sharp corners. 

The general location of the pressure orifices on the spherical 
portion is also illustrated in Fig. 1. The inner five orifices on 
the vertical and horizontal meridians were located primarily 
to obtain the precise pressure readings needed to determine 


5 The results of this experimental investigation were first re- 
ported in an Aveo RAD Aerodynamics Section Technical 
Memorandum by the authors (8). 

6 A list of definitions of symbols appears at end of paper. 
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the stagnation point velocity gradient. The other orifices 
were arranged in a spiral pattern to obtain a complete pressure 
distribution over the surface. Five pressure orifices were 
divided between the top and bottom meridians of the cylindri- 
cal afterbody as shown. 


The key to obtaining experimental values of (du/ds); is 
the extremely accurate measurement of the local pressure in 
the vicinity of the stagnation point. To achieve the required 
accuracy, the model orifices were connected to a multi- 
manometer containing silicone oil, a light liquid with a specific 
gravity of about 0.93. The reference column heights were 
measured to within +0.0005 psia by a micromanometer. 

To guarantee pressure stabilization, an extensive time 
period was allowed prior to recording the data. The manome- 
ter columns were read manually to avoid accuracy losses from 
the usual method of photographing the manometer board. 


fest Facility and Program 


This investigation was conducted in the California Institute 
of Technology Jet Propulsion Laboratory 20-in. supersonic 
wind tunnel. The stagnation temperatures were held to 110 
to 120 F for this test program at facility maximum total 
pressures ranging from 22 psia at Mach 2.01 to 64 psia at 
Mach 4.76. Corresponding free stream Reynolds numbers 
based on the 3-in. model diameter were 1.58 and 1.03 million, 
respectively. 

The test program consisted of ten runs as follows: 

Shape A at Mach 4.76. 

Shapes B, D and E at Mach 4.76 and 2.01. 

Shape C at Mach 4.76, 3.01 and 2.01. 

The models were all aligned with the flow axis. Both front- 
lighted schlieren and spark shadowgraph photographic tech- 
niques were used. The former includes a system for lighting 
the model while simultaneously exhibiting the schlieren. Two 
examples of front-lighted schlierens which also illustrate the 
model installations are shown in Figs. 2 and 3. In Fig. 4, a 
typical example is shown of the shadowgraphs used to deter- 
mine the shock detachment distances® presented later. 


Pressure Distribution 


Surface pressure distributions at Mach 2.01 and 4.76 are 
shown in Figs. 5 and 6 including data taken on the spherical- 
segment forebody and the cylindrical afterbody section. The 
abscissa s/s* represents the surface distance normalized with 
respect to the sonic length. This, because of the sharp corner, 
corresponds to the spherical-segment length except for Shape 
E at Mach 4.76. The sonic point falls ahead of the corner in 


® Shock detachment distance is defined herein as the axial dis- 
tance from the stagnation point to the front face of the bow 
shock. 


Fig. 2 Front-lighted schlieren, Shape D, 
Mach 2.01, knife edge vertical 
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Fig. 3 Front-lighted schlieren, Shape B, 
Mach 4.76, knife edge horizontal 
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Spherical-segment models and geometric relations 


this instance as indicated in Fig. 6. Unfortunately, all of the 
Mach 2.01 Shape E data above s/s* = 0.42 were not usable. 
However, modified Newtonian predicts that the sonic point 
would be aft of the corner if the spherical surface were ex- 
tended at this Mach number so that the sonic point is probably 
at the corner. 

Separation bubble angles aft of the corner on the cylindrical 
afterbody that were visible in the shadowgraphs are noted in 
Figs. 5 and 6. 

All of the forebody results are compared with the modified 
Newtonian prediction. As expected, the pressure level of the 
blunter shapes, other than the concave Shape A, falls well be- 
low the Newtonian values, more so at Mach 2.01 than at 
Mach 4.76. For Shape E, which is essentially a hemisphere 
because the sonic point is ahead of the corner, pressure results 
were about 43 per cent below the Newtonian prediction at 
the sonic point. 


Shadowgraph, Shape C at Mach 
4.76 


Fig. 4 
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Fig. 5 Surface pressure distributions at Mach 2.01 


J(du/ds); Parameter Development and 


The Bernoulli and modified Newtonian approximations 
were combined to produce a stagnation velocity gradient 
parameter in terms of stagnation conditions and body 
The incompressible Bernoulli equation 
stagnation region, assuming constant density, is 


= P + 3psu* 


P/P, 
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Fig. 6 Surface pressure distributions at Mach 4.76 
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ah 7 a velocity distribution data, Shape C at Mach 4.76 


multiplying both sides by r* we obtain 


- 


* 2 
sr 


Ps 


where J is defined as 


and r*/R,, is identical to J (du/ds),§ only for the modified New- 
tonian approximation. 
Equation [3] also can be written 


w 
This is obtained from Cp/Cp,,,, = cos? B = P — Pa/Ps — Po. 
8 The use of this very convenient Sis the was first suggested 
= Zlotnick (7). 
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near the stagnation point. Since r*/R, is a function of body 
geometry only, plots of Ju vs. s, if linear over a measurable 
distance from the stagnation point, directly measure J (du/ds),. 
The velocity from Equation [1] times J is merely 


—P./P, 


is = 


which was the equation used to reduce the raw data. A typi- 
cal plot of Ju vs. s is shown in Fig. 7. The assumption of in- 
compressibility in this region has negligible effect on the 
accuracy. The average of the four meridians was used to ob- 
tain the final J(du/ds), value for each run. Scatter between 


the four meridians was generally less than 5 per cent; the 
maximum was 11 per cent for the concave Shape A. 
Bluntness Parameter and (du/ds), Variation yt 


The parameter commonly used in the past to classify the 
relative bluntness of axisymmetric bodies is the ratio of nose 
radius of curvature to body diameter, R,/D. This appears to 
be a good parameter until one considers a truncated body, 
such as Shape B or any third or higher power body where the 
ratio goes to infinity. Thus, all bodies with a flat or high 
power curve stagnation point, regardless of how rounded they 
might be past the stagnation region, are lumped together at 
infinity. Lees (9) circumvents this problem for flat bodies 
with rounded corners by introducing a second parameter in- 
volving the corner radius. 

However, an all inclusive yet simple parameter suggests 
itself in the form of a ratio of the rectangular coordinates of 
the sonic point, 2*/r* measured from the stagnation point. 
Using this parameter, it is possible to classify all shapes re- 
gardless of stagnation point contour provided that the sonic 
point is known or can be found. The variation of J(du/ds), 
with this bluntness parameter is presented in Fig. 8. This 
figure contains the key information of the investigation. The 
modified Newtonian prediction for the full bluntness range 
at Mach 4.76 is also presented in Fig. 8. It is seen that New- 
tonian theory can be used to approximate (du/ds), for bodies — 
less blunt than 2*/r* = 0.25. 


Variation with M.. and Comparison With Theory 


The variation of J(du/ds), with Mach number is pre- 
sented in Fig. 9 along with two current theoretical predic- 
tions and experimental data by Korobkin and Gruenewald 
(10). They measured (du/ds), on a hemisphere cylinder with 
results about 8 per cent below the values of this investigation. 
The Mach 2.01 equivalent hemisphere (Shape E) results of 
this investigation are slightly lower than they would have been 
if the spherical segment had extended far enough to include 
the hemispherical sonic point. In other words, Shape E was 
effectively blunter than a hemisphere at Mach 2.01 because 
the sonic point was artificially fixed by the sharp corner. 

The close agreement of the Newtonian approximation for 
velocity gradient on a hemisphere with results of this in- 
vestigation is clearly indicated. The variation with Mach 
number is strictly a function of the rearward movement of 
the sonic point as Mach number is decreased. The theory of 


Li (11) 
1 ( 
k 


becomes the Newtonian prediction, where R,, is assumed equal 
to the nose radius. 

Two alternate theories of the stagnation point velocity 
gradient of a truncated cylinder are presented in Fig. 9. The 
Zlotnick and Newman (12) and Probstein (13) results are 
almost identical. Both are intended to predict a hypersonic 
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or in terms of the parameter used herein, J(du/ds), = 0.425. 
It is clear, however, that this theory predicts velocity gradients 
at least 30 per cent highe r than the experimental results. 
The results presented in Fig. 9 indicate that J(du/ds), is 
nearing an asymptotic value as Mach number approaches 5 
Therefore, the Mach 4.76 curve on Fig. 8 probably can be 


value (y > 1). Probstein expresses the result as 


used to obtain a close approximation of the stagnation velocity 
gradient for any shape at hypersonic speeds, neglecting real 
gas effects. 
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Approximate (du/ds), Determination 


po An evaluation of a simple approximate method for calculat- 
\ ? ing the stagnation velocity gradient from conventional pres- 

‘Sis Soe 2 sure data was afforded by this investigation. The method 
consists simply of finding the surface distance s’ where P/P, 
has decayed to 0.95, and applying the following equation 
~ | Pe: ed me The equation was obtained by solving J(du/ds), = constant 
| S| for u, introducing it into the Bernoulli equation, and solving 
ee gl the result for P/P, = 0.95. 
. ; 1 This simple method is effective as long as the density p 
| is nearly constant from the stagnation point to s’. In Fig. 10 
| density is shown to be indeed practically constant for al 


shapes less blunt than about 2*/r* = 0.14. The value o! 
J(u/s) was calculated directly for each point and is seen tc 
converge on a value (after the first few inherently inaccuratk 
points) close to the results of Fig. 8. The equivalent hemi- 
sphere results are in agreement with Korobkin and Gruenewald 
(10) data which showed that density was constant for more 
than 20 deg of arc from the stagnation point (8 = 20 deg is 
indicated in Fig. 10). oe 


Shock Detachment Distance 


\ HEMISPHERE | a |0143 Correlation With Density Ratio 

| The stagnation velocity gradient parameter J(du/ds), was 
“ ninitifincs a shown earlier (Eq. [3]) to be equivalent to r*/R,, (sonic point 
ra sine ei, radius/nose radius) based on the modified Newtonian ap- 
g 02 a proximation. However, experimental results showed an in- 
a  ereasing deviation of J(du/ds), from modified Newtonian 
8 at é theory as body bluntness was increased. This deviation sug- 
“ () = that with increasing bluntness the body has an effec- 
aan tively smaller nose radius than that assumed by Newtonian 
Loy | _ theory. This effective nose radius R.;; would then be defined 
; member is the experimental value. Or, R,,;,; can be considered 
Fig. 11 Shock detachment distance parameter vs. shock density to be the radius of a sphere that has the same stagnation point 


ratio velocity gradient as the body. 

This R, yy Was used as a nondimensionalizing parameter for 
comparing shock detachment distances, obtained from this 
experiment, with density ratio (Fig. 11). This figure shows 
_ that shapes bracketing the full range of body bluntness can be 

correlated on one curve approaching a common value as k 
increases. Data from Lees (9), which extend from a trun- 
cated cylinder to a hemisphere cylinder, are included in the 
figure. 

Theoretical predictions of Serbin (14) for a hemisphere and 
a truncated cylinder are shown to approac th nearly the same 
value as & increases. Serbin’s equations in terms of R,,, for 
the hemisphere and the truncated cylinder are 


a 
© | 387 
o | 492 

© | 5.22 LEES (REF 8) 


Ress} — 1)" \dsJ, 
5 ost \ds/, 
w 
= 
3 The agreement of the latter equation with the results is very 
# 02 good even at the lower density ratios. The approach of Zlot- 
5 nick and Newman (12) is also presented on Fig. 11; their 
3 equation is simply 
[LINE THRU SERBIN 
| Ress cyl ds 


but is farther from the experimental results than the Serbin 
Fig. 12 Shock detachment distance vs. body bluntness theory. 
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Fig. 13 Shock sonic point location vs. body bluntness 


Correlation With Bluntness eatin 

Shock detachment distances nondimensionalized with body 
diameter are presented in Fig. 12 as a function of the blunt- 
ness parameter x*/r*. The variation is surprisingly linear, 
even including the concave Shape A. Note that Lees’ results 
for k = 5.22 are in excellent agreement with the k = 4.92 re- 
sults of this investigation. The linear variation is valid only 
for bodies with sharp corner sonic points, since increasing the 
hemispherical surface length beyond the sonic point lowers 
6/D still farther. The theory of Serbin for the truncated 
cylinder and hemisphere is shown connected by a dashed line 
at each of three values of k. 

An almost linear variation of the shock sonic point coordi- 
nates (x*,,/r*,,) with bluntness was found to exist, as shown 
in Fig. 13. The 2*,, shock coordinate was measured from 
the shock intersection with the axis of symmetry. This is an 
interesting correlation considering the current division of 
blunt bodies into two categories,’ i.e., round bodies where the 
shock curvature is similar to the body curvature, and blunt 
bodies where the shock curvature is largely independent of 
body shape. This correlation (Fig. 13) shows that these two 


categories can in effect be bridged. 
4 
Concluding Remarks 


A summation of the important conclusions derived from 
this experimental investigation” is as follows: 

1 Stagnation point velocity gradient can be correlated 
with a bluntness parameter based on the body sonic point 
coordinates for a range of shapes from concave to an equiva- 
lent hemisphere. 

® Lees (9), for example. 

1° Tt should be mentioned that this investigation was conducted 
using only spherically capped cylinders, and thus the results do 
not necessarily apply to all body shapes. However, the authors 
correlated stagnation point velocity gradients of a considerable 
number of random body shapes in (8) using the approximate 
method mentioned herein, yet found no distinction between 
smooth and sharp cornered shapes. 7 
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2 Modified Newtonian theory has been found to approxi- 
mate the stagnation velocity gradient for rounded bodies. 

3 Current truncated cylinder theories overestimate the 
stagnation point velocity by at least 30 per cent. 

4 Shock detachment distances correlated on the basis of 
an effective nose radius for bodies of varying bluntness tend 
toward a common asymptotic value at high density ratios. 

5 Shock detachment distance and shock sonic point are 
almost linearly related to the «*/r* bluntness parameter. 


Nomenclature 
D = body diameter eee 
f = stagnation point velocity gradient parameter, 
| 2 
ic 

k = shock density ratio, po/po 
M = Mach number 
= static pressure 
r = radial distance from body axis 
Reps = effective body nose radius, r*/J(du/ds)s 
R, = body nose radius 
Rsn = shock radius at stagnation point 
8 = surface distance from stagnation point 
s’ = surface distance where P/P; = 0.95 
u = velocity component in s direction 
(du/ds); = stagnation point velocity gradient 


x = distance along body axis 
B = ray angle, normal to 6 

¥ = ratio of specific heats 

6 = shock detachment distance along body axis 
6 = angle between surface and body axis 

p = mass density 

¢ = leading edge angle of separation bubble 


Subscripts and Superscripts 


2 = behind normal shock 

o = free stream 

s = stagnation point on body 
* = sonic point 
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Technical Notes 


Wall Temperature Due to 


Mass Transfer Cooling With a 
Combustible Gas 


GEORGE W. SUTTON! 


NE POSSIBLE method of heat protection and cooling 
for the hypersonic flight of vehicles in an atmosphere 
is the use of area mass transfer cooling, in which a cool gas 
from a reservoir inside the vehicle is forced through a porous 
external skin, as shown in Fig. 1. This has two effects: First, 
the injected gas thickens the boundary layer and thereby 
reduces the rate of energy transfer to the porous skin; second, 
the gas absorbs energy at a rate equal to the product of the 
mass transfer rate and enthalpy change between its initial 
temperature and its temperature as it passes into the external 
flow. 

Under certain circumstances, it may be desirable to use a 
combustible gas for the mass transfer coolant. For example, 
on a weight basis hydrogen is potentially more effective than 
gases with higher molecular weights. Or one may want to 
increase the electrical conductivity of the surrounding gas by 
taking advantage of the ions which exist in a flame. 

For rapid combustion rates where the flame front is very 
thin, several analyses of the boundary layer profile have been 
recently presented (1 to 6).2. A similar class of problems has 
also been treated by Spalding (7). 

Although certain differences, which are beyond this dis- 
cussion, exist between them, the results of these analyses 
indicate that for low mass transfer rates, all of the enthalpy 


Original version of this paper was presented at ARS Semi-An- 
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Schematic diagram for area mass transfer cooling at an 
axially symmetric stagnation point 


Fig. 1 


of combustion is returned to the porous skin. As the mass 
transfer rate is increased, the flame front is blown away from 
the surface, and further increases in the mass transfer rate de- 
crease the rate of energy transfer. 

In this note, the results of the analysis of (5) for the injec- 


— tion of a combustible gas at an axially symmetric stagnation 


point into a hypersonic laminar boundary layer are presented. 

In that analysis it was assumed that the Prandtl number 
of the gas was unity, that the specific heats and molecular 
weights of each component gas were identical, that the effect 
of the pressure gradient at the stagnation point could be 
neglected, and that the product of density and viscosity in 
the boundary layer was a constant. For this model gas, the 
dependence of the adiabatic wall temperature on the mass 
transfer rate and enthalpy of combustion is shown in Fig. 2. 

As indicated by the dashed line, as the mass transfer rate 
is increased, the adiabatic wall temperature will increase 
above the temperature at the edge of the boundary layer 
The adiabatic wall temperature will be a maximum at the 
intersection of the proper dashed and_ solid curves. 
Further increases in the mass transfer rate will decrease 
the adiabatic wall temperature, as indicated by the solid 
curve of Fig. 2. 

Now for hypersonic flight, the temperature at the edge of 
the boundary layer is of the order of 10,000 F, but the porous 
skin must be kept at some temperature less than 3000 F, de- 
pending on the material, in order that the skin retain its 
structural integrity. Hence we are interested in determining 
the required mass transfer rate for small values of the wall 


5 


| 
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Fig. 2 Effect of heat of combustion and mass transfer rate on 
the adiabatic wall temperature 


Epitor’s Note: The Technical Notes and Technical Comments sections of ARS JourNAL are open to short manuscripts describing 
new developments or offering comments on papers previously published. Such manuscripts are published without editorial review, 
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temperature ratio (7, — 7;)/(T. — T;). This may be ob- 
tained from Fig. 2 by drawing a horizontal line at the desired 
wall temperature ratio; the intersection of that line with the 
proper solid curve gives the required mass transfer rate. It 
is seen that the effect of combustion is to increase the re- 
quired mass transfer rate, but the percentage increase de- 
creases as the required wall temperature decreases. Thus, it 
may be possible to use a combustible gas for mass transfer 
cooling in hypersonic flight with only a small penalty of 
additional coolant weight. 

The analyses of (1 to 6) were based on rapid reaction 
kinetics, but the criteria for “rapidity” were not defined. In 
order to obtain some quantitative information, a combustion 
rate law was assumed that is equal to the product of local 
concentration of combustible gas and an Arrhenius function, 
and divided by a characteristic time 7, thatis 

wy = pre B/RT [1] 

The diffusion equation for the combustible gas was then 
integrated for a particular value of the mass transfer rate 
given by 


= 0.0707 [2(due/dx) |'/? [2] 


The results indicated that for 2/RT. = 0.1, the reaction 
rate was rapid, and the flame front was thin when 2(du./dx),r 
< 10-2, and hence the above analysis applies. However, if 
2(du./dx),r > 10, the reaction was sufficiently slow so that 
the effect of combustion was negligible. For values of 
2(du./dx).r between 10-2 and 10, the flame front has a finite 
thickness, and the assumption of a thin flame front is not 
valid. With E/RT, increased to a value of 10, 2(du./dx).r 
must be decreased to 10~* before it can be assumed that the 
flame front is thin. The full details of these calculations will 
be published at a later date. 


Nomenclature 

Coe = mass concentration of oxygen at edge of boundary 
layer 

Cs = specific heat at constant pressure 

(du-/dzx); = stagnation point velocity gradient 

E = energy of activation for combustion 

m = rate of injection of coolant gas 

Qr = enthalpy of combustion per unit mass of combus- 
tible gas 

R = gas constant 

ro/f = stoichiometric weight ratio of oxygen to combustible 
gas 

= temperature 

we = rate of combustion of combustible gas 

p = density of gaseous mixture = 

M = viscosity of gaseous mixture 7 

= characteristic combustion time 

Subscripts 

f = injected combustible gas 

e = edge of boundary layer 

a = coolant gas reservoir 

w 


= surface of porous skin 
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Effect of Acceleration on the Longitudinal 
we Dynamic Stability of a Missile 


E. V. LAITONE! 
University of California, Berkeley, Calif. 


The effect of high rates of acceleration or deceleration is 
introduced into the linearized equations of motion pre- 
dicting the behavior of the short period longitudinal os- 
cillations of a nonrolling missile having a longitudinal 
plane of symmetry and trimmed to fly at a nearly zero lift 


trajectory. 


iw LINEARIZED equations of motion for the natural 
short period longitudinal oscillations of a hypersonic air- 
craft or missile during large varying rates of acceleration are de- 
rived in (1),2see also (2). The derivation was based upon smal] 
perturbations of Euler’s dynamical equations for axes fixed 
at the body’s center of gravity so as to always coincide with 
the principal axes of inertia as shown in Fig. 1. The resulting 
equations to be solved were 


m(t% + wq) = —(mgeos yO 
= ° S oCn 
— ug) = — w — (mg sin (2] 


2 oa 
puSL 
Iq = (2 - ) w+ o( qt 

oc, 

SL? 

2 af — 
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mit) = (thrust — 3putSCco, — mg sin Yo) [4] 
which are the same as the expressions derived in (3 or 4) ex- 
cept that the large rates of acceleration add the term % to 
Equation [1], and all the terms related to a change in the 
drag coefficient or any variation of aerodynamic coefficients 
with change in forward speed, have been omitted because of 
the restriction to hypersonic speeds. 
Then if the missile is initially trimmed to fly at nearly zero 
lift, the term wq may be neglected in [1], so that 7 is approxi- 
mately defined by Equation [4], and is independent of the 
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w and q perturbations provided the gravity perturbation term 


—(mg cos ¥o)9 = —mg cos Yo 


q(7)dr [5] 


remains negligible in comparison to mtg or any of the terms 
in [4]. 
_ The simultaneous solution of [2 and 3] then yields 


u(t) + + = 


provided that we neglect the gravity perturbation term 
t 
—(mg sin = —mg sin Yo q(r)dr [7] 


Since u is varying rapidly, and its variation may not be 
monotonic, it is generally difficult to develop a stability cri- 
terion from Equation [6]. A more suitable form is obtained 
by transforming to the nondimensional variable ¢ representing 
body lengths traveled along the unperturbed trajectory or 
average, arbitrarily curved, flight path where a: 


so that Equations [6] become 


w"(E) + bew'(E) + cow(&) = 0 


d(uCn 
w(t) = p(E)SL 
2m 


It is important to note how the change in variable has com- 
pletely eliminated the apparently direct effect of acceleration 
on the b term, usually associated with the “damping coeffi- 
cient,” by eliminating the variation produced by u itself. 

Oswald (5) has similarly analyzed Euler’s equations of 
motion for the case of constant acceleration, but includes the 
effects of thrust misalignment and control surface deflections, 
by using the angle of attack (a) as the variable so that 


a = w/u = (w/u — wit/u2) w= (ua + at) [10] 
Introducing a into Equations [6 and 9] transforms them into 
a(t) + + cna(t) = 0 


bu(d) = O(C mg Cme )] = 


u u d(uCwa) 


u i\? “y 
u 


fui- OC ma — HOC + 


wee 


[-0Cng ~ HOC ng + 5 — 
p(t)SL 


case as tig > 0 in [4], to 


u’ “ 
u 


which are in exact agreement with Oswald (5) for the case of 
constant acceleration, no thrust misalignment, and neglected 
gravity forces. Now, however, the effect of acceleration 
has been completely altered. A rapid deceleration can ap- 
parently reduce or even eliminate 6, the equivalent damping 
coefficient for the angle of attack oscillations, whereas de- 
celeration either seemed to aid the damping of w(t), or had 
no effect on the damping of the motion of the center of 
gravity of the body normal to its flight path, w’(é). The 
reason for this apparent anomaly is evident from Equation 
[10], where it is seen that deceleration can increase & even 
though it is decreasing w. 

This case of rapid deceleration would be very important 
during the re-entry of a ballistic missile as investigated by 
Friedrich and Dore (6) and Allen (7) in whose work the de- 
celeration is taken to be 


ws ws 
Coo g sin Yo 3 4Co [13 
m 


Substituting Equation [13] into [11 and 12] and noting that 
Cn. = (Cig + Cr) Co, = Coy {14} 


we obtain expressions which have the same b or ‘“‘dampin; 
terms,” but have the c or “spring constant” terms differen‘ 
from those derived by Friedrich and Dore (6) or Allen (7) 
For example, Equations [12] reduce to 


= wlCrg — Cry — O(Cmg + 


= — poCn + (Cra + Co) 


so the 6 term is identical to Allen’s (7), while the c expressions 
differ by the term 


[16] 


Fortunately, Allen (7) has shown that in the typical hyper- 
sonic case all the u? and yw’ terms may be neglected so Equa- 
tion [15] and Allen’s result both reduce to 


a” + — Cry — + Cre) la’ + 
po(—Cnaa [17] 


In (2) an explanation is given for the difference between 
the results obtained from the linearization of Equations [2 
and 3], and those obtained by Friedrich and Dore (6) and 
Allen (7) who linearized the flight trajectory equations which 
keep the z-axis always tangent to the flight <cauieniliaihes 1) 
[18] 


y= — 
mV = — spV*SCp — mg sin 


mVY = 4pV2SC_, — mg cos [19] 
* 


It is shown in (2) that in the linearization procedure the 
gravity term in [19] has to be neglected in order to obtain a 
linear equation for a. This results in eliminating a perturba- 
tion term that corresponds, in the limiting case as Vo — 0, to 
(mg sin > = —3pV02SCp,(60 — ba) [21] 
On the other hand, in the linearization of Equation [2], the 
gravity term given by Equation [7] has been neglected and 
this perturbation term corresponds, again in the limiting 


Consequently, it is not surprising that the linearized pertur- 
bations based upon different coordinate axes (one set being 
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fixed in the body, while the other set always moves so as to 
be tangent to the flight path), should result in different terms 
being omitted. It is interesting to note that the elimination 
of the uw? terms brings both linearized equations into agree- 
ment, and yu is generally small, being approximately 107° 
at sea level and 10~4 at a 10-mile altitude. Consequently, 
the limiting case given by Equation [17] should be valid for 
hypersonic speeds at high altitudes, and this is essentially the 
expression studied by Allen (7) for the straight line trajectory 
in an exponentia] atmosphere. A useful upper bound for 
Equation (17) for any arbitrarily curved trajectory is given 


by Laitone (8). 


(u, w) = velocity of ve- 
locity vector V, along (z, z) axes, respectively 
(see Fig. 1) 

= d@/dt = angular velocity component about y axis 


Nomenclature 


q 

6 = angle of rotation about y axis = /qdt 

a = angle of attack (see Fig. 1) 

> = flight path angle relative to horizontal 

I = mass moment of inertia about y axis = 
S (x? + 2*)dm 

m = W/g = mass of body 

p = mass density of atmosphere 

m = pSL/2m = nondimensional ratio proportional to 
mass of displaced air to actual mass (m) of 
body 

8, = reference area and reference length, respec- 
tively 

0 = _mL2/I = nondimensional reciprocal of square of 

radius of gyration 


6 oe = standard NACA coefficients for lift and drag 
forces perpendicular and parallel to flight 
path, nondimensionalized by dividing by 


(1/2 
Cm = standard NACA moment coefficient about 
- center of gravity, nondimensionalized by 


dividing by (1/2 pV2SL) 
Cy, Co = standard NACA coefficients for forces normal 
and parallel to principal body axis (x), non- 


= dimensionalized by dividing by (1/2 pV?S) 
(°) 7 = d(_)/dt, derivative with respect to time (¢) 
= d(_ )/dt, derivative with respect to nondi- 
mensional distance traveled along flight 
‘1 path (see Eq. [8]) 
Oa 
Cnq = Cme = /uy where a > 0,q¢—>0, 0 
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Rapid Method for Computing High 
Altitude Gravity 


e Technology Laboratories, Inc., Los Angeles, ( 


HIS note describes a method for computing powered 

flight trajectories of a missile at high altitude in near 
horizontal flight. The method gives closed form solutions 
to the velocity, flight path angle and altitude changes during © 
powered flight under the following assumptions: No drag, — 
near horizontal flight (@ + 30 deg from horizontal), spherical — 
Earth, general initial conditions, and gravity turn (a = 
zero de 

The method is best used to compute the performance of 

final power stages where booster data up through the at- 
mosphere is already available. 


Equations of Motion 


The equations of motion can be written along and normal 
to the flight path as follows 


Fg 
F 
= — =! sin adt + (v.2 — v?) dt 
W r 
Replacing dt with —dW/W, where W = F/J,,, and setting 


a = Ogives 


wi 
dem — + cos 
WF 
vdB = — --sinB [1 - (: ) 
r 


Since 6 is close to 90 deg, sin B & 1; cos B B (w/2) — B. 
Also, it is convenient to use weight ratios rather than absolute 
values 


uf 
Wo 


where Wo is the weight at start of burning. The equations 
dw 


of motion now become 
R 
ly = —Ig — dw 


It is noted that r and v,, which are the radius from the center 
of the Earth to the missile, and the orbital speed at that 
altitude, respectively, are averaged over the altitude incre- 
ment traversed, and are treated as constants. 

The above equations are nonlinear, and cannot be solved 
directly. The technique used here is to obtain approximate 
solutions to each equation in turn by first assuming solutions 
to the alternate equation. For example, the second term on 
the right-hand side of the dv equation is the gravity loss term, 
and is usually small in comparison to the first term. Hence, 
a good approximation to velocity is 


v =v 


This expression is then inserted into the d@ equation, which 
now becomes a function of w only and certain constants. 
As such, the expression can be solved in closed form for ~. 
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Fig. 1 Coordinate System 


Having achieved a solution for 6 in terms of w, 8 can now be 
eliminated from the dv equation, making it a linear function 
of w. A direct solution for v can now be written, which gives 
a reasonably close approximation to the gravity loss term. 

Since v and 6 are known, the altitude can be obtained by 


dh = cos Badt 
h=h + Sf v cos Bat 
The results of this analysis are given below? 


AB = 8, — = Ail —w)+A:Inz—A;wlnw 


Av = Avideat 
1 [4| 
2 ~ 2 In derivi deriving [3], In w was replaced by 2[(w — 1)/(w + 1)] in 
the 1/v term of Equation [2]. 


T T | 
=4(1-w)? 16 | 
97° (In? w - inw +4)] 
= 


Fig. 2 Graphs for use in nn gravity turns 
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Fig. 3 Comparison of method with IBM 704 machine results 


F/Wo _ % Adgrav 


T 


2 
» 
Ay = (72) 
F/Wo 
b(v + Ig) 
Ay = Ao 4) 
Ay = (: *) 
b 
As = — Ao 
Ue 


The functions are given in Figs. 2 (a-d). 


Example 


An example has been worked out to show the authenticity 
of the method. <A typical booster burnout condition was 
selected from a machine computation. Final stage per- 
formance was then computed by the above method, and com- 
pared to that obtained on the machine. The results are 
shown in Fig. 3. Av, A@ and Ah are within 0, 4.4 and 3.3 
per cent, respectively, of the exact machine results. 

It should be noted that if charts of the type illustrated in — 
Fig. 2 are used to aid in the computations, accuracy out to 
the fourth significant figure must be available for the func- 


tions ge and ge. 


b 

PF 

g = gravitational constant, ft/sec? 
h = altitude, ft 
I 

r 
R 


Nomenclature 


= constant = vw + 21g 
= constant = uw — 21g 


thrust, lb 


a 


= specific impulse, sec 
= radius from center of Earth 
= radius of Earth 


t = time, sec 
v = velocity, ft/sec ae) 
w = weight ratio = W/W, — 


W = weight, lb 


z = (a+ b)/(a + bw) 
a = angle of attack, deg 
8 = flight path angle, measured from local vertical, deg 

Subscripts 7 

=! 

0 = conditions at start of burning 
c = circular or satellite conditions (ve = +/gr) : 


See Fig. 1 for additional definitions. 
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| Energy Transfer at a 
Chemically Reacting or Slip Interface’ 


SINCLAIRE M. SCALA? and GEORGE W. SUTTON: 


General Electric Co., Philadelphia, Pa. 


HERE is considerable interest in the determination of the 

heat transfer to a chemically reacting surface. Hence, a 
number of recent investigations (1,2,3)4 have presented ex- 
pressions for the heat transfer rate at a wall with mass transfer 
and chemical reactions. Since the equations given in the 
above references are not complete, an analysis is given here 
which considers the derivation of a more exact expression for 
the energy transfer, starting from first principles. 


Analysis 

The conservation of energy for a multicomponent chemi- 
cally reacting gas may be written equivalently either in dif- 
ferential form (4) 


dE 
pa = + [1] 
— => i 


or in integral form - 
pEdV = — ££ (Q-n)dA — (pEv-n)dA + 
FF (w-n)-vdA + FFF (2) 


The right-hand side of Equation [2] represents the absolute 
energy flux through a stationary control surface, and de- 
pends on the contributions of the relative energy flux vector 
Q, convection, and the work of the pressure, the shear stresses, 


and the external forces. 

In our notation, the energy flux vector which includes the 
effects of conduction, diffusion and radiation may be written 
(4) 

Q = -KVT + 
— (Vi - +Qe 131 


Ne MiDiz o> 


while the pressure tensor yields 


= (—pl + = —pv + [4] 


7 ; If external forces are neglected, the right-hand side of Equa- 
tion [2] becomes 


—$$(Q + pEv — T-v)-n dA [5] 


—> 


Substituting [3 and 4] into [5], and neglecting the Dufour 
term, the absolute energy flux into the control volume is given 
by 


—£$(—KVT + + pEv + pv — 7ij-v)-ndA [6] 
or equivalently by 
-$$| —KVT + + pv (304) dA [7] 


or 


-$$| -KVT + + peo — [8] 
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In general, the chemically reacting surface is in motion, and 
there may be both relative tangential and vertical compo- 
nents of fluid velocity, which result from the convection of a 
liquid phase or surface slip and the interphase mass transfer, 
respectively; and thus, the term involving the shear stress 
does not necessarily vanish. 

When a Prandtl-type order of magnitude analysis is applied 
to Equation [8], one obtains as the resultant absolute heat 
transfer to a surface 


ou 
(—Quw)gas = E Vih: — puH + pu [9] 
rey) 


where the minus sign indicates heat transfer toward the sur- 
face. The first and last terms of Equation [9] have only re- 
cently been considered (5). 

This same equation applies at either side of the surface. 
Thus, looking from a condensed phase side which is devoid of 
diffusion, the net heat transfer into the gas is 


—Quiquia 
or 


ou 
— + puh — pu ‘) [10] 
solid 


oy oy 

At a stagnation point, or if the condensed phase is a solid, 
the last term vanishes. Thus, the energy balance at a wall is 
given by 


[x2 
oy 


where in most practical applications, h, H.. 
siders a vaporizing surface for which 


mut | - mut | [11] 
gw oy cw 


If one con- 


h oo = h cow + L we 12 
(hk)o (hx) [12] 
then one obtains Le 


K + mL [13] 


Note that the term m(h — h,),.. can be omitted only for the 
special case where the enthalpy of the injected gas, species K, 
is identically equal to the enthalpy of the gaseous phase. 

For mass transfer cooling utilizing a gas, if it is assumed that 
the heat conduction and flow within the wall is one-dimen- 
sional, an energy balance requires 


o7 
— ivih =|K 
[ oy | 

such that, usually 
lim (x ) =0 

Hence, if one selects a control volume including a point 
_ within the fluid reservoir, one obtains 


, or 
gw 


For mass injection through a semipermeable wall, such that 
none of the freestream fluid diffuses into the reservoir or reacts 
at the wall, there follows 


oT 
(x 
oy 


Thus, for a single component fluid, e.g., air-air injection, 
one obtains the classical result 


= [14] 


res,w 


hives, [15] 


= Mdhgw — hres,—o) [16] 


oT 
(x2 ) = Ng. w hres, [17] 
Oy J o.w 
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where h,.,, is the enthalpy of the gaseous mixture evaluated at 
the surface temperature. 

It is stressed that for mass transfer into a multicomponent 
stream, Equation [17] is not generally correct, and one should 
utilize either Equation [15] or [16]. For the special case that 
all chemical species have precisely the same enthalpy evalu- 
ated at the surface temperature 7’,,, there follows 


and then E [16] reduces identically (17). 
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Nomenclature 

A = area } gi?, 

Ci = mass fraction of the ith species y 


Di; = binary diffusion coefficient 

D,* = thermal diffusion coefficient 

ei = internal energy of the ith species, including chemica 

E = stagnation internal energy, including chemical = 
e + (4)e? 

F; = external force, acting on unit mass of the 7th species 

hi = enthalpy of the ith species including chemical = 
S cpidT + Ah 

Ah;;,° = heat of formation of the ith species . 

h = enthalpy of the gaseous mixture, }>cshi 

H = stagnation enthalpy, h + (4)v? 

I = unit tensor 

K = thermal conductivity 

L = latent heat of vaporization : 

M = molecular weight 


m = pv = mass transfer rate 
on = unit outward normal 
= 
n = moles per unit volume 
Q = energy flux vector 
R = universal gas constant 
t = time 
= temperature 
u = tangential component of velocity 
v = normal component of velocity 
n = macroscopic stream velocity 
Vi = absolute velocity of the ith species 
Vi = diffusion velocity of the ith species 
V = volume 
x,y = coordinates, parallel and normal to surface 
= pressure tensor 
Tij = viscous stress tensor 
p = density 
= ordinary viscosity coefficient 
= surface integral 


LFF = volume integral 


Subscripts 

c = condensed phase 
g = gas phase - 

t = ith species 

k = injected species 

R = radiation 

res = fluid reservoir 

t = total 

w = surface 
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Multiple Flameholder Arrays: 


Flame Interactions’ 


WRIGHT? 


Jet Propulsion Laboratory, California Institute of 
Technology, Pasadena, Calif. 


RRAYS of bluff body flameholders are used to increase 
flamespreading rates over the rate obtainable with a 
single flameholder. Multiple arrays do, indeed, lead to rela- 
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BLOCKAGE RATIO 
Fig. 1 Peak blowoff speeds for }-in. circular cylinder flame- 
holders 


a Symmetric flame configuration 


tively short combustion chambers. In fact, with several 
flameholders arranged in a single plane normal to the flow, 
‘ach holder operates nearly as though it were a single flame- 
holder in a small duct. The flame quickly fills the small 
duct, and the combustion chamber is short. However, since 
the blockage of the flameholder in this small duct is high, and 
since blockage adversely affects flame blowoff speeds, the 
blowoff speeds for multiple arrays are relatively low. Blow- 
off speeds for multiple arrays may be predicted from results 
showing the influence of blockage on a single flameholder 
(Fig. 1). In this figure the five flameholder point (blockage 
ratio = 0.31) is anomalously low, partly because of extraneous 
pressure and velocity fluctuations in the duct, but also partly 
because of an interesting flame interaction that is exemplified 
in Fig. 2. Two of the flames in this figure are pinched down 
and blow off prematurely. 

Pinching of alternate flames is observed with many multi- 
ple flameholder arrays, particularly for closely spaced flame- 
holders and for conditions approaching blowoff. This flame 
configuration is apparently stable, but interestingly enough, 
another stable configuration is often possible for the same 
average flow conditions: All flames may be equal. A com- 
parison of the two configurations is perhaps more easily made 
by looking at Fig. 3, which shows the flames held on two flat 


Fig. 2 Schlieren photograph of flames held on five }-in flat plate 

flameholders in 3 X 6-in. duct. Flat plates oriented normal to the 

stream and spanning the 3-in. dimension of the duct. Photo- 
graph shows the asymmetric flame configuration 


b Asymmetric flame configuration 


Fig. 3. Schlieren photographs of flames held on two flat plate flameholders in 3 < 6-in. duct 
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Fig. 4 Time exposure photograph of flames held on five flat 
plate flameholders in staggered array in 3 X 6-in. duct 


plate flameholders set normal to the flow. The flames may 
be symmetric as in Fig. 3a, or they may be asymmetric as in 
Fig. 3b. Either upper or lower flame may be pinched off in 
the asymmetric configuration. Any one of the configurations 
may exist for a relatively long time and may shift to another 
configuration without apparent change in average flow condi- 
tions, although shifts from one asymmetric configuration to 
the other are rare, and possibly involve passage through the 
symmetric state as an intermediate step. The flame inter- 
actions are somewhat reminiscent of the jet interactions ob- 
served by Corrsin in isothermal flow.* In his work, adjacent 
jets of a multijet two-dimensional array were found to interact 
if the blocked area between jets was sufficiently large com- 
pared with the jet open area. Similarly, the flames have 
greatest tendency to interact if the open area between flame- 
holders is small compared with the area blocked by the flame- 
holders. 

Pinching of one or more flames also occurs with staggered 
arrays of flameholders (Fig. 4), but in these cases the pinched 
configuration is usually the only one observed. Staggered 
arrays have, in general, higher blowoff speeds than coplanar 
arrays of the same number of flameholders; blowoff depends 
critically on the exact configuration. For example, with 
flameholders in two planes normal to the flow, blowoff speed 
depends strongly on distance between the planes. At small 
separations, the forward flames are pinched as the flow con- 
verges to pass between the aft flameholders, and the forward 
flames blow off prematurely. On the other hand, with large 
spacings between the planes, the flow behind the forward 
flameholders accelerates as the flames spread, and the down- 
stream flames are pinched off. The critical separation deter- 
mining whether forward or aft flames are pinched is about one 
recirculation-zone length, but depends on the relative block- 
ages in the two planes. 


“Investigation of the Behavior of Parallel Two- 
” NACA ACR 4H24, Nov. 1944. 


3 Corrsin, 
Dimensional Jets, 


Effect of Unequal Wall Roughness _ 


HE CASE of a fluid flowing in a duct with walls of unequal 
roughness appears to have received little attention from 
investigators, although it is of importance in practical systems, 
e.g., condensers, rocket film cooling. In a recent paper, 
McManus (1)? reported that in a horizontal two-phase air- 
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on Flow Between Flat Plates 
ae H. N. MCMANUS Jr.! 
Cornell University, Ithaca, N. Y. 


water system undergoing annular flow, the wave height in the 
lower portion of the tube was severalfold that in the upper 
portion. In a previous investigation of velocity distribution 
in an annular system, Krasiakova (2) reported that the 
maximum gas core velocity was significantly displaced from 
the tube center line. Although it is evident from these in- 
vestigations that the unequal roughness felt by the gas core is 
responsible for the skewing of the velocity profile, there ap- 
pears to be at present no method of estimating the magnitude 
of the displacement of the maximum velocity point from the 
tube axis brought about by different degrees of roughness nor 
for the effect upon pressure drop or weight rate of flow. 

That unequal wall roughness can result in substantial dis- 
tortion of the velocity profile with a consequent alteration in 
pressure drop and weight rate of flow can be shown by examin- 
ing the flow of a fluid between flat parallel plates of unequal 
roughness. Although the previously cited reports dealt with 
two-dimensional systems, the one-dimensional case is more 
easily treated, and, importantly, no assumption need be made 
regarding the variation of roughness around the tube pe- 
riphery. In the following treatment, it is assumed that the two 
portions of the full velocity profile are described by the uni- 
versal law and that the Reynolds number is large enough so 
that the wall friction is in the completely rough regime. 

If we locate our coordinate origin as shown in Fig. 1, each 
portion of the turbulent core velocity profile will be described 


by 


= 2.50%» In [30(d2 + y)/k3.] 


where vy is the friction velocity and k, the conventional sand 
roughness. 

For a fluid of constant density, application and integration 
of the continuity equation after 
yield 


W+ = 1 — + 
(1 — d,/d){In(: 30d /ks,)(1 d,/d)|- 


where W*+ = W/(pUd). By considering a force balance 
upon a nonaccelerating fluid element, appropriate substitu- 
tion yields the relation 


AP+ = 0.16{ {In(30d/ks,)\(di/d)| -2 + [In(30d/ks,)(1 — di/d)|-2} 


pU? 


where 


It can be seen that the dimensionless flow parameter and 
pressure drop parameter are functions of the dimensionless 
ratios d/ks,, d,/d and d/ks,. 


Since d/k;, and k;,/ks, are stipu- 


Y 


Diagram of flow system and coordinates | 
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Fig. 2 Distance of maximum velocity line from rough wall 


luted by the channel boundary condition, W*+ and AP* are 
functions of d,/d only. The behavior of each of these func- 
tions is such that as the value of d,/d increases from 0.50 for 
a given set of boundary conditions, each function passes 


through a minimum. The minima for these parametersdo not 
occur at the same value of d,/d. The minimum value for | 


ASP* always occurs at a lower value of d,/d than does the 
minimum value of W*. Indeed the minimum value of Wt 
frequently occurs at values of d,/d > 1.00, which is unreal. 

It is well to reflect at this point that an additional restraint 
upon the system is required if a unique value of d,/d is to be 
arrived at for given boundary conditions. In choosing a 
reasonable restraint, it is possible to advance several condi- 
tions. It can be proposed that the velocity profile will adjust 
so that: (a) AP*+ isa minimum; (b) W* is a maximum, or 
(c) the ratio W*+/AP*+isamaximum. The second restraint can 
be eliminated because the value of W* increases monotonically 
on either side of its minimum point within the range of 
physical interest. 

More reasonable restraints are posed by either the first or 
third condition; in fact, these restraints are essentially 
equivalent because of the slight variation in W+ with change 
of d|/d. The more conservative one is that W*+/AP* possesses 
a maximum in that less distortion of the velocity occurs. 
Consequently, if it is assumed that the flow will arrange itself 
so that this ratio maximizes, the effect of variation of d/ks, 
and k,,/ks, can be examined. The maximum occurs when 


See “9 ae 
= 
is satisfied. 

The effect of unequal roughness upon the value of d,/d is 
shown in Fig. 2. Examination of Fig. 2 shows that unequal 
roughness will produce a substantial skewing of the velocity 
profile. The range of values of k;,/ks, up to 200 are encoun- 
tered in two-phase situations and hence explain Krasiakova’s 
(2) results. 

Fig. 3 shows the relative effect of variations of d/k;, and 
k,/ks, when the variations are due to changes in k;,. Ex- 
amination of this figure shows that any increase in k;, with 
the other boundary conditions remaining constant will always 
cause the flow parameter Wt to decrease and the pressure 
parameter AP* to increase. It is obvious from this figure 
that the effect of variations of d/ks, is severalfold that of the 
term k;,/ks,. Consideration of this figure offers some insight 
into the mechanism responsible in large part for the increased 
pressure drop encountered in horizontal annular two-phase 
flow. 
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Fig. 3 Relation of AP* to W* as affected by system roughness 
ratio and relative roughness of plate 1 


It can be concluded that the above treatment offers a 
method of approximating the alteration in pressure drop and 
weight rate of flow caused by unequal roughness in a duct. 


Nomenclature 
d = distance between plates, ft mn ol 
d; = distance from plate to maximum velocity point, ft 
g = acceleration of gravity, ft/sec? 
vs = sand roughness, ft 
L = length in direction of flow, ft 
AP = pressure drop, lb/ft? 

AP+ = pressure drop parameter, dimensionless, 
U = maximum velocity, ft/sec : 
u = velocity, ft/sec 
Ux = friction velocity, ft/sec 
W = weight rate of flow, lb/sec ft 
W+ = flow rate parameter, W/(pUd), dimensionless 
p = density, lb/ft® 
Subscripts 
1 = lower plate 
2 = upper plate 


References 


1 McManus, H. N., Jr., “An Experimental Investigation of Film 
Characteristics in Horizontal Annular Two-Phase Flow,’ ASME Paper 57-A- 
144, Dec. 1-6, 1957. 


2 Krasiakova, L. I., “‘Some Characteristic Flows of a Two-Phase Mix- 
ture in Horizontal Pipe,’’ Zhurnal Technicheskoi Fiztki, vol. 22, no. 4, 1954, p. 
656. 


145 


he 
er | | 
0.012 
he ust 
m vibe 
n- A 
4 
Is 
= 
le 
0.008 72? 
"sz 100 = 1%, 
1000 
in 
+ 10 
al AP : 
0.004 
ov 40 
= 
100 
= > 
> 1000 
>) te 
> 
4 
| 
5 


W.H, T. LOH! 
Chance Vought Aircraft, Dallas, Texas 


tes 


RBITAL flights are usually analyzed by the diff erential 
d*r 


equations? 
\a) 


(dr) (a6) _ 
de at) \at) ~ 
However, after knowing that the orbit is an ellipse, solu- 


tions may be quickly obtained if the two basic laws of phys- 
ics which govern the mechanics of orbital flight are used 


instead, namely 
Vie 
= — (conservation of energy) [1] 
= (conservation of angular momentum) [2] 


Here V and V, are the orbital velocity and its local Earth 
horizontal component; g is the local gravitational acceler- 
ation; r is the local distance from Earth center, and sub- 
scripts 1 and 2 represent any two points on the orbit. Al- 
though the two simple equations come from ready applica- 
tion of basic physics laws, they also satisfy the differential 
equations, and therefore, offer the same answers. For illus- 
tration, the following problems will be analyzed by using the 
simple equations: 

1 Given burnout altitude 7,, burnout velocity V; and 
burnout angle a; (here V; cos a, = Vi), find orbit, orienta- 
tion and period of revolution. 

2 Given orbit and orientation, find required burnout 
altitude, velocity and angle, so that the orbit produced is the 
one prescribed. 

3 For problems 1 and 2, find orbital velocity and direction 
at all points. 

4 Find the new orbit and its orientation by firing a rocket 
instantaneously when satellite is on its old orbit. 

5 Given two prescribed orbits A and B, transfer from or- 
bit A to orbit B. 

6 Given satellite a on orbit A and satellite b on orbit B, 
have satellite a intercept satellite b. 

7 Find the re-entry from a given orbit (neglect air 
existence). 

8 Find the required amount of impulse, its direction of 
application and amount of propellant required in order to 
solve problems 5, 6 and 7. 

9 Find the minimum amount of propellant (optimum 
trajectory) for solving problems 5, 6 and 7. 

Solutions for each problem will be briefly described below. 
1 Problem 1: Let subscript 1 stand for burnout, sub- 


script 2 stand for perigee or apogee where Ve, = V2. Elimi- 
nating V: from Equations [1 and 2], one obtains - 
>» ‘ 
—gri? = 4/(giri?)? + 2 — ar 

2 
fs = [3] 


= an) 


2 
2 
“ ” 


Here “‘+” indicates perigee r,; ‘““—” indicates apogee rz. 


Received Oct. 13, 1958. 4 4 
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? Assuming a spherical Earth with mass concentrated on Earth 
center and without pertubation effects. This results in an 
elliptical orbit with Earth center situated at one of its focuses. 
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a= — b= 
6 2 


Fig. 1 


Intercept problem 


The semimajor and minor axes can be found from simple 
geometry of the ellipse 


and the elliptical orbit can be drawn according to 


x? y? 

a? b? 
The orientation is obtained by rotating 7; vector about Earth 
center, which is located at one of the two focuses, (ra — rp)/2 
from center of ellipse along the major axis, until it inter- 
sects the ellipse. (This can be done analytically or graphi- 
cally.) The time period of one revolution 7 is simply 

area of ellipse mab 


area swept over per unit time 3rV, 


Note that because of Equation [2], 3rV, may be taken at 
any point on the orbit, say point 1, the burnout point. 

The velocity and direction at any other point on the orbit 
may be obtained directly from Equations [1 and 2] by letting 
subscript 2 stand for the other point, while subscript 1 stands 
for burnout point, which is known. 

2 Problem 2: Problem 1 is to find r, and r, from burnout 
altitude, velocity and direction; problem 2 is the opposite 
problem. It may be solved in a similar but reversed order. 

3 Problem 3: Solved in the last portion of the solution 
to problem 1. 

4 Problem 4: When the additional velocity V, resulting 
from the instantaneous rocket firing, is vectorially added to 
the local orbital velocity, at which the rocket was fired, the 
resultant velocity and its direction, together with its dis- 
tance from Earth center, determine the new orbit in the same 
manner as in the solution to problem 1. 

5 Problem 5: 

(a) Intersecting orbits (coplanar or noncoplanar): A.t the 
intersecting point, the velocities and directions along both orbit 
A and orbit B are known from the solution to problem 3. 
In order to transfer from orbit A to orbit B, a rocket may be 
fired instantaneously at the intersecting point to give the 
satellite an additional velocity V¢ so that the vectorial sum 
of V4 and V¢ is exactly equal to Vg. This problem then re- 
duces to problem 4. 

(b) Nonintersecting orbits (coplanar or noncoplanar): 
An intermediate orbit C may be drawn to intersect both 
orbits A and B, and two instantaneous rocket firings may be 
made at the two intersections. This problem then reduces 
to problem 5(a). 


6 Problem 6: Use problem 5 solution to transfer satellite 


afrom orbit A to B. Determine the time difference At between 
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satellite a on orbit B and satellite 6 on orbit B. Here At is 
equal to the shaded orbital area on Fig. 1, between a and b, 
divided by 3rV,; the value of 3rV, may be taken at any point 
on orbit B. In order to have satellite a intercept satellite 
b on orbit B, it is necessary to find an orbit D which passes 
through point N and has a period equal to rp = tz — At. 
Here rz is the period of orbit B. After elapse of time rz — At, 
satellite b along orbit B and satellite a along orbit D will 
reach point N at the same instant and make the exact inter- 
cept. Using ra, rp from Equation [3], one obtains 


1 
2 Nh 


In this expression 7,, At, rv, gy are known, so the only un- 
known, Vy (for orbit D), can besolved. For minimum propel- 
lant, orbit D could be made tangent to orbit B. This requires 


Here (Vy,/Vwy)ze is known from problem 3 solution. With 
}yand Vy, on orbit D and ry known, orbit D is completely 
determined by problem 1 solution. There are many other 
ways to make the intercept. The above example takes more 
than one rocket firing. Let us now briefly illustrate another 
example where only one instantaneous rocket firing is used: 
At time ¢, satellite a at point # on orbit A and satellite 6 at 
point F on orbit B. In order to make the intercept, it is 
necessary to determine the locations of the corresponding 
pair of points # and F on orbits A and B, respectively, so 
that an orbit D can be determined in such a way that: 

(a) orbit D, with Earth center as its focus, passes through 
both points Z and F, and 

(b) on orbit D, the orbital area between £ and F divided by 
its ($rV,) is exactly equal to rs. The second requirement 
specifies that, after an elapse of a time equal to 7z, satel- 
lite a along orbit D (from £ to F) and satellite b along orbit 
B (from F to F) will both reach point F at the same instant and 
therefore will make the exact intercept (such solutions exist 
in certain cases). 

7 Problem 7: When the existence of air is neglected, the 
problem becomes one of determining an elliptical orbit with 
Earth as its focus that passes through two points, # and F. 
E is the location at which instantaneous rocket firing for 
transfer occurs, and F is the point where landing is desired. 
This reduces to problem 6. 

8 Problem 8: Problem 8 can simply be determined from 
the law of action equal to reaction.’ Action is MpV;, and 
reaction is (Mop — Mp)V. Here Mpand Mare the propellant 
mass being ejected instantaneously and the satellite total 
mass before instantaneous propellant ejection, respectively; 
V is the additional velocity required for the transfer, and 


referred to as the specific impulse “7” of the rocket). 


(Mo — Mp)V = MeV; = Ig. Mp [4] 
Therefore the amount of rocket impulse required is ioe 

Mo — Mp)V = Mo{ 1 ] V 


° If propellant ejection is not instantaneous but takes time dt 
then Me — Mp) = M; dM = —dMp; Equation [4] reduces to 
MdV = —V;dM and results in V = Igo ln [Mo/(Mo — Mp)). 
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\, is the exhaust velocity of the propellant (V ;/go is usually — 


and the amount of propellant required is 


V 
Mp = M, | ——— 


The direction of application, of course, should always be in 
the direction of the additional velocity required to accomplish 
problems 5, 6 and 7. 

9 Problem 9: Whenever and wherever possible and ap- 
plicable, the orbit transfer should always be made in the 
direction tangential with the original orbit, so that the addi- 
tional velocity V required for transfer will be a minimum, 
and consequently, the propellant required will be a minimum. 
For transfer to a smaller orbit or re-entry, the minimum propel- 
lant transfer always occurs at apogee, mainly because there 
the original orbital velocity is a minimum; consequently, the 
additional velocity required to reduce the original orbital 
velocity is a minimum. Similarily for transfer to a larger 
orbit, such as transfer to an orbit with higher apogee or even 
escape (from gravitational field), minimum propellant transfer 
will always be necessary at perigee, mainly because there 
the original orbital velocity is a maximum, consequently the 
additional velocity required to increase the original orbital 
velocity isa minimum. Similarily for the purpose of increas- 
ing the life of a satellite with minimum propellant, the in- 
stantaneous rocket firing should always be made tangentially 
at the apogee, so that the increase of perigee altitude will be 
a Maximum. 

10 Other Problems: Many other satellite orbital me- 
chanics problems not covered here may also be solved in a 
similar manner by the two simple equations used here. Such 
problems include escape mechanics (parabola or hyperbola 


trajectory when V; s 


4/2g:r; is the so-called escape velocity; it is obtained from 
Equation [1] by setting escape conditions: V2 = 0, rz = @. 


This results in: V,; = V 29:71 = V 


Terminal Phase of Satellite Entry 
Into the Earth’s Atmosphere 


ELLIOTT D. KATZEN! 
Ames Research Center, NASA, Moffett Field, Calif. — 


NONLIFTING satellite entering the Earth’s atmosphere 

at a small or zero inclination of the flight path eventually, 

in the terminal phase, reaches a relatively high flight path 

angle. Under these conditions, the equations of motion can 

be simplified to provide a tractable solution for the satellite 

trajectory. This solution is particularly useful since the 

trajectory in the terminal phase is difficult to obtain from less 
restrictive solutions, such as that of Chapman (1).? 

In Chapman’s nomenclature the two component equations 


_ of motion are 


dty dv wD. 
dt? dt g r m — 
dy w D 
dt r m ee [2] 


During the terminal pnase of the trajectory the terms u?/r 


Received Oct. 24, 1958. 
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and uv/r may certainly be neglected. With the further as- 
sumption sin g  —1, and introducing the drag coefficient, 


the equations become 
_a /2 | 
"| 


du cos 


[4] 


du 


so that the two: voneiniiil equations are independent. 

For the vertical component, Linnell (2) gives a solution to 
Equation [3] involving use of a tabulated function of a single 
variable. The solution for the horizontal component is the 
same as that given by Allen and Eggers (3) for the vertical 


component with gravity neglected mation 


where jp is defined from the exponential approximation for the 
variation of atmosphere density 


| 


and the subscript 7 denotes initial velocity and altitude. 

In Fig. 1 are shown flight path angles computed by combin- 
ing the vertical component of velocity from (2) with the hori- 
zontal component from Equation [5]. Comparison is made 
with machine calculations made using more complete equa- 
tions (4), and it is seen that agreement is good. In making 


= exp( —By) 


this comparison, the two methods of computation were 


matched at an altitude of 80,000 ft. ae 


1 Chapman, Dean R., ‘‘An Approximate Analytical Method for Studying 
Entry Into Planetary Atmospheres,’’ NACA TN 4276, 1958. 

2 Linnell, R. D., “Vertical Re-entry Into the Earth’s Atmosphere for 
Both Light and Heavy Bodies,” Jer Propusion, vol. 28, 1958, pp. 329-330. 

3 Allen, H. Julian and Eggers, A. J., Jr., ““A Study of the Motion and 
Aerodynamic Heating of Missiles Entering the Earth’s Atmosphere at High 
Supersonic Speeds,’’ NACA TN 4047, 1957. 

4 Nielsen, Jack N., Goodwin, Frederick K. and Mersman, William A., 
“Three-Dimensional Orbits of Earth Satellites, Including Effects of Earth 
Oblateness and Atmospheric Rotation,’’, NASA Memo 12-4-58A, 1958. 
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- Surveillance of Solid Propellant 


Rockets 
CARL BOYARS! 


Research and Development Department, U. S. Naval 
Propellant Plant, Indian Head, Md. 


HE ORIGINAL meaning of “surveillance’’ for munitions, 

i.e., inspection to assure that the material was safe to re- 
tain and use, has been somewhat broadened. The purpose 
of a rocket surveillance program is to ascertain the length of 
time that rockets of the type under investigation will remain 
in storage and will continue to function after storage and the 
exposure to conditions incident to use, in accordance with the 
requirements of the propulsion system. When the rockets 
are those which have been accepted for service use and are in 
storage, periodic inspection and test may be called “quality 
surveillance.” If the rockets under investigation are the re- 
sults of a development rather than production program, ani 
it is desired to learn what the safe and useful life expectancy is 
under various storage conditions, such a program is often 
called “type life testing.” Studies of the nature and effects of 
the degradation processes, which should be carried out in the 
early stages of research and development of a new rocket, are 
generally considered investigations of “stability.” 

Definition of terms relative to propellant and rocket life is 
required. The safe storage life is the time that a propellant or 
rocket can be retained without undergoing auto-ignition in 
storage. The safe use life is the time that it can be retained 
and still function without undue hazard to the user. The 
useful life is the time that it will meet the ballistic require- 
ments it was designed to meet. 


Deterioration Processes 


A propellant is intrinsically an unstable system, for it must 
be capable of releasing a large quantity of energy in the form 
of hot gas. However, it is also necessary that a propellant be 
stable under all conditions of storage up to the time that it is 
called upon to release its energy. The usual solution to this 
apparent inconsistency in requirements involves taking ad- 
vantage of the activation energy of the degradation reactions. 
The rates of the degradation processes due to this necessary 
instability are temperature dependent in accordance with the 
Arrhenius equation. The results of degradation processes in 
propellants are usually observed as some malfunction, some 
measured change in physical properties, a measurement of rate 
of gas evolution from propellant at constant temperature, 
etc. It is generally incorrect to compare such measurements 
at two temperatures and speak of an “activation energy” 
derived from such data. Activation energy is a term which 
should be restricted to reactions whose mechanism involves a 
single step or a single rate-controlling step. In propellant 
degradation, we have many reactions occurring simultane- 
ously and sequentially, all having their individual activation 
energies. The fallacy of considering the end result of a group 
of such reactions as having an activation energy is dangerous 
because of the tendency to extrapolate data beyond the ex- 
perimental temperature limits by use of this pseudo activa- 
tion energy. 

In addition to the kinetic considerations of instability, 
there are equilibrium factors which are also temperature de- 
pendent. When a propellant is stored under other than con- 
stant temperature conditions, it is subject to phase transi- 
tions, which may also result in undesirable side effects in 
heterogenous systems. One example is the 32 C phase transi- 
tion of ammonium nitrate; the difference in density between 
the rhombic and monoclinic forms of ammonium nitrate re- 
sults in an expansion of the binder, which is not completely 

Received Oct. 21, 1958. 
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ideally it should consist of flight firing tests carried out in suf- 
ficient quantity to give statistically significant data. Un- 
fortunately, such a program is often prohibitively expensive, 
especially for larger missiles. In addition, although flight 
testing will reveal failures, it is not always likely to reveal the 
cause of failure. Static ballistic testing can generally provide 
more information about causes, but this is still quite expen- 
sive for large missiles. Often, delving into the cause of failure 
requires chemical analyses, physical properties testing and 


elastic, and overall growth of the propellant is a consequence. 

The conditions of manufacture of a propellant, e.g., casting, 
molding or extrusion, often result in a product which is under 
stress. On aging, stress relief processes occur resulting in 
change of shape and dimensions. Variation of density with 
val temperature even in the absence of phase transitions can also 
lead to difficulties which are magnified when the propellant 
grain is very large or the tolerances small. 

Another important equilibrium consideration for hetero- 


Pes geneous systems is migration of components. For example, stability testing. Appropriate use of these latter techniques 
se plasticizers like nitroglycerin readily migrate from regions of in the early stages of development can establish and lead to 
of high concentration (actually activity) to low. This can be elimination of potential difficulties. 

in serious where a propellant surface is coated with a plastic Stability testing depends upon selecting for measurement 
he inhibitor to restrict burning or where bipropellant systems a parameter directly related to the degradation rate of the 
he are used consisting of segments of propellants of different com- propellant and, in general, accelerating the rate of degradation 


by increasing the temperature of the propellant. Stability 


position formed into a single unit. Vapor phase migration of 
tests may involve the detection of some decomposition prod- 


plasticizer can affect portions of the rocket system which are 


ty not in direct contact with the propellant, such as O-rings, elec- uct, such as oxides of nitrogen from propellants containing 
= trical insulation and igniter cases. nitrate esters or ammonia from propellants containing am- 
nd Auto-ignition during storage was, historically, the first monium salts or amine linkages. Detection may be visual or 

3 recognized manifestation of instability in propellants. Until olfactory, or may rely upon the use of sensitive indicator 
en the development of propellant grains for rockets and the con- papers. Other stability tests are based upon detection of 
of sequent use of larger webs, this was the only stability problem thermal effects within the propellant, upon measurement of 
he viewed with alarm by the armed services. For rockets, the rates of gas evolution and absorption, or upon the detection 
re useful life and safe use life are generally so much shorter than of fissures developing in small cylinders or rectangular prisms 


of propellant. Tests of physical properties under stress may 


the safe storage life, that the latter is not considered a problem 
also be considered stability testing, although the rate of degra- 


because the propellant is normally discarded while still safe to 


Is 

or store. However, it is quite possible to develop propellants dation is accelerated by stress rather than increased tempera- 
in which have a shorter safe storage life. Auto-ignition is a ture. 

od] phenomenon relatively easy to understand. The overall reac- 


tions in propellant degradation are exothermic. When the 
mass of a propellant and its thermal conductance are such that 
the energy thus released is not dissipated sufficiently to the 
surroundings, the temperature within the propellant con- 
tinues to go up, accelerating the exothermic degradation 
further. The process continues until the propellant ignites or 
detonates. Detonation is expected under such circumstances, 


Accelerated Aging 

Accelerated aging of propellants and rockets is an 
attempt to compress the time scale by storage at elevated 
temperatures so that predictions covering long periods of 
normal storage can be made after a relatively short time of 
testing. Accelerated aging may also be considered to include 


cycling through temperature extremes outside those antici- 


n 
= temperature — self-heating is pated as storage conditions for the rocket and cycling 
is len through temperature limits simulating daily or seasonal 
is eterioration can result in loss of propellant strength by variations at a more rapid rate than nature’s. The validity 


- scission of bonds in the polymer network or by affecting the of conclusions about storage behavior at temperatures far re- 


Wri 


binder-oxidizer interface. Decomposition reactions yielding 
gases within the propellant matrix can build up internal pres- 
sures if gas diffusion is slower than evolution. The eventual 
result of such processes can be deformation and internal fissur- 
ing, with consequent serious effects on burning rate. Weaken- 
ing of bonds between liner (inhibitor) and propellant or liner 
and case cause similar difficulties. Deformation of propel- 
lants may also be caused by scission of polymer bonds under 
stress and reformation in different positions. 

To reduce the rate of deterioration of propellants, purifica- 
tion procedures are introduced into the manufacture of in- 
gredients to eliminate trace impurities that would accelerate 
breakdown. An example is the elaborate purification by boil- 
ing and washing of nitrocellulose. Stabilizers are incorporated 
to react with primary decomposition products and thereby 
inhibit further degradation of the propellant by these prod- 
ucts. Stability testing is used to check on the quality of the 
ingredients and the propellant. 

Although the propellant is usually blamed for any failure of 
a rocket to function properly, the prime cause may often be 
one of the other components. A surveillance program must 
therefore ultimately be carried out on the complete loaded 
rocket, whenever possible, for assurance of reliability. Con- 
sidering the variety of possible modes of failure, it is obvi- 
ous that scaled-down model motors are of limited use in a 
surveillance program because of the impossibility of scaling all 


Stability Testing 


Because a surveillance program is meant to provide infor- 
mation on the probability of proper functioning of a rocket, 


Fesruary 1959 


moved from those of the accelerated aging is subject to criti- 
cism as the discussion of “activation energy” above points out. 
Even less valid is the practice of drawing conclusions with re- 
spect to behavior at various temperatures from tests carried 
out at only one temperature. 

The egg provides an example of the hazards of extrapolation 
beyond temperature limits. At around 100 C, a hard-boiled 
egg is produced in minutes. Activation energies have even 
been computed. Yet at somewhat more moderate tempera- 
tures, the product is a baby chick rather than solidified pro- 
tein. Presumably, one could compute an activation energy 
for the production of a chicken from an egg, but the results 
of storage of eggs at normal ambient temperatures are 
well-known. 

Thus accelerated cycling aging is subject to the criticism 
that (a) at temperature extremes beyond those anticipated in 
actual use, deleterious phase changes may occur which would 
never occur otherwise, and (b) rapid cycling subjects the 
system to severe thermal shock, which could result in effects 
never encountered in ordinary storage. 

Nevertheless, accelerated aging is a necessity in order to 
forecast, even with limited accuracy, the probable future be- 
havior of the rocket. When rocket development reaches the 
stage of late pilot production, accelerated aging programs 
under controlled temperature conditions must be instituted. 
Representative samples are randomly selected, stored and, 
at regularly scheduled intervals, withdrawn for test firings, 
visual and X-ray examination, chemical analyses, determina- 
tion of physical properties, etc. 

Accelerated aging programs have proved their value in the 
past. Such a program on rockets showed that, prior to a seri- 
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ous increase in malfunction frequency, certain phenomena 
occur. The accelerated aging program indicated that ma- 
terial in service, while still useful, would give a signal that 
closer observation and monitoring was necessary. In the case 
of another rocket, such a program revealed in a few weeks 
that the generation of hydrogen gas within the metal igniter 
case resulted in a dangerous defect. The igniter was therefore 
redesigned and replaced. 


Long Term Aging 

Long term aging studies generally consist of routine quality 
surveillance testing of rocket motors stored in various parts 
of the world. Rocket motors are selected from the storage in- 
ventory so that the collected sample is representative of the 
entire population of a particular rocket motor. Results of 
ballistic and other tests are analyzed, using the proper statis- 
tical techniques. 

Useful information is gained from such programs. For ex- 
ample, certain air-to-ground rockets have been found to be still 
safe and serviceable at a time when the material might other- 
wise have been discarded because of age. Quality surveillance 
of still another type of rocket revealed an unsafe condition in 
which the round blew up near the end of burning. The reason 
for this behavior was found to be shrinkage of the grain 
diameter due to a change in the inhibitor. A reworking pro- 
cedure has been developed to salvage these rockets. 

High temperature storage is not always accelerated aging. 
In cases where propellant or gas generator systems are sub- 
jected to supersonic heating in use, high temperature storage 
and testing designed to simulate use conditions is more prop- 
erly considered as long term aging rather than accelerated 
aging. In addition, if a propellant is to be subjected to 
stresses, such as those due to flight acceleration, in use, a 
proper ballistic testing program carried out as part of an 
aging study must simulate these conditions. Pe ne 


Determination of Elements of an Elliptic 
Orbit From the Orbital Velocity Vector 


H. MUNICK! 


HIS paper considers a space vehicle which is traveling in 
a certain elliptical orbit. Suddenly an impulse is applied, 
causing the instantaneous velocity to change in magnitude 
and direction. The space vehicle goes into a new elliptical 
orbit. The elements of the resulting elliptical orbit are 
determined in this paper. 

The problem of finding the elements of an elliptic orbit, 
given the orbital velocity vector has been investigated by 
Vertregt (1),? and by Schiitte (2). The paper by Vertregt 
gives an interesting, detailed solution of the entire problem. 
The following discussion gives an alternate solution to this 
problem using a polar coordinate system. 

The velocity vector is at a distance r from the center of the 
Earth (see Fig. 1). From the magnitude of the velocity 
vector, and from 7, we get 


From Kepler’s law we can write 


(r6)r = /R*Gra(1 — 


Received Dec. 8, 1958. 
1 Research Mathematician, Dynamics Group. 
? Numbers in parentheses indicate References at end of ‘Paper. 
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Fig. 1 Elliptical orbit 


Considering the component of the velocity vector pe:- 
pendicular to the radius vector 


ré = vsin B [3] 


From Equations [1, 2 and 3] the eccentricity follows 
v? sin? Br? 


- 
R°Gya 


Finally, @, the angular displacement with a, known , is give: 
by both Vertregt (1) and Schiitte (2) 


1+ (1 —f) tan? B 
The parameter f is given by f = 2 — (r/a). ww ABs 
Nomenclature 
r = distance of vehicle from the center of the planet 3, 
a = semimajor axis of ellipse _ 
e = eccentricity of ellipse ; 
R = radius of the planet bry o 
G = acceleration of gravity at the surface ofthe planet = ; 
@ = angular displacement 
8 = angle between radius vector and orbital velocity vector 
v = magnitude of orbital velocity vector 
References 


1 Vertregt, M., “Die Bahnbestimmung aus Dem Vektor der Bahn- 
geschwindigkeit,”’ Astronautica Acta, vol. IV, fase. 2, 1958, pp. 1385-137. 

2 Schiitte, K., ‘‘Die Bahnbestimmung aus Dem Vektor der Bahn- 
geschwindigkeit und der Einfluss einer Aenderung Desselben auf die Bahn- 
elemente,”’ in “Probleme der Weltraumforschung,’"’ Laubscher and Co. 
1953, p. 89. 


Seale of Separation Phenomena in 
Liquids Under Conditions of Nearly 


Free Fall 


Dougias Aircraft Co., Inc., Santa Monica, Calif. 


An estimate is given of the relative importance of the 
capillary (adhesion and cohesion) forces as compared to 


_ the gravitational or the (equivalent) inertial forces acting 


on a liquid. 


The results indicate that extremely small 


_ lead factors are necessary for the occurrence of capillary 


phenomena, such as conglomeration of liquids in spherical 
masses or their attachment to the walls of the con- 


Received October 29, 1958. 
‘Group Physicist, Theoretical Physics and Applied Math- 
ematics. 
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tainer. It is also shown that the liquid has to be in a 
gravity-free or unaccelerated environment for a compara- 
tively long time before such phenomena develop. The 
influence of viscosity upon the development of capillary 
effect was found to be negligible except for liquid masses of 
microscopic size. 


HEN a liquid is removed from the action of gravity, or 

more generally, when the inertial forces (qualitatively 
equivalent to gravity) acting upon the liquid are of suf- 
ficiently low intensity, the intermolecular forces (cohesion or 
surface tension and adhesion) will prevail in determining the 
shape and the free surface of the liquid. Such phenomena 
must obviously be taken into account in the design and opera- 
tion of missiles and space ships; in addition, their physiologi- 
cal consequences on the human body are obviously of great 
in.portance for aero- and space medicine. Quantitative results 
on this subject (if existing at all) did not apparently find their 
way into the more readily available literature. Therefore, the 
present elementary analysis was undertaken for the purpose 
o! providing a survey of the expected scale of such phenomena. 


Inception of Capillary Phenomena 

A basic question which occurs in connection with the prob- 
lems described above is the range of accelerations for which 
capillary phenomena become important. For this purpose, 
consider a portion of a liquid of density p, surface tension a, 
occupying a volume of linear dimensions LZ. The intensity of 
the capillary forces? acting upon it will be of the order of aL, 
whereas, if this portion of liquid is accelerated with an ac- 
ccleration ng, it will be subjected to a force of inertia of in- 
tensity having the order of ngpL*. The (nondimensional) 
ratio 

= 

measures the relative importance of the capillary forces in an 
accelerated liquid. It can therefore be expected that capillary 
phenomena will have to be considered when St > 1, i.e., for 
liquids subjected to accelerations corresponding to values of n 
smaller than 


a/pngL? 


o/pgL* [1] 


ne = 


The same conclusions can also be reached by considering an 
actual physical case, e.g., the capillary rise (or depression) h of 
a liquid in a tube of radius r, subjected to a longitudinal ac- 
celeration of intensity ng. According to Jurin’s law, we have 


h = o/pngr 


Detachment or attachment of liquids from or to the walls of 
the tube can be anticipated when n is small enough so that h > 
r,i.e., forn < o/pgrh. 

The surface tension of most liquids of interest in missile or 
space ship engineering is of the order of magnitude of 25 
dyne/cm, with the exception of water, whose tension is about 
75 dyne/cm (the exact value depending upon its purity). The 
adhesion (or interfacial tension) between liquids and solid 
surfaces is of the same order of magnitude. We shall accord- 
ingly assume oa to be of the order of magnitude of 50 dyne/cm. 
Taking further p & 1 gr/em?, and g 1000 cm/sec?, we ob- 
tain from Equation [1] 


nx = 0.025/L? (ZL in em) = 3.9 X 1073/L? (L in in.) 


lor example, capillary phenomena in tanks of the size of the 
order of 10 in. will become important if their acceleration falls 
below 8 X 1075 g. 


2 We shall denote with this term the cohesive (or surface ten- 
sion) as well as the adhesive (or interfacial tension) forces. For 
the purposes of a preliminary investigation, distinction between 
these forces (all of which originate from intermolecular actions, 
and are accordingly of the same order of magnitude) appears un- 
necessary. 
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Time Scale of Separation of Phenomena 


Another important question to be answered is the minimum 
duration for which a liquid is to be exposed to conditions of 
weightlessness or near weightlessness in order for conglomera- 
tion or separation phenomena to occur. For this purpose, 
consider a certain volume of liquid which is neither accelerated 
nor under the action of any gravitational field, i.e., a freely 
falling liquid. If the shape of the surface boundary of the 
liquid does not initially possess the equilibrium configuration 
corresponding to a minimum of potential (surface) energy, 
the combined action of the forces of surface tension and ad- 
hesion will cause the surface boundary to assume such a shape. 
(In cases in which the surface tension forces predominate, the 
configuration of equilibrium can be expected to be a sphere, 
this figure possessing the smallest surface bounding a given 
volume.) Indicating as before with L the linear dimensions 
of the liquid body, its mass will be of the order of pZ*. Since 
the forces causing its deformation are of the order of oL, the 
deformation will proceed with an acceleration > aay 


a/pL* [2] 


If the discrepancy*® of the initial shape of the surface from 
that of the equilibrium configuration is of the order of L, the 
time required for the deformation to take place will be of the 
order of 


a= 


Taking as before 50 dyne/em, r L’/:/7 sec. For 
example, the separation of a liquid from the walls of a tank 
having linear dimensions of the order of 10 in. would take 
place about 18 see after removal of the force of gravity. _ 


Influence of Viscosity 3 

The fact that the retarding effects of viscosity have been 
neglected merits a few words of justification. The order of 
magnitude of the speed v of flow is given by V/al, so that we 
obtain from [2] 


Vo/pL [3] 


Indicating with v the kinetic viscosity, the density of the rate 
of dissipation of energy due to viscous friction will be given 
by pyv?/L?, so that for the rate of the total dissipation of 
energy we shall have, using Equation [3] 

The rate at which the surface tension forces perform work can 


be estimated by the expression 


U = (d/dt)\(oL?) 


That is, again using [3] 


U = V0'L/p 
From [4 and 5] we obtain the estimate = 7 


U/w = (1/v)VoL/p 


for the relative importance of capillary as compared to viscous 
processes. Since for liquids (other than oil) vy ~ 10~? em?/sec, 
and ~/a/p can be taken to be of the order of unity (except for 
water, for which this quantity is about 8), it follows that the 
effects of viscosity can be neglected as long as one deals with 
liquid masses whose linear dimensions are much greater than 
1 yu. This means that viscosity can be disregarded as long as 
macroscopic separation effects are dealt with. 


3 This discrepancy can be measured, e.g., by the differences be- 
tween the maximum dimensions of the initial configuration and 
a dimeusion of the configuration of equilibrium. 
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Turbojet aircraft engine with thrust aug- 
mentation (2,857,740). J. R. Hall and 
D. E. Singelmann, Snyder, N. Y., assign- 
ors to Bell Aircraft Corp. 

Jet Propulsion engine with primary and 
auxiliary combustion chambers. Valves 
cause a gas generator to drive the auxiliary 
turbine to in turn drive a compressor while 
the discharge air is directed to flow into 
the auxiliary combustion chamber. 
Temperature control device for an engine 
(2,857,742). J. A. Drake, Van Nuys, 
Calif., assignor to Marquardt Aircraft Co. 

Heat added to the inlet of a combustion 
chamber accelerates the hot fluid to sonic 
velocity, venting it through a choked out- 
let Adjustment means control the actu- 
ating means responsive only to a total tem- 
perature and to static and total pressures 
of the fluid at the inlet. 

Balloon having a reinforced seal (2,858,- 

090). O. C. Winzen and V. H. Winzen, 
Mendota Townshi Dakota County, 
Minn., assignors to Wi inzen Research, Inc. 

Reinforcing and load carrying member of 

flexible thermoplastic material having a 
part heat-sealed to one side of a fin formed 
by seams and gores turned outward from 
the envelope. 
Fuel tank with expellent bag (2,859,808). 
R. Youngquist, D. F. Ferris and H. S. Bell 
Sr., Whippany, N. J., assignors to Thiokol 
Chemica Corp. 


Flexible bag mounted in a fuel tank to 
seal off the fuel from a pressurizing gas 
acting on the opposite side of the bag. 
Pressure extends the bag to force fuel 
from the tank through a discharge port, 
the bag being sealed off from the fuel until 
pressurization is effected. 


Hypergolic fuel and method of propelling 
rockets (2,859,587). C. R. Scott and A. 
L. Ayers, Idaho Falls, Idaho, assignors to 
Phillips Petroleum Co. 

Method of applying immediate thrust 
to a mass. Separate streams of an oxi- 
dant and fuel are introduced into contact 
with each other in a combustion chamber 

a spontaneous combustion. 
Fuel consists of an oxidizer and tert-butyl 
oo at a ratio of not more than 1 to 

We 


Flameholder structure for ramjet combus- 
tor (2,859,588). H. I. Wilson, Beacon, 
N. Y., assignor to The Texas Co. 
Plurality of holder elements in the shape 
of hollow truncated cones with bases open 
downstream. A splash plate, upon which 


Epriror’s Note: Patents listed above 
were selected from the Official Gazette of 
the U.S. Patent Office. Printed copies of 
patents may be obtained from the Com- 
missioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. 
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fuel impinges, disperses the fuel to provide 
a localized fuel-air mixture downstream. 
Fuel metering device for ramjet engines 
(2,859,589). D. W. Gunnarson, Seattle, 
Wash.., assignor to Boeing Airplane Co. 

Air pressure from a pitot tube outside 
any shock wave of an aircraft maintains a 
constant ratio of fuel and mass airflow into 
the engine throughout supersonic speed 
range. 


Apparatus for causing intermittent com- 
bustion of fuel in a chamber (2,860,484). 
P. Schmidt, Munich, Germany. 


caning fuel valve.. 


combustion chamber 
press button 
starting air~ “starting fuel 


Charges of gaseous fuel formed in peri- 

odic sequence in a combustion chamber 
having a controlled inlet and outlet. 
Energy of the gases can be used selectively 
as a pressure or driving force, and for 
heat energy. 
Speed control system for gas turbine en- 
gines (2,863,283). H.R. Schmider, E. C. 
Palmenberg and W. J. Dietz Jr., Hillside, 
N. J., assignors to Bendix Aviation Corp. 

Servo motor connected to a valve on 
the combustion chamber for adjustment 
of the chamber and exhaust gas nozzle to 
vary the impingement of exhaust gases on 
the turbine. 

Hydraulic servo system (2,863,284). R. 
Scheib Jr. and R. N. Cronemeyer, Sea- 
ford, N. Y., assignors to Sperry Rand 


2,800,484 


rp. 

for stabilizing a controlled 
member from a controlling member by 
means of a hydraulic motor. A signal 
generating means is electrically coupled 
to a solenoid-actuated reversing valve for 
actuating the valve upon change in direc- 
tion of motion of the controlling member. 
Rocket structure (2,862,447). G. A. 
Lyon, Detroit, Mich. 

Shell secured upon the top bulkhead of 
a motor casing, shell and bulkhead having 
engaging thrust shoulders. The shell 
shoulder comprises an indented annular 
rib bearing against the margins of the 
bulkhead. 
Structure for diverting gases of high ve- 
locity (2,860,663). W.J. Kroeger and G. 
S. Schecter, Sommerton, Pa. 


Hollow chamber used as the diverting 
structure for a sub-supersonic or super- 
sonic velocity fluid stream. The stream 
enters through at least one duct and exits by 
means of at least one other duct angularly 
disposed relative to the original path of the 
stream. 


George F. McLaughlin, Contributor 


Water injection system for gas turbine 
power plant (2,863,282). B. N. Toreil, 
Wethersfield, Conn., assignor to United 
Aircraft Corp. 

Valve for controlling the injection of 
coolant into a power plant. Position of 
the valve is determined by a temperature 
sensing element, and the pressure drop 
across the valve is determined by a re;:- 
ulating valve which meters the coolant «s 
a function of inlet temperature and am.- 
bient atmospheric temperature. 
Aerodyne remote control system (2,863,- 
619). L. Servanty, Paris, France, assi.- 
nor to SNCA du Sud-Ouest. 

Detecting means on board an aerodyie 
for measuring angular velocity around 
each of the axes of roll, pitch and yaw. 
Measurements are transformed to ele: 
trical effects and multiplex radio transmi' 
ters emit the effects to a remote radio rv- 
ceiver. 

Jet propelled aircraft (2,863,620). A. ©. 
R. Vautier, Paris, France, assignor to 
Sud-Aviation SNCA. 

Monoplane with two jet engines on 
supporting structures positioned laterally 
outside the fuselage at the same level as 
the wings and horizontal tail surfaces. 

Fuel supply means for combustion appara- 
tus (2,861,424). P. L. Jurisich, Torrance, 
_ assignor to Douglas Aircraft Co 
ne. 


2,861,494 


An injector of fuel and air located in the 
upstream zone of the casing, so profiled 
and arranged as to have low drag to pro 
duce negligible wakes. The injector dis- 
charges a multiplicity of fine fuel-air 
streams that fill the cross section of the 
casing with a continuous, uniform sheet of 
combustible fluid. 

Vertical and horizontal flight aircraft 
(2,863,621). J. W. Davis, El Paso, Texas 

Reaction motor with a nozzle facing 

rearward for ejection of exhaust gases 
horizontally within the fuselage. A tur- 
bine is connected to the exhaust for rota- 
tional movement about the vertical axis 
responsive to impingement of gases ejected 
from the nozzle. 
Means for control of the air inlet opening 
of a jet propulsion unit (2,864,236). K. 
Toure and J. R. Lauzely, Paris, France, 
assignors to SNECMA 

Unit designed for supersonic speeds. 
Flow of atmospheric air into the casing 
is controlled by modifying the effective 
leading edge, and projecting fuel forward 
of the unit and outside the casing to form 
a continuous fluid wall around the leading 
edge. 
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BAKER & ADAMSON® 


Storable and Powerful 


Rocket Fuel Oxidizer 


Chlorine Trifluoride (CIF, 


‘Much is a of Chlorine Trifluoride as a rocket fuel 


oxidizer. That’s because it combines good handling and 


storage properties with very high performance. To be 
specific: 


CTF is not difficult to handle! Boiling point of 
Chlorine Trifluoride is 53.15°F. Freezing point is 


-—105.4°F. Vapor pressure at 140°F is 80 psia. 


CTF may be stored conveniently! Chlorine Tri- 
fluoride is storable over a very wide range of temperatures. 
You can count on years of storage life! 


CTF has excellent stability! Shock resistance of 
pure CTF is very high. And it is thermally stable to high 
temperatures. 


CTF offers high performance! The high density of 
Chlorine Trifluoride (1.825 at 68°F) leads to outstanding 
density impulse values. An important plus—CTF is hyper- 
golic with hydrogenic fuels over a wide range of pressures 
and temperatures. 


Fluorine and other Fluorine-Based Oxidizers 


CTF is one of several fluorine-based chemicals produced 
by General Chemical which are now considered as excel- 
lent rocket-fuel oxidizers for various missions. Another is 
Bromine Pentafluoride which shares many of CTF’s desir- 
able physical and chemical properties. General Chemical, 
the sole producer of elemental liquid fluorine, has exten- 
sive technical data available on Fluorine and Halogen 
Fluorides. Write today for further information on these 
high-energy oxidizers. 


Products 


Typical test stand of Rocketdyne, a division of North American Aviation, Ine. 


* 
GENERAL CHEMICAL DIVISION 


40 Rector Street, New York 6, N. Y. 
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Philadelphia, March 
New York, 1957. 


Method and apparatus for rail transporta- 
tion (2,864,318). H. A. Toulmin Jr., 
Dayton, Ohio, assignor to The Common- 
wealth Engineering Co. of Ohio. 

Land vehicle for traveling upon spaced 
rails. Runners on trucks engage the 
rails and slidably support the vehicle. 
Flanges limit the vertical movement of 
the runners. Jet propulsion units propel 
the vehicle along the rails, and airfoils 
elevate the vehicle above the rails when a 
predetermined speed is reached 
Pressure ratio senser (2,864,393). J. A. 
Drake, Van Nuys, Calif., assignor to Mar- 
quardt Aircraft Co. 

Control device having an output pro- 
portioned to the ratio between two com- 
pressible fluid pressures. A convergent- 
divergent nozzle receives the higher of the 
two pressures, and has sonic flow at its 
throat. Supply pressure may be modified 
to obtain a regulating pressure propor- 
tioned to the ratio of the two pressures. 
Method of generating electrical power 
from solar energy (2,864,879). H. A. 
Toolmin Jr., Dayton, Ohio assignor to 
Basic Research Corp. 

Solar generator consisting of thermo- 
electric elements forming a hollow cylin- 
der. A lens system concentrates heat 
rays of the sun into a beam within the 
cylinder. Means at the cylinder bottom 
reflect the beam against the inner wall to 
further energize the elements. 

Target tracker (2,864,959). J. J. Nas- 
tronero, Port Washington, N. Y 

Device for tracking intermittent targets. 
A substitutive signal supplied by a mem- 
ory means impels the scanning means, 
during periods when the target radiation 
is not received, to continue along the 
tracking course it was following prior to 
the cessation of radiation. 

Variable nozzle for jet engine (2,865,165). 
R. Kress, La Mesa, Calif., assignor to 
Solar Aircraft Co. 

Swinging of interlocked channels to- 
ward or away from the center of the nozzle 
opening causes the size of the opening to 
decrease or increase. The shape of the 
opening remains circular for any size. 
Reaction motor-powered hydraulic and 


electric power unit (2,865,168). J. H. 
Zillman and H. M. Steele Jr., Pacific 
Calif., assignors of one-third 
each to A. Deutsch, L. Deutsch and B. 
Weingart. 

Bearing supporting a plenum chamber 
for rotation about a fixed axis, and a jet 
nozzle in the wall of the combustion cham- 
ber disposed to cause rotation of the cham- 
ber. speed sensitive braking means is 
rotatively coupled to the chamber to reg- 
ulate its speed of rotation. 


Jet exhaust nozzle (2,865,169). G. F. 
Hausmann, Glastonbury, Conn., assignor 
to United Aircraft Corp. 

Annular sleeve surrounding a casing 
with sections which are perforate and 
imperforate. Perforate sections comprise 
two sets of passageways directed down- 
stream and discharging to the free air- 
stream for silencing the exhaust flow. 
Another set is directed upstream, dis- 
charging to the free airstream for reversing 
the exhaust flow. 


Method of jet deflection (2,865,170). M. 
Kadosch, Paris, France, assignor to 
SNECMA. 

Guides positioned outside of the normal 
undeflected path of the jet, and control- 
lable means in the tail pipe increase the 
temperature of an eccentric part of the 
stream. Dissymmetrical heating estab- 
lishes a transverse temperature gradient 
across the stream. 


Rotary reaction engine (2,865,172). T. 
Bodde, Philadelphia, Pa 

Means for forcing a mixture of air and 

slugs of water into a combustion chamber 
at a position adjacent to the radial outer 
end of the chamber. Separate means 
force fuel into the chamber at a position 
located radially inward of the point of en- 
try of the air and water. 
Screech reduction in combustion chambers 
(2,865,174). A. W. Blackman Jr. and 
G. D. Lewis, Manchester, Conn., assignors 
to United Aircraft Corp. 

Perforated wall in a burner downstream 
of an imperforate burner, and forming a 
tuned chamber. Fuel and air pass from 
outside to the inside of the burner veatiaeais 
the perforations in the side wall. : 


Jet propulsion system for maintaining air- 
craft in a vertical attitude (2,865,579). 
G. Caillette, Nevitty-Sur-Seine, France. 

Aerodyne with first jet propulsion 
means positioned to develop thrust to 
coincide with the longitudinal axis at all 
times. Second propulsion means develop 
thrust to coincide with the vertical 
axis. The second means may be varied in 
accordance with the magnitude of the 
angle between the two axes. 


Oil cooling and drag reducing system 
(2,865,580). E. V. Marshall, Dallas, 
Texas, assignor to Chance Vought Air- 
craft, Inc 

Fuel-oil heat exchanger transferring 
heat from the oil to the fuel. A valve re 
sponsive to a Mach meter reduces aero 
dynamic drag caused by ram air on the 
air supply duct inlet system to provide oi! 
cooling and drag reduction. 


Balloon gondola (2,865,581). H. E. Froe 
lich, New Brighton, Minn., assignor t 
General Mills, Inc. 

Gondola for high altitude manned flight 
comprised of a rigid body for supporting a 
passenger. Suspension lines from th¢ 
lifting cell to the gondola permit changing 
the angular orientation between the bod) 
and cell in a vertical plane in response to 
movement of a passenger between selec- 
tive positions. 


Rocket motor liquid propellant combina- 
tion (2,865,727). J. W. DeDapper, Los 
Angeles, Calif., assignor to North 
American Aviation, Inc. 

Fuel comprised of not less than 80 per 
cent by weight turpentine and the re- 
mainder an ethereal liquid in which up to 
4.5 per cent of lithium aluminum hydride 
is dissolved, and which is miscible with 


turpentine. 


Autopilot steering system (2,866,146). 
C. F. Rodriguez Jr., Los Angeles, Calif., 
assignor to Aerojet-General Corp. 
Precession of the gyro rotor correspond- 
ingly changes capacities of four conden- 
ser plates fixed at equal angular distances 
from each other. Changes in the capaci- 
tances produce an output voltage at the 


frequency discriminators. 


Advances in Nuclear Engineering, vol. I, 
Proceedings of the Second Nuclear 
Engineering end Science Conference. 

1957, Part 1. 


se 7 Reviewed by J. P. FRANKEL 
University of California 


The six nations of Western Europe that 
comprise the European Community for 
Coal and Steel have estimated that they 
will need approximately 24 million kilo- 
watts of additional electrical producing 
capability by 1967, and that of these, 15 
million, or about 60 per cent, will be atomic 
installations. 

At the same time, the most optimistic 
estimates of the need for atomic energy in 
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a+ 2-2 4% 


Ali Bulent Cambel, Northwestern University, Associate Editor 


the United States place the fraction of new 
1967 capability that will be nuclear at less 
than 5 per cent. 

From these estimates it is safe to con- 
clude that the need for rapid atomic en- 
ergy development in the United States 
is primarily for the foreign market, a fact 
that is amply brought out in the lead-off 
papers of the Proceedings of the Second 
Nuclear Engineering and Science Con- 
ference, held in Philadelphia in March of 
1957. (Incidentally, this is at least the 
fifth such conference, although the second 
to be sponsored by the same set of Socie- 
ties. ) 

Since the Conference was held, the Sput- 
nik era dawned, and engineers and scien- 
tists, as well as the rest al the thinking 


world, learned the tremendous propaganda 
and psychological value of large-scale 
scientific achievement. This, added to 
the fact that the only early market for any 
American atomic industry is the foreign 
market, makes it obvious that the real 
motivation for governmental support of 
atomic energy development must be from 
the fact that atomic energy, as satellites, is 
instrument of foreign policy. 

For an understanding of the bases of 
such conclusions there is no better ref- 
erence than the papers contained in these 
proceedings, especially the excellent report 
by Nelson and Keagy of Battelle on ‘The 
Economic Background for the Competitive 
Development of Nuclear Power’’ starting 
on p. 210. (Continued on page 158.) 
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OF A SERIES 


The “tuning out” of excessive shake and vibration 
by Oldsmobile engineers produces a comfortable, 
balanced ride that adds thousands of miles to the 
life of an automobile. 

One of the most critical areas of engineering in today’s 
automobile is “ride”. It is critical because an unsatis- 
factory ride means an unsatisfactory automobile. To pro- 
duce an over-all balanced and smooth ride, free from 
harshness and fatiguing vibration, Oldsmobile engineers 
begin the complex task of “tuning” the car in the early 
stages of a new model program. Not only is ride important 
from the comfort standpoint, but an improperly “tuned” 
car can literally shake itself apart after several thousand 
miles, 

The tuning operation is a series of intricate tests that 
determine a car’s “shake” characteristics—where and how 
much the metal bends and twists. To produce beaming 
and torsional moments, a mechanical oscillator is at- 
tached to the frame and vibrates the car in a frequency 
range of 714 to 15 cycles per second. To measure the dis- 
placement of the metal, a velocity pick-up is attached 


OLDSMOBILE > 


Fesruary 1959 


directly to the area under study. As the metal vibrates, 
a signal is produced by the pick-up and is fed to a vibra- 
tion meter where it is integrated. The resulting signal is 
then transmitted to an X-Y plotter that instantly converts 
it into a continuous magnitude-vs.-frequency trace. 


With this valuable information, refining can begin by al- 
tering the structure of the various component parts. A 
complex network of infinite variation must be analyzed 
intensively to produce the mark of quality that stamps 
every Oldsmobile. 


Over the years, Oldsmobile’s reputation for quality manu- 
facturing and precision engineering has grown, step by 
step, until today it is a car of recognized distinction—in 
a class by itself. Oldsmobile’s durability and long service 
life is further attested to by its continued leadership in 
resale value. You owe it to yourself to first examine, then 
test-drive, a truly outstanding automobile—the 1959 Olds- 
mobile. Visit your Local Authorized Oldsmobile Quality 
Dealer as soon as possible. 


OLDSMOBILE DIVISION, GENERAL MOTORS CORPORATION 


Pioneer in Progressive Engineering 
... Famous for Quality Manufacturing 
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SOLAR SAILING: Space travel with the aid of solar radiation 

pressure —an area of advanced research at Lockheed. Vehicle 

would employ a sail that would be raised and lowered in flight. 

The artist has depicted Magellan’s ship “Trinidad” to symboliz = i 

man’s great voyages of discovery. 
Lockheed Missile Systems Division is engaged in all fields of “4 
missile and space technology —from concept to operation. a 
Advanced research and development programs include—man in space: . a 
space communications; electronics; ionic propulsion; nuclear and solar 
propulsion; magnetohydrodynamics; computer development; oceanography; 
flight sciences; materials and processes; human engineering; electromagnetic 
wave propagation and radiation; and operations research and analysis. 


The successful completion of programs such as these not only encompasses 
the sum of man’s knowledge in many fields, but requires a bold and 
imaginative approach in areas where only theory now exists. _— 
The Missile Systems Division programs reach far into the future. 7 7 
It is a rewarding future which men of outstanding talent and inquiring 
mind are invited to share. Write: Research and Development Staff, 

Dept. B-25, 962 W. El Camino Real, Sunnyvale, California, or 7701 
Woodley Avenue, Van Nuys, California. For the convenience of 

those living in the East or Midwest, offices are maintained at 

Suite 745, 405 Lexington Avenue, New York 17, and at 


Suite 300, 840 N. Michigan Avenue, Chicago 11. 


advancement of the art of missiles and astronautics.” 


NATIONAL MISSILE INDUSTRY CONFERENCE AWARD 


Lockheed / 


MISSILE SYSTEMS DIVISION 


SUNNYVALE, PALO ALTO, VAN NUYS, SANTA CRUZ, SANTA MARIA, CALIFORNIA 
CAPE CANAVERAL, FLORIDA - ALAMOGORDO, NEW MEXICO 
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It is interesting to consider the apparent Fluid Dynamics and Heat Transfer, by 


growth of industry participation in the James G. Knudsen and Donald L. 
technical papers of these proceedings. Katz, McGraw-Hill Book Co., Inc., 
In 1953, at the Berkeley Conference, as New York, 1958, 576 pp. $12.50. 


far as this reviewer can remember, all the Reviewed by W. T. Snyper 

_ technical work reported was paid for Se North Carolina State College 

with government money; in the present _ nS 

proceedings, 19 out of 50 technical papers =— The present edition of ‘Fluid Dynamics 
were not obviously supported by govern- and Heat Transfer’ is an enlarged and re- 
ment funds, and therefore presumably vised version of the authors’ book which 
represent outlays of private funds. But originally appeared in 1954 as Bulletin 37 
what kinds of work are reported by indus- of the Engineering Research Institute at 
ew exceptions ese papers are Oo e of the book, as state y the authors, is 
Survey, or State of the Art class involving “to present the fundamentals of fluid dy- 


no capital expenditure. Practically all namics which are basic to an understanding 
the experimental data come from the of convection heat transfer.’’ The book is 
government laboratories, or from indus- divided into three parts with the following 
trial firms spending government money. headings: Part One (44 pp.) Basic 


It is to be noted, however, that the work Equations and Flow of Nonviscous Fluids; 
reported in these proceedings was per- Part Two (248 pp.) The Flow of Viscous 
formed in 1956. Presumably, later pro- = Fluids; Part Three (122 pp.) Convection 


ceedings will show a marked increase in Heat Transfer. 
expenditure of private industry funds for Included also is an appendix with 152 
research. problems, definitions of the commonly 
This volume of the Proceedings, in ad- used vector quantities and a brief but 
_ dition to the general papers, contains five clear discussion of complex variables and 
a sections of technical papers: Fuel Cycles, conformal mapping as applied to the the- 

Plant Containment Concepts and Design, ory of ideal fluid flow. 

=A _ Plant Components, Waste Disposal and In the first part of the book after pre- 
_ Protection and Safety Measures. senting the fundamental definitions es- 
__‘The technical papers are uniformly ex- sential for any fluid dynamics develop- 
cellent, being fairly concise, yet packed ment, the authors develop the fundamen- 
_ with valuable information. Every li- tal conservation equations using a dif- 
brary, in and out of the government, ferential derivation. The derivations are 


should have this for reference. standard and are presented in sufficient 
detail to be followed by most students 
familiar with the Taylor series expansion. 

In the last chapter of Part One, the classi- 


SEZ pp. $9. cal equations of ideal fluid flow theory are 
re Reviewed by S. F. SINGER presented including the concepts of veloc- 
3 University of Maryland ity potential and steam function. Sev- 
___ The title of this book may be mislead- eral examples of two-dimensional, in- 
ing. Natural Aerodynamics is not op- compressible, irrotational flows are pre- 
posed to unnatural aerodynamics but re- sented. 
fers rather to aerodynamic problems which Part Two of the book discusses thor- 
occur in nature, particularly the aerody- oughly Poiseiulle flow, laminar flow in 
namics of meteorology. The book is annuli, and laminar flow in noncircular 
written on an intermediate level and can conduits. Friction factor is defined, and 
therefore be read with profit by the non- friction factor expressions for several lam- 
specialist. The first few chapters deal inar flow configurations are derived. A 
with flow problems including the impor- short chapter on turbulence is presented, 
tant one of motion over a rotating Earth. which provides a satisfactory introduction 


Vorticity and viscosity are then discussed, to the subject. In the chapter on tur- 
leading into boundary layers. The math- bulence, the Reynolds stresses are derived, 
ematical treatment is kept quite simple, the Prandtl mixing length theory is dis- 


although vector notation and differential cussed, and definitions of the quantities 
operator notation is used throughout. most commonly used in a statistical theory 
The author then discusses turbulent flow of turbulence are defined. The chapter 
in a very simple manner. This is followed concludes with a short discussion of the 
by other important atmospheric effects, principles of hot wire anemometry along 
such as buoyant convection, plumes and with some not too recent turbulence 
jets, airwaves produced by explosions and measurements. 
other sources. There is a short discussion A rather complete chapter on turbulent 
of topical interest on clouds and fallout. flow in tubes and annuli is presented, in- 
Throughout the book there are very clear cluding numerous experimental velocity 
diagrams and excellent photographs. I distribution and friction factor data. 
think, really, that even if one is not a The authors are to be commended for dis- 
professional meteorologist one might want cussing in considerable detail the concept 
to read this book because it explains so of the laminar sublayer in turbulent flow. 
many things: The shape of plumes from The notion of the laminar sublayer is us- 
smokestacks, the shape of clouds, why ually given a superficial treatment, and it 
heavy cold air doesn’t necessarily fill up a is rewarding to have this elusive concept 
depression, and others. placed in a reasonable perspective. After 
From a technical point of view this discussing the inherent difficulties involved 
book should be valuable to the engineer in experimentally detecting a laminar 
who is concerned with the problem of ex- sublayer and citing the most significant 
plosion waves, the dispersal of toxic gases in attempts at such a detection, the authors 
an atmosphere, rocket exhaust dispersion suggest that the laminar sublayer may 


and air pollution. best be considered as a statistical concept. 


Part Two concludes with discussions on 
flow near the entrance of tubes, external 
flows over flat plates, cylinders and bodies 
of revolution, and a chapter on the shell- 
side flow in heat exchangers. 

Part Three of the book parallels rather 
closely Part Two in the arrangement of 
material. The various flow configurations 
considered in Part Two as isothermal flows 
are treated in Part Three with heat trans- 
fer present. Dimensional analysis is :ip- 
plied to the energy equation to produce the 
usual dimensionless groups arising in 
forced convection heat transfer. A very 
comprehensive discussion of the analogy 
between momentum transfer and heat 
transfer is presented. Included also is a 
chapter on heat transfer with liquid metils, 
a very pertinent topic for the modern 1u- 
clear age engineer. 

The book contains numerous useful ref- 
erences and also contains many appropri- 
ately chosen example problems throughoitt. 
A striking characteristic of the book is the 
good balance between handbook data and 
theory. Because of the degree of com- 
pleteness possessed by this book, it will 
be a very profitable addition to the library 
of practicing engineers and teachers, as 
well as students seeking insight into con- 
vective heat transfer phenomena. 


Introduction to Heat Transfer, third ed., 
by Aubrey I. Brown and Salvatore M. 
Marco, McGraw Hill Book Co., Inc. 
New York, 1958, 332 + xvi pp. $6.75. 

Reviewed by R. L. Youne 
University of Tennessee 

This is the latest edition of a text which 
has been widely used for a first course in 
heat transfer for third and fourth year 
undergraduate engineers. Major changes 
in this edition are the addition of a chapter 
on fluid flow, inclusion of material on 
numerical solutions of steady and transi- 
ent conduction problems, modernization 
of the information on forced convection 
and more complete consideration of units 
and dimensions. 

The chapter ‘Applications of the Prin- 
ciples of Heat Transfer to Design Prob- 
lems” has been expanded. As in earlier 
editions, this section makes the beginning 
student well aware of the engineering 
judgment that must be used in the design 
of heat transfer apparatus. A generous 
selection of values for equilibrium and 
transport properties is given in convenient 
graphical and tabular forms. 

A few years ago, the reviewer selected 
this book as the text for an undergraduate 
course. At that time, certain aspects of 
the book caused trouble and the reviewer 
is surprised that this revised edition re- 
tains some of these bothersome features. 
Among them are: 


1 In discussing transient conduction, 
the expressions “heat content’’ and “heat 
stored’? are used. Admitting that the 
usage of “heat’’ in heat transfer termi- 
nology is not in accord with that of thermo- 
dynamics, still, terms of this sort should 
be avoided; particularly in a text designed 
to closely follow the student’s first exposure 
to thermodynamics. 


2 In radiation, the authors retain the 
F, (configuration factor) Fe (emissivity 
factor) concept. This is misleading, for 
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This is the cabin interior of Lockheed’s test Electra airplane. The 
engineer adjusts control panel for oscillographic recorder at lower right 


"WE USE DUPONT LINO-WRIT 4 = 
FOR OUR MOST CRITICAL TESTS” oe: 


—Dynamics Flight Test Laboratory, 
Lockheed Aircraft Corporation, Burbank, California 
_ Here’s a report from a Lockheed Group Engineer. heed: it can be used as an intermediate in making — d 
: “The photorecording papers we use have to deliver diazo reproductions; it is strong and durable, capable — ; 
dependable performance while combining a high of withstanding stresses incurred in processing; no— 
enough degree of emulsion sensitivity to handle quality problems with either exposed or unexposed 
special tests. We started using Lino- Writ 4 early in rolls because of high temperature and humidity — 
1957 and now use about 2000 feet per week. storage conditions—even the oldest records are still 
“The first advantage we noticed was Lino-Writ’s clear and thoroughly legible. 
( extra thinness which allows us to load almost 100% Isn’t it time you investigated the Lino-Writ line 
more footage in a 12-inch magazine than we could for your oscillographic installation ? 
- before. This meant that we could run longer tests E. I. du Pont de Nemours & Co. (Inc.), Photo ; 
Canada: Du Pont of Canada Limited, Toronto. 
re Other advantages of Lino-Writ 4 cited by Lock- tpl 
= This advertisement was prepared 
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Better Things for Better Living . . . through Chemistry 
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SUBSIDIARY OF CONTINENTAL MOTORS CORPORATION 


-Development 


«Production 


RAMJET 


@ Continental Aviation & Engineering Corp. is exceptionally 
well qualified, both by experience and by facilities, for work 
on the weapons systems of tomorrow. Our background 
embraces not only a half-century of internal combustion 
engine experience, but also years of pioneering in gas turbine 
engine development, and a full decade of work in the field 
of solid fuels for ramjet propulsion of missiles and target 
drones . . . Continental is staffed and equipped for a wide 
range of assignments, military and commercial. The Detroit 
Division Research and Development Department is supported 
by our modern-to-the-minute Component Testing Laboratory 
complete with environmental facilities located at Toledo. The 
Toledo Production Division now producing various turbine 
engines in volume is capable of supporting diversified pro- 
grams... The CAE record of achievement is one of which many 
a larger company might be proud. Inquiries are invited from 
those having propulsion problems, on the ground, on the 


aad 


CONTINENTAL AVIATION & ENGINEERING CORPORATION 


12700 KERCHEVAL AVENUE, DETROIT 15, MICHIGAN 


water, in the air. ; 


RECIPROCATING 


newer considerations clearly show that 
configuration and emissivity effects are not 
separable. 


3 As the book is written for a student 
of very limited mathematical background, 
adequate treatment of transient conduc- 
tion is difficult. The authors present the 
differential equation of transient conduc- 
tion without derivation and without men- 
tion of the fact that the equation assumes 
no internal heat generation and constint 
properties. The reviewer believes tliat 
the student of such mathematical back- 
ground would benefit from a careful deri- 
vation of this equation. Failure to note 
restrictions on this and other equations 
might be defended on the grounds that it 
leads to less confusion on the part of the 
beginning student, but such omissions 
may be serious as the sudent uses the book 
for reference purposes. 

Forms of the solutions for transient 
conduction in simple situations are ob- 
tained by dimensional analysis. Admit- 
tedly, this additional acquaintance with 
the method of dimensional analysis is bene- 
ficial, but the reviewer doubts the wisdom 
of employing this technique in a situation 
where it is so clearly not needed. How- 
ever, the procedure of using the resulting 
expressions and graphs is made quite 
simple. 

4 As was the case in earlier editions— 
the reviewer believes that Figure 14-la 
illustrating a map of conduction through 
an object of irregular shape is incorrect. 
Surely it must confuse the student to 
read that ‘‘constant-heat-flow lines inter- 
sect the isothermals at right angles,’’ and 
then observe that they do not in Figure 
14-la. Once again it is not clearly 
pointed out that all of this is based on 
constant thermal conductivity. 

In reviewing this book, one must keep in 
mind its intended use. Quoting from the 
preface: ‘The authors’ purpose. . .has 
been to present the fundamentals of heat 
transfer in a manner readily understand- 
able to engineering students and engineers 
in practice who have completed the courses 
in physics and mathematics usually re- 
quired in the first two years of engineering 
curricula and an elementary course in heat 
power.” 


Having determined to write a book for 
students of such limited background in 
mathematics and transport phenomena, 
the authors’ approach has had to be 
largely empirical. Thus, mathematics 
beyond simple integral calculus has been 
avoided (perhaps too much so in light of 
recent advances in second year mathe- 
matics) and the treatment of convection 
is primarily based upon dimensional anal- 
ysis-experimental considerations. 

Though the reasons for the empirical 
approach are clear, one may question the 
wisdom of giving a heat transfer course to 
students of this background. Perhaps at 
this level their time might be better spent 
in additional mathematics or in courses 
leading to a better understanding of the 
fundamentals of mass, momentum and 
heat transfer. But such questions are in 
the province of curriculum committees, 
and should they decide that a heat trans- 
fer course is needed at this level, this book 
deserves careful consideration as a possible 
text. 
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he- | ,.. Requires no refrigeration ... Easier to store than liquid oxygen N,0, Oxidant for liquid 
Performance of Nitrogen Tetroxide as an oxidizer for many fuels is comparable rocket propellants 7 
to that of liquid oxygen. In convenience it is far superior to liquid oxygen. Molecular weight 92.02 a 
ical N2O4 can be shipped easily and stored indefinitely without refrigeration in Boiling Point 21°C | 

the | Ordinary carbon steel containers. It is a dense, mobile liquid that is noncorrosive pr: kage —11.3°C 

.to | if kept dry. The quantity of oxygen contained per unit volume of N20Qx is 1.01 Vaporization 99 cal/gm @ 21°C 
‘at | Kg/liter at 20°C. Critical Temp. 158°C 


ent N20x is available in tank car quantities from Allied’s Hopewell, Virginia plant. gp ecitic Heat of Liquid 0.36 cal/gm 
ses For experimental purposes, N20. is —10 to 20°C 
available in 125 lb. cylinders and 2000 lb. 1.45 at 20°C 
lied containers. Write for technical data, and atl atm 
ie 


information on prices and delivery. Vapor Pressure 2 atm at 35°C 


NITROGEN DIVISION 
Dept. NT 4A-1: 4A-12-4, 40 Rector St., New as N.Y. 
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= 
DECONTAMINATION 
BOOTH 


Instantly ready for 
vital FIRST AID 


Miscues and accidental exposure to dan- 
gerous propellants and other chemicals can 
occur with shocking suddenness. Adequate 
water irrigation is an important key to min- 
imizing such injuries and subsequent 
claims. HAWS Decontamination Booth is 


the answer...a complete safety station for 

immediate first aid. 

HAWS MODEL 8600 2 | 
DECONTAMINATION BOOTH 
is made of durable, lightweight reinforced 
fiberglass plastic, and features Haws Eye- 
Face Wash Fountain, eight lateral body sprays 
and overhead spray unit. All are simultane- 
ously activated by weight on the base- 
mounted foot treadie! Contaminated victims 
are instantly “covered” with water that 
floats away foreign matter from body and 
clothing. 
At aeronautical and astronautical installa- 
tions everywhere, HAWS Safety Facilities are 
important in boosting air-age safety pro- 
grams. Find out what this equipment can 
mean to your operation. Full details sent on 
request, with no obligation. 


DRINKING FAUCET CO 


1443 FOURTH ST. © BERKELEY 10 « 


CALIFORN 


EXPORT DEPARTMENT * 19 COLUMBUS STREET 
SAN FRANCISCO 11, CALIFORNIA 
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| Calorimeter, by 


Technical Literature 


| Instrumentation, Tele- 


metering, Data Recording 


Telemetry: Key to Missile Operation, 
by Raymond M. Nolan, Missiles and 
Rockets, vol. 3, June 1958, pp. 67-71. 

An Adiabatic Calorimeter for the Range 
30° to 500° C, by E. D. West and D. C. 
Ginnings, J. Res., Nat. Bur. Standards, vol. 
60, April 1958, pp. 309-316. 

Satellite ‘‘“Eyes’’ to View Earth Weather 
Conditions, by William Strouse, Missiles 
and Rockets, vol. 3, June 1958, pp. 72-75. 

Cameras: Top Testing Tool for Mis- 
siles, by Peer Fossen, Missiles and Rockets, 
vol. 3, June 1958, pp. 134-136. 

Automatic Calibrator Speeds Missiles, 
by Donald E. Perry, Missiles and Rockets, 
vol. 3, June 1958, pp. 144-145. 

Electroanalogic Methods (Part VI), by 
Thomas J. Higgins, Appl. Mech. Rev., vol. 
11, May 1958, pp. 203-206. 

Some High-pressure, High-temperature 
Apparatus Design Considerations: Equip- 
ment for Use at 100,000 Atmospheres and 
3000° C, by H. Tracy Hall, Rev. Sci. Instr., 
vol. 29, April 1958, pp. 267-275. 

Cryoscopic Assembly for Precise Meas- 
urements under Controlled Atmospheres 
at Temperatures up to 500°C, by C. 
Solomons and G. J. Janz, Rev. Sci. Instr., 
vol. 29, April 1958, pp. 302-304. 

Instrumentation Plays Vital Role in 
Rocket Development, by Kurt R. Stehling, 


Aviation Age, vol. 30, July 1958, pp. 20-21, 
| 76-81. 


A Fluidless Direct-drive Moving Bomb 
W. A. Keith and H. 
Mackle, Trans. Faraday Soc., vol. 54, 
March 1958, pp. 353-366. 
—— Effects in Interferometry, by 
. Kahl and D. B. Sleator, Aberdeen 
Proving Ground, Ball. Res. Labs., Rep. 
1040, Nov. 1957, 21 pp. 


Guidance Systems and 
Components 


Satellite Guidance, Antennas, Doppler: 
Highlights of Dayton Confab, by James 
Holihan, Aviation Age, vol. 30, July 1958, 
p. 118. 

Mace’s Atran Guidance Resists Jam- 
ming, by Philip J. Klass, Aviation Week, 
vol. 68, June 23, 1958, p. 53. 

Environment Burdens Satellite Control, 
Aviation Week, vol. 68, June 30, 1958, pp. 


| 37-39. 


Project Vanguard Report No. 21: Mini- 


| track Report No. 2: The Mark II Mini- 
| track System, by Roger L. 


Easton, J. 
Brit. Interplanet. Soc., vol. 16, May-June 
1958, pp. 390-419. 

Guided Flight Trials, by R. W. M. Bos- 
well, Roy. Aeron. Soc., vol. 62, June 
1958, p. 408. 


Eprtor’s Note: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged. 


Digest 


_ M.H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


Servomechanisms and 
Controls 


On the Theory of Hydraulic Servomo- 
tors without Feedback, by K. C. Kolesni- 
kov, Akademiia Nauk SSSR, Izvestia, 
Otdelenie Tekhnicheskikh Nauk, no. 2, 
Feb. 1958, pp. 126-128. 


Flight Vehicle Design and 
Testing 


Developmental State in 1957 of Un- 
manned Airplanes, Supersonic Rockets 
and Space Vehicles, by Eugen Singer, 
Forschungsinstitut fur Physik der Strahlai- 
triebe, Mitt. 12, Oct. 1957, 143 pp. (in 
German). 

Monitor Designed to Air Jet Takeoffs, 
by Philip J. Klass, Aviation Week, vol. 68, 
June 23, 1958, p. 65, 69-71. 

Swiss Build Mobile Anti-aircraft Mis- 
sile, by David A. Anderton, Aviation 
Week, vol. 68, June 30, 1958, p. 40-45. 

Beyond the Air Age, by John B. Me- 
daris, J. Franklin Inst., vol. 265, May 
1958, pp. 363-370. 

Generalized Variational Approach to the 
Optimum Thrust Programming for the 
Vertical Flight of a Rocket, by Angelo 
Miele and Carlos R. Cavoti, Leitschrift 
fiir Flugwissenschaften, vol. 6, April 1958, 
pp. 102-109 (in English). 

Optimum Thrust Programming along 
Arbitrarily Inclined Rectangular Paths, by 
Angelo Miele and Carlos R. Cavoti, 
Purdue Univ., School of Aeron. Engng., 
Rep. A-57-8 (AFOSR-TN 58-48; ASTI A 
AD 148088), Dec. 1957, 19 pp., 6 figs. 

The Long-range Glide Rocket, by T. R. 
F. Nonweiler, Aeronautics, vol. 38, May 
1958, pp. 36-38. 

Two Approaches Used in First Produc- 
tion Nose Cones, by Michael Yaffee, 
Aviation Week, vol. 68, May 12, 1958, pp. 
52-55, 59, 61, 63, 65, 67, 69. 

The Problem of Maximum Range in a 
Uniform Gravitational Field, by B. Fraejis 
de Veubeke, Astronautica Acta, vol. 4, 
no. 1, 1958, pp. 1-14 (in French). 

Interplanetary Ballistic Missiles—a 
New Astrophysical Research Tool, by S. F. 
Singer, Astronautica Acta, vol. 4, no. 1, 
1958, pp. 59-69. 

Work of the German Rocket Collectives 
in the Soviet Union, by Helmut Gréttrup, 
Raketentechnik und Raumfahrtforschung, 
vol. 2, no. 2, April 1958, pp. 58-62 (in 
German) 

Missile Development in a Hurry: 
Navy’s Polaris IRBM, by Victor DeBiasi, 
ye Age, vol. 29, June 1958, pp. 20- 
21, 79-89. 
Step Vehicles Approach 
Ideal Rocket Performance, by Hugh R. 
Wahlin, Aviation Age, vol. 29, June, 1958, 
pp. 166-172. 

Vanguard Components to Get Wide Use, 
by Evert ¢ Clark, Aviation Week, including 
Space Technology, vol. 68, June 16, 1958, 
pp. 131-136. 

Vital Role Set for New Missile Range, 
by Russell Hawkes, Aviation Week, in- 
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cluding Space Technology, vol. 58, June 
16, 1958, pp. 197-209. 

Miniature Missile Range: Testing 
Models at 20,000 g’s—2,000 mph, by 
D. L. Duff, Canadian Aviation, vol. 31, 
June 1958, pp. 42-43. 


Space Flight ev 
An Approximate Analytical Method for 
Studying Entry into ne Atmos- 
pheres, by Dean R. Chapman, NACA TN 
4276, May 1958, 101 pp., diagrs., tab. 

Human Factors in Space Flight, by 
Eugene B. Konecci, Aero-Space Engng. 
(formerly Aeron. Engng. Rev.), vol. 17, 
June 1958, pp. 34-40, 48. 

Exit and Re-entry Problems, by G. V. 
Bull, K. R. Enkeuhus and G. H. Tidy, 
Aero-Space Engng. (formerly Aeron. Eng- 
ny. Rev.), vol. 17, June 1958, pp. 56-62. 

New Engineering Regime: Orbit Me- 
chanics, by K. Evan Gray, Aviation Age, 
vol. 29, June 1958, pp. 174-181. 

Research Reveals New Problems in 
Space, by James A. Fusca, Aviation Week, 
including Space Technology, vol. 68, June 
1958, p. 76. 

NACA Shifting to Space Agency Role, 
by Robert Hotz, Aviation Week, including 
Smee Technology, vol. 68, June 16, 1958, 

. 79-82. 
ARPA Shapes Military Space Research, 
Aviation Week, including Space Technol- 
oqy, Vol. 68, June 16, 1958, pp. 83-85. 

USAF Space Effort Based on Research, 
Aviation Week, including Space Technol- 
ogy, Vol. 68, June 16, 1958, pp. 86-88. 

Navy Moves to Establish Space Mis- 
sion, Aviation Week, including Space 
Technology, vol. 68, June 16, 1958, pp. 
89-90. 

Army’s Mission in Space is Expanding, 
Aviation Week, including Space Technol- 
ogy, Vol. 68, June 16, 1958, pp. 91-92. 

Industry Faces Stiff Space Age Rivalry» 
Aviation Week, including Space Technol- 
ogy, vol. 68, June 16, 1958, pp. 93-94. 

Man-in-Space Timetable Still Debated, 
by Irving Stone, Aviation Week, ine luding 
Space Technology, vol. 68, June 16, 1958, 
pp. 97-103. 

Medicine Paces Man’s Flight in Space, 
by Craig Lewis, Aviation Week, including 
Space Technology, vol. 68, June 16, 1958, 
pp. 105-115. 

X-15 Dyna-Soar Will Put Man into 
Space, by Richard Sweeney, Aviation 
Week, including Space Technology, vol. 
68, June 16, 1958, pp. 117-126. 

Army Gaining Vital Space Assignments, 
Aviation Week, including Space Technol- 
ogy, vol. 68, June 16, 1958, pp. 137-144. 

IGY Research into Space Continues, by 
Michael Yaffee, Aviation Week, including 
Space Technology, vol. 68, June 16, 1958, 
pp. 153-159. 

New Test Tools Needed for Space Age, 
Aviation Week, including Space Tech- 
nology, vol. 68, June 16, 1958, pp. 179- 
185. 

Canaveral Supports Space Exploration, 
Aviation Week, including Space Technology, 
vol. 68, June 16, 1958, pp. 187-195. 

Space Navigation Challenges Designers, 
by Philip J. Klass, Aviation Week, includ- 
ing Space Technology, vol. 68, June 16, 
1958, pp. 217-225. 

Soviet Lead in Space Due to Foresight, 
Aviation Age, including Space we 
vol. 68, June 16, 1958, pp. 285-288 

Europe’s Space Research is on Paper, 
by David A. Anderton, Aviation Week 
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Our Advanced Design Section 
offers these career positions 


Engineering Specialists. Equivalent of PhD in physics or electrical 
engineering. Seven years experience in analysis and synthesis of 
dynamic systems and controls or related fields. Previous experience 
with rocket engines desirable but not necessary. 

Senior Engineer. Equivalent of MS in physics or electrical engineer- 
ing and four years experience in analysis and synthesis of dynamic 
systems and contrajs or related fields. Rocket engine experience 
desirable but not necessary. 

Dynamic Analysis and Control. MS or equivalent. In addition to 
being well founded theoretically, he should have at least four years 
of experience in dynamic analysis—preferably applied to liquid rocket 
engine control. He must be able to synthesize pump and pressure fed 
rocket engine systems, performing both linear and non-linear analysis 
of the steady-state and transient system operation. He must also have 
the ability to guide and utilize analytic investigations of system stabil- 
ity, and to interpret the results in terms of practical solutions for new 
engine designs. This job requires both sound engineering judgment 
and the ability to direct the efforts of others. 

Mission Analysis. MS, with PhD desirable. Experience in systems 
employing advanced liquid rocket engines. He must have the ability 
to direct the efforts of others in the analytic solution of space flight 
equations, as well as the experience needed to interpret the effect of 
the solution on the propulsion system. Experience in the use of an 
IBM 704 will be useful. He will investigate the rocket engine require- 
ments of future space missions, trajectory and flight optimization, and 
system attitude control requirements. 

Rocket Engine System Analysis. MS or equivalent—plus at least 
four years of engineering experience in system anaylsis and evalua- 
tion of propulsion systems. He should also have a deep insight into 
the solution of gas dynamic and thermodynamic problems of liquid 
rocket engines. He will direct the translation of mission requirements 
into optimum designs of propulsion systems. 

Analytical Turbopump Specialists. MS or PhD, with experience in 
the analysis of turbomachinery design and integration into high per- 
formance engine systems or feed systems. Emphasis and interest in 
advanced concepts and unconventional design practice are particularly 
desirable. Experience in axial flow machinery is preferred. 
Analytical Thermodynamics and Heat Transfer Specialists. 
MS, PhD, or equivalent experience in muitidimensional analysis of 
complex systems and components of rocket engines. Considerable 
experience with approximate solutions and familiarity with unconven- 
tional materials, properties, and design concepts is desirable. 

Senior Engineers. BS in engineering and experience in preliminary 
design and development of advanced rocket propulsion systems for 
customer liaisons and studies relating to applications of rocket engines. 


Please write to Mr. A. B. Jamieson, Engineering Personnel Dept., 
6633 a Avenue, Canoga Park, California. 


ROCKETDYNE F2 


A DIVISION OF NORTH AMERICAN AVIATION, INC. 


FIRST WITH POWER FOR OUTER SPACE 
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ENGINEERS 
AND 
SCIENTISTS 


Here is your opportunity to grow 
with a young, expanding subsidiary 
of the Ford Motor Company. Out- 
standing career opportunities are 

in Aeronutronic’s new RE- 
SEARCH CENTER, overlooking the 
Pacific at Newport Beach, and the 
facility in Glendale, California. 
You will have all the advantages of a 
stimulating mental environment, 
working with advanced equipment 
in a new facility, located where you 
can enjoy California living at its 
finest. 


PhD and MS RESEARCH SPECIAL- 
ISTS with 5 to 7 years’ experience in heat 
transfer, fluid mechanics, 
an 
thermoelasticity. To work on theoretical 
and experimental programs related to re- 
and advanced rocket 
propulsi assignments are open 
in design, high temperature 
materials studies, boundary layer heat trans- 
fer with chemical reaction, thermal stress 
analysis and high temperature thermo- 
dynamics. 
PROPULSION ENGINEERS with 5 years’ 
experience in liquid and solid rocket design 
and test. Familiarity with heat transfer 
problems in engines desirable. To work on 
program of wide scope in R & D of advanced 
concepts in rocket engine components, and 
for missile project work. 
ADVANCES AERODYNAMIC FACIL- 
ITY DESIGNER. Advanced degree de- 
sired. To supervise work in design and in 
instrumentation of advanced aerodynamic 
test facilities such as shock tubes, shock 
tunnels, plasma-jets, and hyper-velocity 
guns. 
STRUCTURAL ANALYSIS SECTION 
SUPERVISOR with 8 to 10 years’ experi- 
ence, including supervision, in the missile 
field. Graduate degree for design and 
analysis required. Will be req 
apply knowledge of high temperature ma- 
terials and methods, thermal stress, dynam- 
ics, etc. to advanced hypersonic vehicles, 
re-entry bodies, and space vehicles. 
FLIGHT TEST & INSTRUMENTATION 
ENGINEERS with 5 to 10 years’ experi- 
ence in laboratory and flig nt test instru- 
Will d p tech- 
niques utilizing ‘advanced instrumentation 
associated with space vehic! 
THEORETICAL AERODYNAMICIST. 
Advanced degree and at least 5 years’ ex- 
ae in high-speed aerodynamics. 
owledge of viscid and inviscid gas flows 
required. To work on program leading to 
advanced missile configurations. Work in- 
volves analysis of the re-entry of hypersonic 
missiles and space craft for determining 
optimum configuration. 
DYNAMICIST. Advanced degree, ap- 
plied mathematics background, and ex- 
perience in missile stability analysis de- 
sirable. Work involves re-entry dynamics 
of advanced vehicles and dynamic analysis 
of space craft. 
ENGINEER or PHYSICIST. With ex- 
perience in the use of scientific instru- 
ments for making 
Work related to o~% test and facility 
instrumentation. Advanced degree de- 
sired with minimum of 3 years of related 
experience. 


Qualified applicants are invited to send 
résumés and inquiries to Mr. K. A. Dunn, 
Aeronutronic Systems, Inc. 


AERONUTRONIC 


A subsidiary of Ford Motor Company 


1234 Air Way, Bidg. 19, Glendale, Calif. 
CHapman 5-6651 


Newport Beach, Glendale, mg Ana and 
Maywood, Californ 
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including Space Technology, vol. 68, June 
16, 1958, pp. 289-290. 

Torques on a Satellite Vehicle from In- 
ternal Moving Parts, by R. . Roberson, 
J. Appl. Mech., vol. 25, June 1958, pp. 
196-200. 

The Equations of Motion of a Tumbling 
Re-entry Body, by E. F. Dobies, Calif. 
Inst. of Techn., Jet Prop. Lab., Progr. Rep. 
20-339, Nov. 1957, 12 pp. 

Solenoid Satellites, by W. B. Klemperer 
and E. T. Benedikt, Astronaulica Acta, 
vol. 4, no. 1, 1958, pp. 25-30. 

On Relativistic Rocket Mechanics, v4 
J. M. J. Kooy, Astronaulica Acta, vol. 4 
no. 1, 1958, pp. 31-58. 

The Communication Satellite, by R. P. 
Haviland, Astronaulica Acta, vol. 4, no. 
1, 1958, pp. 70-89. 

Instrumentation of Artificial Satellites, 
by F. I. Ordway III, Astronautica Acta. 
vol. 4, no. 1, 1958, pp. 90-110. 

Our Philosophy of Space Missions, by 
Krafft A. Ehricke, Aero-Space Engng., 
(formerly Aeron. Engng. Rev.), vol. 17, 
May 1958, pp. 38-43. 

Scientist Compares US-Red Satellites, 
by S. B. Kramer, Aviation Week, vol. 68, 
May 26, 1958, pp. 50-51, 53, 55, 59. 

Some Results of the Observations of the 
First Russian Satellite, by H. K. Paetzold, 
Raketentechnik und Raum fahriforschung, 
vol. 2, no. 2, April 1958, pp. 50-54 (in 
German). 

_ Explorer I (1958) Satellite, by Dietrich 
E. Kélle, Raketentechnik und Raumfahrt- 


forschung, vol. 2, no. 2, April 1958, pp. 


62-64 (in German). 

Preliminary Orbit Information for US- 
SR Satellites 1957 a 1 and a 2, by G. F 
Schilling and T. E. Stern, Smithsonian 
Contrib. to Astrophys., vol. 2, no. 10, 1958, 
pp. 191-198. 

Additional Orbit Information for USSR 
Satellites 1957 a 1 and 81, by J. S. Rine- 
hart and G. F. Schilling, Smithsonian 
Contrib. to Astrophys., vol. 2, no. 10, 1958, 
pp. 199-205. 

Glossary of Astronomical Terms for the 
Description of Satellite Orbits, by J. 
Ashbrook, G. F. Schilling and T. E. Sterne, 
Smithsonian Contrib. to Astrophysics, vol. 
2, no. 10, 1958, pp. 211-217. 

Soviet Orbit Predictions and Orbital 
Information for Satellites 1957 al, a2, 
and Beta, by G. F. Schilling and FE. S 
Fergusson, Smithsonian Contrib. to Astro- 
phys., vol. 2, no. 10, 1958, pp. 219-244. 

Visual Observations of Satellite 1957 al 
Made by Moonwatch Stations, by L. 
Campbell Jr. and J. A. Hynek, Smithsonian 
Contrib. to Astrophys., vol. 2, no. 10, 1958, 
pp. 245-274. 

Soviet Orbit Information for Satellites 
1957 a2 and £1, by G. F. Schilling, Smith- 
sonian Contrib. to Astrophys., vol. 2, no. 
10, 1958, pp. 281-284. 

Basic Orbital Data for Satellite 1957 61, 
by L. G. Jacchia, Smithsonian Contrib. to 

Astrophys., vol. 2, no. 10, 1958, pp. 285- 
286. 

Processed Observational! Data for Satel- 
lites 1957 Alpha and 1957 Beta, by R. M. 
Adams, N. McCumber and M. Brinkman, 
Smithsonian Contrib. to Astrophys., vol. 2, 
no. 10, 1958, pp. 287-337. 

Successive Revisions of Orbital Ele- 
ments for Satellite 1957 Beta, by L. G. 
Jacchia, Smithsonian Contrib. to Astro- 
phys., vol. 2, no. 10, 1958, pp. 339-341. 

A Chart for Finding a Satellite’s Dis- 
tance and Elevation, by J. W. Slowey, 
Smithsonian Contrib. to Astrophys., vol. 2. 
no. 10, 1958, pp. 343-344. 

The etna of Space Vehicles, | vy 


C. A. Cross, Aeronautics, vol. 38, May 
1958, pp. 39-41. 

Interplanetary Orbits, by M. Vertregt, 
Brit. Interplanet. Soc. J., vol. 16, March- 
April 1958, pp. 326-354. 

The Satellite Telescope, by F. A. Smith, 
J. Brit. Interplanet. Soc., vol. 16, March- 
April 1958, pp. 361-367. 

A Simple Method of Plotting the Track 
of an Earth Satellite, by R. A. Anderson 
and C. 8. L. Keay, J. Brit. Interplanet. Sov., 
vol. 16, March-April 1958, pp. 355-361. 

Perturbation of Elliptic Orbits by At- 
mospheric Contact, by T. R. F. Nonweilcr, 
Brit. Interplanet. Soc. J., vol. 16, Mare|i- 
April 1958, pp. 368-380. 

Spatiography: Geography for Space, | y 
Dr. Hubertus Strughold, Missiles ard 
Rockets, vol. 3, May 1958, pp. 106, 108. 

Red Sputniks for War? by Dr. Raymond 
L. Garthoff, Missiles and Rockets, vol. 3 
May 1958, pp. 134, 136. 

Interplanetary Communication and 
Navigation, by Dr. P. A. Castruccio, We: 
inghouse Engineer, vol. 18, May 1958, p). 
88-92. 

Mars, Moon Bases Foreseen in 20 
Years, by Irving Stone, Aviation Weei., 
vol. 68, June 30, 1958, pp. 20-21. 

Space Communications Techniques 
Ready, by James A. Fusca, Aviation Week:, 
including Space Technology, vol. 68, June 
1958, pp. 229-233. 

Space to Spark Avionics Revolution, 
Philip J. Klass, Aviation Week, including 
Space Technology, vol. 68, June 16, 1958, 
pp. 243-263. 


Astrophysics, Aerophysics 


Some Preliminary Values of Upper At- 
mosphere Density from Observations of 
Satellites, by T. Sterne and G 
*. Schilling, Smithsonian Contrib. to As- 
town vol. 2, no. 10, 1958, pp. 207-210 

An Interim Model Atmosphere Fitted 
to Preliminary Densities Inferred from 
USSR Satellites, by T. Ek. Sterne, B. M 
Folkart and G. F. Schilling, Smithsonian 
Contrib. to Astrophys., vol. 2, no. 10, 1958, 
pp. 275-279. 

Cosmic-ray Intensity Variations at Sea 
Level During Magnetic-storm Periods, by 
G. Ramaswamy and 8. D. Chatterjee, 
Canadian J. Phys., vol. 36, May 1958, pp. 
635-637. 

Apparent Temperatures of Some Ter- 
resterial Materials and the Sun at 4.3- 
millimeter Wavelengths, by A. W. Strai- 
ton, C. W. Tolbert and C. O. Britt., J. 
Appl. Phys., vol. 29, May 1958, p. 776. 

Probes Will Explore Cislunar Space, by 
Irving Stone, Aviation Week, including 
Space Technology, vol. 68, June 16, 1958, 
pp. 165-173. 

On Two Parameters Used in the Physi- 
cal Theory of Meteors, by L. G. Jacchia, 
Smithsonian Contrib. to Astrophys., vol. 2, 
no. 9, 1958, pp. 181-187. 

Measurements of Fluctuations of At- 
mospheric Temperature, by A. W. Adey 
and W. J. Heikkila, Canadian J. Phys., 
vol. 36, June 1958, p. 802. 

Transient Radio-frequency Ground 
Waves over the Surface of a Finitely Con- 
ducting Plane Earth, by J. R. Johler, J. 
Res., Nat. Bur. Standards, vol. 60, April 
1958, pp. 281-286. 


Nuclear Propulsion 


On the Mechanics of Photon Jet Pro- 
pulsion, by Eugen Sanger, Forschungsin- 
stitut fiir der Strahlantriebe, 
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Mitteilungen no. 5, Jan. 1956, 92 pp. (in 
German). 

On the Radiation Physics of Photon 
Jet Propulsion and Radiation Weapons, 
ly Eugen Sanger, Forschungsinstitut fiir 
Physik der Strahlantriebe, Mitteilungen 10, 
June 1957, 173 pp. (in German). 

Elementary Approximaticns the 
Theory of Neutron Diffusion, by P. R. 
Wallace and J. LeCaine, Atomic Energy of 
Canada, Lid, Chalk River Project, no. 
336, Aug. 1943, 172 pp. (Reprinted, 
1956). 

Calculation of Internal Pressures in the 
Fuel Tube of a Nuclear Reactor, by B. M. 
Rosenbaum and G. Allen, NACA Res. 
Mem. E52B28, July 1952, 32 pp. (De- 
classified from Secret by authority of 
NACA Res. Abstr. 127, p. 24, 6/5/58.) 

Effect of Capture on the Slowing-down 
Length of Neutrons in Hydrogenous Mix- 
tures Containing Uranium, by H. C. 
Volkin and L. Soffer, NACA Res. Mem. 
1353B05, April 1953, 4 pp. (Declassified 
from Secret by authority of NACA Res. 

ibstr. 127, p. 24, 6/5/58.) 

Control Rod for Temperature Patterns 
Simulating Two Reactor Operating Con- 
ditions, by Tibor F. Nagey, NACA Res. 
Mem. E53B26, April 1953, 18 pp. (De- 
classified from Secret by authority of 
VACA Res. Abstr. 127, p. 24, 6/5/58.) 

Effect of Drawbar Upstream Location 
on Air Velocity Distribution at the Inlet 
Face of Reactor Segment Designed by the 
General Electric Company, by T. F. 
Nagey and E. W. Sams, NACA Res. 
Mem. E52L22, Jan. 1953, 11 pp. (De- 
classified from Secret by authority of 
NACA Res. Abstr. 127, p. 24, 6/5/58.) 

Proceedings, Fluid Fuels Development 
Conference (2nd), Oak Ridge Natl. Lab., 
April 17-18 1952, Atomic Energy Comm., 
CF-52-4-197 (Rev.) 434 pp. 

‘‘Hydrox’’ Cells and Nuclear Energy for 
Spacecraft APUs, by Charles F. Drexel, 
Aviation Age, vol. 30, July 1958, p. 140. 

On the Statistics of Plasmas (The Dy- 
namic Basic Equations for a Classical 
Statistics of Plasmas), by P. 8. Putter and 
F. Sauter, Annalen der Physik, series 7, 
vol. 1, no. = 1958, pp. 4-15. 

Use of Computers in Reactor Design, by 
D. S. Billingsley, W. S. McLaughlin Jr., 
N. E. Welch and C. D. Holland, Ind. and 
Engng. Chem., vol. 50, May 1958, pp. 
741-752. 

A Directory to Nuclear Data Tabulations, 
by R. C. Gibbs and Katherine Way, 
Atomic Energy Comm., Jan. 1958, 185 pp. 

Nuclear Problems of Non-Aqueous 
Fluid Fuel Reactors, by Clark Goodman, 
John L. Greenstadt, Robert M. Kiehn, 
Abraham Klein, Mark M. Mills and Nun- 
zio Tralli, Atomic Energy Comm., MIT- 
5000, Oct. 1952, 265 pp. 

The Fast Critical Assembly, by B. C. 
Cerutti, David Okrent, R. E. Rice, F. W. 
Thalgott and H. V. Lichtenberger, Alomic 
Energy Comm., ANL 5513, Jan. 1956, 28 
pp. 

The Fast Exponential Experiment, by 
F. H. Martens, Atomic Energy Comm., 
ANL-5379, Nov. 1955, 116 pp. 

Neutron Energy Distribution Inside the 
Fast Reactor, by Norris Nereson, Atomic 
Energy Comm., LA 1192, Dec. 1950, 49 pp. 

Description of Developmental Fast 
Neutron Breeder Power Reactor Plant, 
Atomic Power Dev. Assoc., Inc., Alomic 
Energy Comm., AP DA 108, Sept. 1955, 
104 pp. 

Estimation of Fission Product Spectra 
in Discharged Fuel from Fast Reactors, 
by Leslie Burris Jr. and Ira G. Dillon, 
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Plasma Propulsion 


at Republic Aviation 
RR 


Space-Time Trace: With space as ordinate and time as abscissa, photograph shows develop- 
ment of pinch effect in plasma, followed by shock waves. Picture was obtained with special 
streak camera — part of the instrumentation devised for Republic’s experimental Plasma 


Propulsion program. Each space at top measures an interval of 10 microseconds. 


An experimental Plasma Propulsion System under 
test at Republic Aviation gives promise of a power 
plant ideally suited to space vehicles. The system 
generates plasma from a heavy gas and subjects it 
to magnetic acceleration to produce thrust at high exhaust velocity. 

Research and Development in Plasma Propulsion and in a 
number of branches of Hydromagnetics and Plasma Physics is 
being sharply expanded as part of Republic’s new $35,000,000 Re- 
search and Development Program. Investigations currently in prog- 
ress include studies of plasma generation of electricity and the 
application of Hydromagnetics to Hypersonics. 


Opportunities to Lead Theoretical and Experimental Research 


The Scientific Research Staff welcomes the affiliation of scien- 
tists and engineers of stature in the following fields: 


HYDROMAGNETICS PLASMA PHYSICS 
GASEOUS ELECTRONICS COMBUSTION AND 
DETONATION 
HYPERSONICS AND 
INSTRUMENTATION 


SHOCK PHENOMENA 
HIGH POWER PULSE 


PHYSICAL CHEMISTRY ELECTRONICS 


Salaries commensurate with the high degree of talent and creativity 
required. You work with stimulating associates in a laboratory 
atmosphere. $14,000,000 of additional facilities now being built for 
Republic’s new Research Center in suburban Long Island. 


i Write in confidence directly to 
, Dr. THEODORE THEODORSEN, Director of Scientific Research 
STEP 


FARMINGDALE, LONG ISLAND, NEW YORK 
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< POLARIS A-2 
NIKE ZEUS 


IND MANY OTHER CLASSIFIED PROJE 


organization geared to handle 


Send flight tests in the shortest possible time. 


Call or write 


CELCO DEVELOPMENTS 


MILL STREET «© PHONE 101 
SILVER CREEK, NEW YORK 
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Atomic Energy Comm., ANL-5742, July 
1957, 56 pp. 

The Coupled Aspects of a Fast Thermal 
Critical: ZPR-V, by B. J. Toppel, Atomic 
Energy Comm., ANL 5775, Oct. 1957, 20 
pp- 

Hazards Summary Report on the Oxide 
Critical Experiments, by W. C. Redman, 
J. A. Thie and L. R. Datew, Atomic 
Energy Comm., ANL5715, April 1957, 86 
pp. 

Fission Rate Measurements in_ the 
Experimental Breeder Reactor (CP-4), 
hy Samuel Untermyer, Atomic Energy 
Comm., ANL4811, April 1952, 14 pp. 

A Physical Determination of the Con- 
version Ratio of the Experimental Breeder 
Reactor, by C. D. Curtis, 8. H. Klein, D. 
Okrent, W. C. Redman and 8. Untermyer, 

{tomic Energy Comm., ANL5222, Aug. 
1954, 113 pp. 

History and Status of the EBR, by 
Warner E. Unbehaun, Atomic Energy 
Comm., AECD3712, April 1953, 61 pp. 

Estimation of Fission Product Spectra in 
Fuel Elements Discharged from the Power 
Breeder Reactor and the Experimental 
Breeder Reactor No. 2, by I. G. Dillon and 
Leslie Burris Jr., Atomic Energy Comm., 
Rep. 5334, Oct. 1954, 36 pp. 

Design of Fast Neutron Breeder Test 
Pile, by W. H. Zinn, Atomic Energy 
Comm., TID-10085, Jan. 1946, 30 pp. 

Minimum Critical Mass and Uniform 
Thermal Neutron Core Flux in an Experi- 
mental Reactor, by J. W. Morfitt, Alomic 
Energy Comm., Y-1023, Dec. 1953, 168 pp. 

LAMPRE;; a Molten Plutonium Fueled 
Reactor Concept, by R. M. Kiehn, Atomic 
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FLORIDA 
Openings 
for 
ENGINEERS 


_ Expanding efforts in connec- 
tion with the Boeing Bomarc | 
supersonicair-defense missile 

_ have created a number of — 

rewarding, long-range oppor- 
tunities for engineers. These [ 

are Florida assignments, in a 

_widerange of fields,including: 


Missile Test and Launching 
Modification of Ground 
t Support, Operational Test and 7 
Special Test Equipment 7 Ww 
4 Planning and Scheduling 
of Test Operations 
Evaluation and Coordination 
of Test Operations and 
System Demonstration 
Development of Drawings for — 
Weapon System Test 
Installations 
Planning Engineering 
Requirements for Missile 
Base Implementation 


Requirements are experience 
_ in any of the above fields, and 
a degree in ME,CE, EE, Math 
Physics. 
You will find these Boeing 
assignments interesting and 
challenging. They offer ex- 
cellent opportunities for ad- 
_vancement and long-range 
career growth in a field with 
a truly unlimited future. In 
addition, you'll enjoy the 
many advantages of Florida 
living, with world-famous 
beaches and a host of year- 
round outdoor recreational 
facilities only minutes from 
Boeing’s Florida locations. 


For full details, write Mr. 
Robert H. Glass, Boeing Air- 
Company, Dept. FL-3, 
?. O. Box 97, Melbourne, Fla. 
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ENGINEERS 
Research 


27 
Creative young engineers will find challenging re- 
search projects in experimental hydrodynamics at 
Armour Research Foundation. 
B. S. or M. S. in Mechanical Engineering with a strong 
background in Physics and two years’ experience in 
experimental hydrodynamics required. 
This is a long established research organization 
located in a metropolitan area offering cultural and 
educational advantages. Extensive employee bene- 
fits including a liberal vacation policy. Please send 
resume to: 


LP. Bloch 


RESEARCH 
FOUNDATION 


Illinois Institute of Technology 
10 West 35th Street 


Chicago 16, Illinois 
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SELF-ALIGNING BEARINGS 


PLAIN TYPES ROD END 
TYPES 
PATENTED U.S. A. 


World Rights Reserved 


CHARACTERISTICS 
ANALYSIS RECOMMENDED USE 


: For types operating under high 
Stainless Steel Ball and Race { temperature (800-1200 degrees F.). 
Chrome Alloy Steel Ball For types operating under high radial 
and Race ultimate loads (3000-893,000 Ibs.). 


Bronze Race and Chrome { For types operating under normal loads 
Steel with minimum friction requirements. 


Thousands in use. Backed by years of service life. Wide variety of 
Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in similar 
size range with externally or internally threaded shanks. Our Engi- 
neers welcome an opportunity of studying individual requirements 
and prescribing a type or types which will serve under your demand- 
ing conditions. Southwest can design special types to fit individual 
specifications. As a result of thorough study of different operating 
conditions, various steel alloys have been used to meet specific 
needs. Write for Engineering Manual No. 551.Address Dept. JP-59. 
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